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16380 RELIABILITY AND HYDRAULIC ENGINEERING be 


KEY WORDS: Earthwork; ‘Failure; Hydraulic Hydraulic 
engineering; Hydraulic structures; Hydrology; Levee failures; Levees; <i 
ABSTRACT: | The 1 reliability ‘of two hydraulic en engineering “systems with random 
resistance and/or loading is estimated. The first system is a flood levee reach that is . 
subject to four failure modes: overtopping, boiling, slope sliding and wind wave 

_ erosion. The second system is a flood levee system along a confluence reach. For this 
levee, failure by overtopping depends on the bivariate distribution of water stages in 
the tributary and the main river, linked by a backwater curve. The failure probability, 
which is the complement of reliability, is estimated analytically for the first system, 
and then by hydraulic calculations and simulation for the second system. In each case, | 
estimated reliabilities agree well with experimental observations. “Wert, 
REFERENCE: Duckstein, Lucien (Prot, 1 Dept. of Systems Industrial Engrg., 
Univ. of Arizona, Tucson, Ariz. 85721), and Bogardi, Istvan, “Application of 
Reliability Theory to Hydraulic Engineering Design,” Journal of the Hydraulics 
ASCE, Vol. HY7, Proc. Paper 16588, 1981, 


KEY WORDS: Failure; Flood ‘control; water; Hydraulic structures; 
Hydrology; Mine water; Mining; Probability theory; Reliability; Simulation; ig dldarboey wo. 


ABSTRACT: The reliability of an underground hydraulic engineering system, which is 
characterized by both resistance and loading being random, is estimated. This system 
consists of a mine water control system in which failure may be caused either by 
excessive underground inrushes (or floods) or by equipment breakdown. The failure — 
probability, which is the complement of reliability, is estimated by a combination 1 
failure tree analysis and simulation. The estimated reliabilities nage ro to — well 
REFERENCE: Duckstein, Lucien (Prof., Dept. of Systems and Industrial Engrg., 
Univ. of Arizona, Tucson, Ariz. 85721), Bogardi, Istvan, and Szidarovzky, = 
a: ‘Reliability of Underground Flood Control System,” Journal of the Hydraulics 
_ Division, ASCE, Vol. 107, No. HY7, Proc. Paper 16381, July, 1981, pp. 817-827 
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16372 MODELING GRADUAL DAM BREACHES 
KEY WORDS: Channel erosion; Channels; Dam breaches; Dam stability; 
_ Earth dam performance; Earth dams; Flood control; Mathematical models; _ 
Models; Numerical analysis; Open ‘channel flow; Sediment transport; 
' _ ABSTRACT: A simulation model of the gradual failure of an earth dam has been +. 
"formulated, developed and tested with real-life data. A significant feature of the model 
is its ability to account for the growth of the breach and the eventual draining of the 
_ reservoir. Concepts of water and sediment routing are used in conjunction with a 
channel geometry descriptor to arrive at a self-contained mathematical model of the _ 
~ breach enlargement and the ensuing flood wave. Unsteady flow elements of the 
- simulation are an implicit numerical solution of the complete Saint Venant equations 
- coupled with a sequential sediment routing technique. This approach provides = 
increased rational basis for determining outflow hydrographs from postulated earth 


REFERENCE: Ponce, Victor Miguel (Assoc. Prof. of Civ. Engrg., San Diego State 
‘Univ., San Diego, Calif. 92182), and Tsivoglou, Andrew J., “Modeling Gradual Dam © 
; | 4 Breaches,” Journal of the Hydraulics Division, ASCE, Vol. 107, No. HY7, Proc. Paper 
i. 16372, July, | 1981, PP. 829- 838 
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16364 BOUNDARY ‘SHEAR IN CHANNELS ana 
KEY WORDS: Beds under water; Boundary shear; Channels (oaterweye)s *, 
Depth; Flow measurement; Flumes; Friction; Open channel flow; Rivers; insert 
Roughness (hydraulic); Shear stress 
5 ‘i ABSTR ACT: A series of flume experiments are reported in which the walls ont alleles 
- were differentially roughened, normal depth flow set, and measurements made of the - 
- boundary shear stress distribution. An empirically derived equation is presented giving — 
_ the percentage of the shear force carried by the walls as a function of the breadth/ 
_ depth ratio and the ratio between the Nikuradse equivalent roughness sizes for the bed 
+ and the walls. The results are compared with other available data for the smooth — 
- channel case and some disagreements noted. The systematic reduction in the shear 
force carried by the walls with increasing breadth/depth ratio and bed roughness is 
illustrated. Further equations are presented giving the mean wall and bed shear stress _ 
_ variation with apsect ratio and roughness parameters. Although the experimental data - 
_ is somewhat limited, the equations are novel and thought to indicate the general 
_ behavior of such open channel flows with some success)§ 
REFERENCE: Knight, — Donald W. (Lect., Civ. Engrg. Dept., Univ. of Birmingham, 
P.O. Box 363, Birmingham, B15 2TT, England), “Boundary Shear in Smooth and 
4 _ Rough Channels,” Journal of the Hydraulics Division, ASCE, Vol. 107, No. HY7, 
Proc. Paper 16364, July, 1981, pp. 839-851 


properties; Models: Hydrologic, dat; Statistical 4 
analysis; Stochastic models; Stream flow; Time factors 

ABSTRACT: The identification of streamflow stochastic models is explored. The 
_ identification of the type of stochastic model is made, based on a conceptual physical 
representation of natural watershed. The identification of the form of the model is — 
made, based on the recently developed R-functions and S-functions. ARMA is the © 
autoregressive moving average. For a ARMA(p,q) precipitation input, the ground-water 
storage is an ARMA(p+1, q) process, and the streamflow is an ARMA(p+1, q+1) 
process. In general, such ground-water and streamflow processes belong to the class of 
restricted ARMA processes in the sense that their parameter space is a subspace of © 
that corresponding to the general ARMA models. The form or order of the ground- 
water and streamflow ARMA processes for given historical time series can be uniquely 
identified by using the R-functions and S-functions. An coemgie is given as a 

REFERENCE: Salas, Jose D. (Assoc. Prof., Dept. of Civ. Engrg., Colorado State 

_ Univ., Fort Collins, Colo. 80523), Obeysekera, Jayantha T.B., and Smith, Ricardo A., | 
“Identification of Streamflow Stochastic Models,” Journal of the Hydraulics — 
ASCE, Vol. 107, No. HY7, Proc. Paper 16368, July, 1981, pp. 853- albino aaiaac 


KEY WORDS: Exponential functions; Forecasting; Lakes; Mathematical 


.s models; Observation; Predictions; Statistical analysis; Statistical distributions; 
Time factors; Time Series analysis; Water Water 


management wit. 
ABSTRACT: A Statistical univariate forecasting called Exponentially 


_ Weighted Moving Average (EWMA) is used to obtain the estimates of future water 
: 3 levels of a large lake. The characteristics of the time series data consisting of average x 
monthly lake levels is examined, and the parameters of the EWMA model are 
determined. The ex-post forecasts generated by this model are compared with the 4 
actual observations of lake water levels and with results obtained earlier with other 
stochastic methods. EWMA yields forecasts which are statistically indistinguishable 
REFERENCE: Koppula, Sivajogi D. (Section Head, Design and Construction Div., 
Alberta Environment, 16403-102 Street Edmonton, Alberta, TSX 2G9, Canada), 
“Predicting Lake Levels by Exponential Smoothing,” Journal of the ai 
Division, ASCE, Vol. 107, No. HY7, Proc. Paper 16384, July, 1981, pp. 867-878 
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16365 THREE HYDRAULIC VARIABLES 


KEY WORDS: Alluvial channels; Boundary shear; Erosion; Mathematical 
- models; Open channel flow; Scour; Secondary flow; Sediments; —. Jt 


ABSTRACT: mathematical modeling and computational technique is 
a which is capable of treating and dealing with interactions among hydraulic variables is 
open channel flows, such as the distribution of primary flow velocity, channel cross — 
section, discharge rate, secondary flow, shear stress, and sediment concentration. A. 
- shear stress formula is presented that can include the effect of secondary flow. The 
_ formula gives peak values of boundary shear that occur on the channel bottom = 
| the corners and on the side walls. Such peak values of boundary shear are owing to a : 
momentum transfer to the region by secondary flow. The simulation technique can be 
used as a decision-making tool in selecting engineering designs and controls for 
channel protections. The main role of secondary flow in sediment transport seems to 
E - be to transport sediments in the transverse direction, whereas the gravity brings 
sediments downward and the diffusion brings sediment particles upward. The vertical 
_ component of secondary flow seems to help gravity bring sediment down in most ail 


the area in a transverse cross section of a channel. 


REFERENCE: Chiu, Chao-Lin (Prof. of Civ. Engrg., Univ. of Pittsburgh, Pittsburgh, 
Pa.), and Hsiung, David E., “Secondary Flow, Shear Stress and Sediment Transport,” fey. 
7 Journal of the Hydraulics Division, ASCE, Vol. 107, No. HY7, Proc. Paper 16365, 


16363 FLOW RESISTANCE IN GRAVEL BED RIVERS 

_ KEY WORDS: Boulders; Boundary shear; Channels (waterways); Cross 

sections; Field “tests; Flow resistance; Geometry; New Zealand; Rivers; 

Roughness | (hydraulic); Sediments; Velocity distribution 2 


ABSTRACT: Flow. tine deduced from data on coarse gravel bed 
reaches of New Zealand rivers. Measurements of hydraulic variables were made during 
_ reasonably uniform subcritical flows at cross sections regular in confirmation. Channel 7 
_ bed and bank roughness are not distinguished, and surface bedmaterial is characterized 
_ by the median size. When the channel boundary is rigid the Darcy-Weisbach friction , 
factor is only a statistically significant function of the ratio of hydraulic radius to 7 
median size of surface bedmaterial. A model of this form explained 60 percent of the 
friction factor variance inherent in the field data set. With active bedload transport _ 
conditions and changing bed topography, the Froude Number based on median surface _ 
bedmaterial size was statistically significant, and 40 percent of friction factor variance 


REFERENCE: Griffiths, bs, George A. | A. (Scientist, am and Soil Div, Ministry of W orks , 
and Development, P.O. Box 1479, Christchurch, New Zealand), “Flow Resistance in | 
Coarse Gravel Bed Rivers,” Journal of the Hydraulics Division, ASCE, Vol. 107, No. 
HY7, Proc. Paper 16363, July, 1981, pp. 899-918 


16385 STANDARDS FOR COMPUTER-BASED ) DESIGN STUDIES 

KEY WORDS: Computer applications; Computer models; 

systems hardware; Engineering education; Model studies; Municipal 

engineering; Project planning; Standards; Storm water; management 
ABSTRACT: As a result of the large variety of computer- ee studies available, their 
rapid evolution and their complex structures, these programs are becoming increasingly 
difficult to manage. The following activities could be specified in detail when study 
bids are requested: problem review, requisite accuracy, review of available programs, 
program — selection — criteria, model verification, data preparation and output — 
interpretation, and documentation of the modified program actually used. Many of — 
these activities will be implemented by a reputable engineer in any case, but by 
including these terms of reference, all consultants proposing services will be assured ¥ 
that their estimates will be met, and clients will be ensuring their investment. Each of 
‘these points is explored in general terms; the examples used are appropriate 7 
‘stormwater management modeling studies; however, the principles could be applied to a 
REFERENCE: James, William (Prof., Dept. of Civ. Engrg. and Engrg. Mech., 
-McMaster Univ., 1280 Main Street West, Hamilton, Ontario, Canada, L8S 4L7), and 
Robinson, Mark A., “Standards for Computer-Based Design Studies,” Journal of the 
Hydraulics Division, ASCE, Vol. 107, No. HY7, Proc. Paper 16385, July, 1981, pp. Zz 
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(US. S. Customary-S! CONVERSION ‘Factors 


_ In accordance with the October, 1970 action of the ASCE Board of Direction, which stated 


that all publications of the Society should list all measurements in both U.S. Customary and 


7 


to compute the SI unit values of measurements. A complete guide to the SI system and its 7 


SI (International System) units, the following list contains conversion factors.to enable readers | 


ae All authors of Journal papers are e being asked to prepare their papers in this dual-unit format. 


_ To provide preliminary assistance to authors, the following list of conversion factors and guides 


t tt M t 


yards (yd) meters (m) 
square inches (sq in Square millimeters (mm*) 
square feet (sq ft) square meters (m ) 
‘Square yards (sq yd) "square meters (m 
square miles (sq miles) 
acres (acre) 
- cubic inches (cu in.) 
cubic feet (cu ft) 
cubic yards (cu yd) meters (m 
pounds (Ib) mass kilograms (kg) 
tons (ton)mass kilograms 
kilogram force (kgf) newtons(N) 
pounds per square foot (psf) ( 


acre-feet (acre- -ft) cubic meters (m’*) | 
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APPLICATION OF RELIABILITY THEORY TO HypRAuLic 


_BNGINEERING DESIGN™” 

Bay om By Lucien Duckstein' and Istvan Bogardi6 


= boiling, slope sliding, wind wave erosion) may cause failure. The second 7 
system considered is a confluence reach of a tributary where different flooding . 
_ conditions (flood in the main river or or the ‘tributary, | or both) may result ir in 


During the last two decades, there has been | a rapid development in the 7 


“that structures designed by using deterministic (fixed) values of load and resistance 
= often not be safe enough; in other Pare pition were underdesigned. On 


4 were undertaken to estimate the failure probability. Two approaches to this 
_ problem can be distinguished. The first one is characterized by the work of 
Freudenthal (16,17), who used classical probability computations to Peon all 
the distribution function (DF) or the probability density f function (pdf) of dl 


when load is greater ‘than r resistance ce, ie e., when the safety. factor v < <1, the 

probability of failure P, can be calculated by taking the value of the DF for 

_ which v = 1. Since the structure generally consists of more than one element, 
4 and each element generally has different properties, several modes of failure 

are Possible. The re resistances to mare v various modes are not ot independent 1 pantom 


described (16, 17), has been to simple structures; the s safety 
_of complex structures, or systems, has received relatively little attention (27). 
__The second approach to investigating failures statistically is exemplified by 


2 Prof., Dept. of Systems and Industrial Engrg., Univ. of Arizona, Tucson, Ariz. 85721. 
_ ?Sr. Systems Analyst, Mining Development Inst., 1037 Budapest III, Mikoviny u. 2-4, 


_ Note.—Discussion open until December 1, 1981. Separate discussions should be submitted 7 
for the individual papers in this symposium. To extend the closing date one month, ‘1 
a written request must be filed with the Manager of Technical and Professional Publications, : 
ASCE. Manuscript was submitted for review for possible publication on July 18, 1980. 
This paper is part of the Journal of the Hydraulics Division, Proceedings of the American . 
Society of Civil Engineers, © ASCE, Vol. 107, No. HY7, July, 1981. ISSN 0044- 
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the work of Cornell (10) and Lind (23). These intact 
ree in present design codes but suggest introduction of probabilistic con- ce 

__ siderations into the codes. This approach is of great practical importance in 
‘ as the design of concrete and steel structures where the relatively small randomness ‘ 


of load moderates the effect of the deterministic 


effort to combine both approaches. 
Reviasiuity oF Fiooo Levee Sysrems 
Description. _the probability of failure, P 
; % investigated. The reach refers to a levee line of several kilometers that protects 
ae. given area in the flood plain. The levee reach behaves as a structure. If 
there is a levee rupture at any point of the reach, a structure fails om os. 
whole area under protection becomes inundated. 


f load is a random variable whiee distribution may be estimated using obeervelions 
=m taken at a gaging station along the reach. If there is no effect from 7 ae 
id or wind waves along the reach, _ this load is the Same in every cross ‘Section 


every ¢ cross section a a different pdf of load» should be estimated (4,30). The 
4 Fi pdf of load is estimated by known probabilistic or stochastic methods. Naturally, 
‘ the pdf parameters bear the uncertainty imposed by finite samples (18). 
The approach of Freundethal (16,17) cannot be applied in a simple case like 
ve slope stability analysis of earthworks because the simplest earth structure should _ 
a be regarded as a system subject to other 1 modes of failure (11). The assumptions — : 
~ generally used in earthworks analysis are: (1) The load is assumed to be constant; 
(2) the load and resistance are independent, i.e., statistically rigid strectures: 
; are assumed; (3) relatively simple, often approximate methods of soil mechanics _ 
i are used to estimate resistance; and (4) spatial changes of soil properties are 
q described by linear regression m models. In the case of hydraulic earthworks, 
especially along natural waterways, the ; assumption of constant load cannot 
a In the case of large flat rivers, failure of the levee reach in any cross section ~ 
_ can be caused by one of the following ener (Fig. 1): (1) Overtopping caused 
5 due to adverse soil conditions; (3) failure of the levee slope due to seopage 
and wetting; and (4) wind wave effects (e.g., erosion). The resistance or the | 7 
protection in one cross section should be determined for each of the four effects 7 
i and should be expressed in terms of the flood parameter that triggers the failure - a 
i Generally, there is no single flood parameter governing all four effects: the 


_ maximum flood level, the duration of the flood, the intensity and duration of 
levee soaking, characterized by the value of the so-called flood 
or another parameter may trigger the adverse phenomena. 


_ Resistance values can be determined by direct measurements (e.g I levee crest 


a Omnic a 2 ry > jam he nresen nane enrecen 
} 
| 
| | 


| location and time of adverse Re es and on the corresponding flood 00d load). 
Specifically, re resistance values were estimated in the following way. 
Resistance against - Overtopping, H, .—The triggering flood parameter is the 


peak flood level, and the resistance is the levee height minus 30 cm-80 cm 


from one common n gaging station within the reach have ‘random errors. 
‘Mowers, , the variance e of H, is generally quite small as compared to ‘that of 
_ Resistance against Sand Boiling or Subsoil Failure, H,.—This value can also 
| expressed by the peak flood level or alternatively by the flood exposure | 
or duration. Theoretical methods and empirical equations (34) generally consider 
_ the flood level as the triggering value. However, there have been data showing» 
that given values of flood exposure or duration may cause sand boiling (29). - 4 
7 Since the mechanism of this failure phenomenon is rather complicated and there | 
is no motel acceptable for cases, the resistance H, should 


mag 


we, 


be a random variable with considerable variance (5,6). |The varian I 
ean only be estimated in a systematic way if historical data are available ae 
- the actual flood levels at the time when boiling started; otherwise, only the 
expected value of H,canbeestimated. 
A Bayesian approach (12) offers a way to connect in-situ practical experience © 
- on boiling with regional boiling data. In that case, regional boiling data can 
_ be regarded as prior information (28,35) and observation data of past floods — 
_ can serve for the likelihood function. The Bayes formula yields a posterior 
_~ipdf using the prior information and likelihood function mentioned. Foloyan, 
= Hoeg, and Benjamin (14) used such a Bayesian approach: to determine a posterior 


29410 


and slope sliding may occur. Experience shows that this failure phenomena — 
is influenced not only by the peak flood level but also by the flood duration. | 
_ This complex load can be fairly well characterized by the flood exposure w 
a (3) which is the area of the stage hydrograph (m x day) above bankful stage. 
Classical methods of slope stability analysis cannot regard this flood exposure 
as a load on the levee; therefore, experience from me flood data is used 


HY7 HYDRAULIC ENGINEERING DESIGN™ 
3 
reflecting a spongy layer under the crest. This resistance, H,, is a random ‘ 
: = for at least two reasons: (1) There could be a random fluctuation in “q 
| 
| 


_is a lower threshold water level H, , below which no sliding occurs veguetions = 
of the magnitude of the flood exposure. . This value, H,, may correspond to 
the height of a very old levee core or to the pressure necessary to wet the | 
landslide slope from the bottom up. Theoretically, there exists some function 
_ between the geometry and soil properties of the levee, on the one hand, and 7 
the resistances, W, and H,, on the other hand. The approach applied to determine — 
the DF of resistance H, can be used to estimate the DF of W, and H,, : considered — = 
as random variables. A Bayesian approach provides again an appropriate way 
to estimate W, and H, by combining prior information and observation data. = 
Resistance against Wind Wave Failure—This dynamic failure is triggered 
by the sum of the flood level and wave effect (wave height and run-up) (4). 
In the case of flat, » large rivers, the distance between the levees can be as — 
as km- -8 km. Naturally, adverse wind conditions (large 


_ 20 cm-30 cm. The stochastic model of levee system failure also considers this 
cerns mode (30). However, the example Presented in that paper deals with © 


lues of 
ee cross sections may govern the performance of the whole system ; 
(11). In fact, as cross sections are taken closer and closer, the subreaches 
_ become shorter and the spatial correlation increases. In the example shown, _ 
E bf the distances between c1 cross sections are re large enough so that the spatial correlation a 
between 
been verified by the ‘first order 
I 4 _ between the estimated resistances. No significant correlation has been found — 
according to the Durbin- Watson test (19) on the 5% level of significance. However, 
in other cases, the spatial correlation might take on significant values and woven, 
te neglected. The proposed model can also be used in such cases as it was 
shown by Szidar (33) in a paper examining the effect of random errors 
_ subsequent values of economic loss and cost functions. Between a Pr 
: errors, a strong correlation exists, as it might between neighboring regions in 


7 stochastic character of flood load; Q) the v various ous degrees of protection ~ home 
~ along the levee reach; and (3) the different modes of failure. The model is 
_ presented in two parts: (1) P, in a given cross section; and (2) P, of a levee 
— ”* in Given Cross Section. —There is a failure in the cross section if either: 
(1) A > H,: overtopping; (2) h > H,: subsoil failure; (3) 4 > H,)N (w> W): 
. - slope stability failure; and (4) h + x > X: wind wave attack, erosion in which | 
the random loads are h = peak flood level; w = flood exposure, i.e., the 
area of the stage hydrograph (m x day) above bankful capacity; and x = wave 
height and run-up (5). The random resistances are: es ee 
. H,: the highest flood level allowable with respect to 
bs H,: the highest flood level with respect to boiling. 


, _ resistance against this composite effect can be taken as the levee height minus _ ; 
| 
2 
| 


Hy is ne since a flood exposure Ww, pertain a low flood level 


h < H,, say 10 cm above bankful level for a few weeks, _ which certainly a 
4. X: the highest dynamic water level | eons Static level + wave + run- =) 


(h+x)>X] Uh>H, )U H, )U (A, (w> W)) 
‘The probability. of that is, probability of failure, can be ‘computed 
; a by dividing the right side of Eq. | into disjoint events. As a result, the following 
expression for P, can be reached under the assumptions that x and Gow w) are 
— 


x The relations be transformed into 


in which g= F, x(X — h)f(h, w) dw dh]. The last term 
3 of Eq. 2 accounts for ‘the effect of waves. Thus, in the cases when the distance © . 
_ between the levees is small, say less than | km, considerable waves cannot — 
=a be generated for lack of enough fetch, the first three terms of Eq. 2 are sufficient 7 _ 
to calculate P,. Such will be the case for the practical example presented in — 
- the next section. On the other hand, _ the wave erosion can be a critical factor _ 
‘ : aa Pep of Levee Reach. —For each failure mode, ‘the | levee reach is divided int 


subreaches, within which the mean value and variance of resistance for. eac 
failure mode is constant. Assuming that every subreach is characterized by 
one cross section, the failure event for the whole reach is the union of the 
failure events T,of the subreachesi=1,...,.m: 


H = min (H,,,H H,= min H,,; 1h) min 


in which , and W, are the resistances of subreach i. The system 
Pr, , without considering wind wave attack, can be calculated as: pet ms Yo St 


Note that P,.is considered as a random 1 variable itself because of the uncertainties _ 


It is assumed that the resistances follow normal distributions with means _ 
equal to the deterministic values and W and Standard deviations 


= ts Often there is a spatial dependence between subsequent resistances. This — 


spatial dependence can be characterized by a significant first order correlation | 
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_ , 3. H,, W,: respectively, the smallest necessary flood level and the highest 
| 
a 


Naturally i in cases, +, Ea. 3 not hold and the syst 
- for each failure mode should be represented with an n-dimensional normal > 

Eq. 3 can be solved in! four different onthe various uncertainties 


1. No uncertainty in in ‘the 1 resistances is assumed. mean values of 
are substituted into Eq. 3 which asingle 


2. Uncertainty in the resistances is assumed. Based on estimated means, 
assumed distributions of the ‘Tesistances, a number 


variance system can he also a or continuous ‘distribu- 
tion can be fitted to the ee 
in the estimated stochastic properties of the load are con- 
i sidered: namely, parameters of the joint pdf /(h, w) are estimated from a finite 
_ sample size. From a Bayesian viewpoint, these parameters may be regarded 
« random variables (2); on this basis, equally likely samples and parameters — 
can be generated. Again, with the help of Eq. 3, possible values of P, can — 
be calculated using the generated parameters. Statistical properties of the random 
values of P, can be determined as in the preceding case. 
4. Uncertainties in the resistances and in the load parameters are considered; 
- then one generates simultaneous independent sets of resistances and load 


reach, 14 km in length, of the Sebes-Koros River, Hungary, is sought. A flood 
event occurs each time the stage climbs above 500 cm. For each event, flood 

_ load is given in terms of peak stage and flood exposure data taken in the : 

1901- 1970 at the gaging station, which is located i in the 


‘that his lognormally distributed cannot t be rejected o on n the basis of a a Kolmogorov- 
Smirnov test at the 0.05 level. A similar result is found for w. The sample 
_ mean M and standard deviation D of h and w, respectively, are as follows: | 
@ M(h) = —0.5132; (2) D(h) = 1.1683; (3) M(w) = 0.7197; and (4) D(w) 
= 1.8417. The correlation coefficient between A and w is r= 0.97, which is 
a fairly high value. _ As a result of this high correlation, the difference between 
the P, of the system and that of the weakest section is expected to be small, _ 
since failure caused only by h also causes failure by w, and vice-versa. Note 
_ that this strong linear correlation between h and w cannot be taken for granted 5 
Resistance values H,,, H,,, and W,, are presented in Table 1 and 


= P H, , is a geometric measurement, ‘taken as the height of the levee ‘crest 
_ minus 40 cm; these 40 cm account for a layer of spongy soil just under the 
crest. On the basis of the existing record, the coefficient of weer a V of 


_H, has been taken as 2% for every cross section. 


7 
_ 
| 


TABLE 1.—Resistance Values H,, H,, L, and H, of Subreaches with Associated 
1H 


guilstie without Yo 

limits, in i- jin centi-| in centi- | in 

‘Number kilc meters meters | xX day ; | tainties 

d nad 

phe sues 


0. 0594 


“Coefficient of variation = =2%. 
... P,. for the levee reach without soil-mechanical uncertainties = 0.0643; Mean ; 


i. of P, for the subreaches = 0.0520; E(P,.) for the levee reach with soil-mechanical — 
uncertainties = 0.0865; D(P,.) for the levee reach with soil-mechanical uncertainties . 
0.0161; E(P, -.) for the levee reach with parameter uncertainties = 0.0676; D(P,) for — 
the levee reach with parameter uncertainties = 0.0276; E(P,.) for the levee reach with a 
- soil-mechanical and parameter uncertainties = 0.0887; and D(P,) for the levee reach 
with soil-mechanical and parameter uncertainties = 0.0395. 


af 
| 
4 
} 
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iin resistance against boiling H, is a quantity computed by the formulas _ 
developed for the Mississippi River levees (34), and verified 
during those flood events which actually caused boiling. The coefficient of a 


based on existing data, is taken as as 10% for eve: ry cross s 


variation V of 2» | 


(3. The ‘resistance against sliding W, and the samedi H, are taken as 
the flood exposure and the peak stage of the flood wave observed at the time -_ 
of the critical phenomenon, namely, sliding during the 1970 flood (36). The | 
coefficients of variation V of W, and which reflect the } of soil 
_ parameters, have been taken as 10% for every cross section. ee 
AS shown in Table 1, there are four subreaches, numbered 2, 4, 17, and : 
21, where dangerous boiling has occurred and three subreaches, numbered 8, 
17 and 21, where slope sliding has happened. No entries have been made for 
the values of H, and H, that exceeded the value of H, for any given ‘subreach. 
In this ‘eanngle. the s reach is divided into the same equispeced subreaches for 
—_ failure mode. For the estimation of P, of the system, the value of min 
. rs [H,,H,] was found at subreach 8 (3.5 km-4 km), while W, and H, ae 
‘Nocated at subreach 21 (10 km-10.5 km). Table | also contains the ‘system P,. 


‘than ‘the mean wie of the P,s of the subreaches, found to ) be 0.052, but 
it is also greater than the P, of 0. 0622 for the “weakest link’* (subreach 8). 


loads hi and w w, were weaker, then the difference between system P, and ‘ ‘weakest — 
_ link’? P, would be even greater. This statement can be justified by replacing 
SS w) in Eq. 3 by a joint lognormal DF, which contains r explicitly. in 
_ 2. With soil-mechanical uncertainties, the mean and standard deviation of 
system P, estimated from 100 samples are, respectively: E(P,) = 0.0865 and 
. 4 D(P,)= 0. 0161. Note the considerable change in P ,; in fact, the soil-mechanical 
rer a increase the expected P, result in a non-negligible variance of | 
e | With load parameter (or hydrologic) uncertainties, the mean and standard 
deviation of system P,, also estimated from 100 samples are, respectively: 
E(P,) = 0.0676 and D(P,) = 0.0276. Note in this case the large variance — 
of P, but also a smaller « change in the mean of P, than in case 2. wagh ay 
With both soil-mechanical and hydrological uncertainties, there is a 


_ change i in the mean value of P, (0.0887) than i ‘in case 2, and the standard Gveten 


of P, (0. ).0395) i is al: also so larger than i jan in case 
System Reuasuity Avonc Conriuence Reach 
gf Problem Description.—One of the most complicated problems it in water resources 
engineering is to design a hydraulic structare under uncertain beckwater effect, 


example of this situation is the levee system of a tributary along the confiuen — - a 


| 
= the cal¢ 
| 
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HYDRAULIC ENGINEERING DESIGN 
purpose herein | is to a methodology to for the effect 
of uncertain | backwater curves on the Teliability of a "structural: system. The 
~ "structural: ‘system consists of a levee along the tributary reach just adjacent 
to the main river. Note that the methodology can be used for estimating the 
4 reliability of any system operating under similar conditions. Other hydraulic _ 
- engineering examples where backwater effects may strongly influence the 


reliability of the system may be found in the design of irrigation channels and 


- a flood flows down from the tributary. This is a common working hypothesis 
me. justifiable for flat rivers where the duration of peak flow may be 


The differential equation governing the nonuniform seedy flow in “if tributary 


? 


in which z = . the height of the water level; x = the horizontal distance; 4 
“the average flow velocity through a cross section located at x; Q = the 
discharge of the tributary; and K(z,x) = the water carrying capacity of the 
tributary. Eq. 4 can be solved by a common step method, such as a finite oe 
Available ¢ data include time series at a gaging station on the main river near 
the confluence _ and at - another si station in the tributary just upstream of the 
7 end of the backwater curve; also known is the stage-discharge relationship = 
the tributary observation point. From the two time series, partial duration series = 
be constructed and simultaneous stages h and H can be calculated. 
Generally, the random variables, and H, are not independent; statistical analysis 
is used to ) determine a best- fit | joint D DF. ‘The proposed methodology to estimate 
flood protection reliability « can be applied when this fitted DF is either of the: 
Reliability Model.—The reliability model is presented in two parts: 
. P, of a given levee cross section is obtained by partitioning the set of backwater 
curves into two subsets, i.e., _ the curves passing above the existing levee and 
A _ those not passing above it; and (2) the P,. of the whole levee system is calculated | 
by combining failure events in a finite number of cross sections. In this section, _ } 
‘ _ only levee failure by overtopping is considered; the parameter of the flood — 
-_- wave that triggers this phenomenon is the peak height. If other modes of failure 
= - were considered, then other flood parameters (e.g., flood exposure or duration 7 
‘ q of the flood 1 wave) must be brought into tl the picture (see the section “Reliability 
of Flood Levee Systems”’), but the assumption of steady condition would no- 
—_— be valid. Dynamic routing would have to be used to determine the variation 
of discharge in time and space along both the tributary and the main river. 
Forsimple situations, different routing methods have been elaborated (1, 15,22, 26), 
a however, in general, dynamic routing is rather complicated, even for the case 
f b a single flood moving down a rectangular channel of constant width; then, 
the present case, many floods (h, H) are considered. = 


_ Probability in One Levee Cross Section.—Given joint it sample elements 


| 
| 


¥ th, A) a for the main river and the tributary (i = k, 2, ..., m), backwater curves _ 


can be computed using Eq. 8 and water stage sample elements z,, are obtained bd 
for each cross section. Let Z, and R, denote, respectively, the levee height 
and DF of the water le level in cross k; then the probability of 


>z, 
- On the other hand, the reliability of any given cross section means the a ae 
that all backwater curves in that cross section fall lower than the height of — 
ie levee. Dropping the index sake of brevity, this reliability be 
determined in the following way. 


_A sample value of A is taken the height, Z, in the cross section, 


7 discharge Q and stage H can be computed by Eq. 4. Repeating this procedure, 

_ pairs of values of h and H can be determined, which yield a function, H(h). 

To each pair [A, H(h)), there belongs a backwater curve which has a height, a 
z, in = cross section. The reliability of the cross section is ‘then by definition: 


R= 


the reach considered the water level z, is greater than the levee height, y A 


_ Using the union-intersection principle of eT theory yields (24): di 


= U @,.>Z,)|=1-P 


=1-P 
q 


pa 
Since there are several cross sections along th ‘the 1 reach, ‘the p procedure previously 4 : 

outlined for one cross section cannot be used. A direct method to compute 
P, will now be developed. The joint sample elements (h,,H,) and the values 

a of z, given by Eq. 4 are used to calibrate a predictive model of z, for each 4 
cross section. The following linear regression model is used to estimate th this" 


of 2 asaf of b 


in which of and are least- squares estimates. ‘Note 
that, for the first cross bacrsaady gad = h, and fer the last one, z, = H, so 

that a, = 6, = 1 and b, c, =a, = 0. ‘Eq. 8 can be substituted into 

_ Consider the region of the (h,H) plane defined by the inequalities a,h + 
b,H +c, = Z,. This region is bounded ~— the piecewise linear function n T(h) 


"defined by the set of m segments 
k= 1,2, (10) 


| 
| 
d | 
| 
in which f(h, H) = the joint density function of the pair 
™ _ Probability of Failure of Whole Levee System.—The probability of failure 
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ENGINEERING DESIGN 
r the whale system is ne 


in which f = the joint of as 
The integral in Eq. 11 is teadily computed when the pdf f(h,H) is a joint — 
~ normal pdf. In case of lognormal variates, a logarithmic transformation of all 
*, The proposed methodology recognizes ines this existence of the uncertainties that _ 
are generally present in engineering  ie., natural, sample, and model 
The natural uncertainty caused double randomness of floods is 
~ taken into consideration by using the joint fr S(h, H) for the reliability calculation. 
2. Sample or parameter uncertainty originates from the fact that the joint — 7 
y “pdf f(h, H) is estimated from a finite sample, i.e., if another historical sample 7 
were at hand, the value of P, would be different. A Bayesian method (2) may — 
be used to analyze the effect of this sample uncertainty on the value of P, 
as follows. The possible parameters a ¢ of the joint (h,H) are regarded 


the wales of P, isa function of the parameters a: 


= n Tth\a) 


3 
in 1 which T(h\a) = = the function H = T(h) ‘corresponding. to the “parameter, 4 
— and f(A, sidan the joint density function of the pair (h, H) given parameter 
‘The -Tight- -hand side of Eq. 12 defines P, as a function of the parameter, 2 
a, whose pdf is known. Random values of a are simulated and introduced into | 
m a Eq. 12 to compute random values of P,. An empirical distribution of P, is 
thus obtained and may then be fitted to a continuous pdf using standard statistical F 
procedures. Using this empirical pdf with an appropriate objective function, — 
Bayes decision theory may be applied to account for uncertainty and 
obtain an optimum design v value of P,. (12, 13). 
from several assumptions in the hydrological model 
[type of the joint pdf f(h,H)], hydraulic computations (steady flow, lumped > 
_ friction parameters), and system definition (finite number of cross sections, _ 
it may be noted that simulation sua alternative procedure for finding 
a levee design, as shown in Ref. 32 for the same case study problem. q 
Results are about the same as with the analytical method. 
Numerical Example and Consideration.—The method is now applied to estimate _ 
; - the reliability of the levee system along the confluence reach of the Zagyva 
i River, in Hungary. A levee reach of 37.5 miles (60.4 km) betw 
4 


- gaging station and the mouth in the Tisza River at Szolnok is _ analyzed Fig. } ; 


2). Along this reach, wetted cross-sectional areas and hydraulic radii for both 7 


the probability of failure fo | 
7 


sztelek 


ne. 2.—Geographic Sketch of Confluence Levee Reach Investigated 
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tan dar eviatio 


| ze FIG. 3. —Variation of of Mean and Standard Deviation of P Peak iit s along Confluence 


Station is also known. The roughness coefficients for the bed and the flood 
plain has been determined by a special calibration algorithm (30). 
Statistical samples for h and were obtained from observation data of peak 


i 
| 


sages at the duo gaging stations of ‘es: stage at the 
* ‘station. Between the years of 1932 and 1967, the number of sample pairs (h,H) : 
is 68; thus 68 surface curves can be computed. The hypothesis that f(h,H) 
‘is joint normal cannot be rejected by the Kolmogorov-Smirnov test at the 5% 
level. Parameters of the pdf are E(h) = 33.33 in. (850 mm); D(h) = 1. 05 
ia @0 mm); E(H) = 35.78 in. (910 mm); D(H) = 0.54 in. (14 mm); and 4 =a 
= 0.4889. The use of Eq. 12 to calculate the reliability of cross section k = 
= PZ 0.986 or Pp (6) =1 10 
the pre 4 of failure of the whole system is de determined. The 
_ linear regression model of Eq. 8 was estimated for every cross section. The. < 
variation of the mean and standard deviation of z along the confluence reach 
a is shown in Fig. 3. A minimum of the standard deviation occurs between 12. 4 
miles (20 km) and 18.6 miles (30 km). The maximum n residual variance among 
-Tinear regression equations for the cross sections is 0.256 and the smallest multiple 4 
correlation coefficient in ‘the same cross section is 0.821; however, 90% of 
TABLE 2. —Effect of Sample Uncertainty on P, 


yon P,x10?| Remark 
highest P,* 
45 samples were simulated; the mean P, was 3.79 x 10~*, and the standard deviation, 


me Note: Convert meters to feet by dividing by 0. eel 


the correlation coefficients greater r than 0.90, which 

_ Function T is composed in this case of 49 line segments (49 cross sections). 

Eq. yields a probability of failure of 3.61 x 10-* for the whole system. 
_ The effect of sample uncertainty is then evaluated using 45 simulated joint 

samples of the pair (h, =) of size 68—the same size as historical 
data. For each of the 45 ‘samples, a P, is obtained from Eq. . 12, which was 
integrated numerically t by a simulation method. Runs resulting in the highest 


and lowest P, values are shown in Table 2. The mean P, with sample uncertainty af 


found to be 3.79 x 107 which is somewhat higher than 3.61 x without 
sample uncertainty, but the remarkable feature is a non-negligible variance of 


_Concuusions ‘AND RECOMMENDATIONS 
1. In the case of earth works such as flood levees, statistical approximations 


be used to estimate the of the of the failure 


3 


7 
 § 2. Both load and resistance ‘should be considered as random \ variables; the 
- probability of failure is then a function of the ratio of two random variables. ie 
_ 3. The stochastic model of levee failure may yield fairly high values of the 
‘probability of failure P,; these high values are in agreement with experimental - 
4A stochastic approach can also be used in the case of bivariate hydraulic 
ncertainty. The reliability | of a structure under backwater effect, with both 
ends of backwater curve being at a random stage, is used as an example. In * 
‘this bivariate case, both sample and model uncertainty are high, oe 
a 5. Model uncertainty in hydraulic systems can be handled by use of either 


analysis o or r Bayesian methods. le of 
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APPENDIX TATION 
The following sy symbols ai are used in this pape paper: 


of f estimated oom sample; 


acceleration of gravity; 


"resistance: against overtopping; M 
resistance against sand boiling or subsoil failure; 
lower threshold water level for slope failure; 
_ peak flood level and Stage at beginning of | backwater cur curve; 
sample element index; 
= sample mean; 
number of sample 
pdf = = probability density function; 
boundary i in plane (h, 4); 
re time; 
mean velocity; 


ne resistance slope failure; 


DF = (cumulative) distribution func 


iy union of sets; 
= intersection of sets; sand sone Ft 
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The reliability model of an underground hydraulic engineering system is 

oa presented. Specifically, the case of a system designed to protect underground — 

7 spaces (such as mines) against inflows or inrushes is considered. This paper 
may be viewed as a complement of the companion paper (5), which deals with 
reliability models of surface hydraulic systems. In that paper, a brief review * 
of approaches to reliability investigations of civil engineering structures is given. 4 
A more thorough review of the reliability of complex systems may be ven, bd 
in Ref. 9. The importance of investigating — risk and reliability at the ign 
stage has been pointed out in Ref. 13. 
_ Mining operations, construction of underground spaces ; such as san 
tunnels, must often be undertaken under water hazard. In case of karstic aquifers, j 
the forecasting of inflows or inrushes into the underground spaces can oy A 
be given in a statistical sense (1) since the fault structure of the aquifer exhibits vA 

random variations. If the underground space to be constructed or operated 
is located under karstic water level, some inflow control, called heretofore an = 
underground flood control system, should be provided. The elements of this — 

_ system are the drainage equipment, artificial sealing facilities, the water con- 

veyance, sediment settling and removal equipment, and the pumping station(s). 
Since the system load and also its resistance or capacity are likely to be stochastic, 

methods for estimating system reliability will be used. 

a _ Kesseru (8) has pointed out the necessity of using reliability theory i in mine 

water control systems in Hungary, where mining operations generally take place i 
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It is appropriate to distinguish between the a analyses for protection of human 
- life and for protection of the facilities of a mine. In this paper, only the protection : 
_ of facilities is considered; in other words, equipment and structures for life 
_ protection, such as rescue routes, are outside of the scope of f the present reliability — 
_ In summary, an estimate will be found for the probability that the underground 4 
flood control fulfills i its task within a given time interval. 


- cut of faces; (3) block water cut; (4) mine water cut; (5) sediment settler and 
removal equipment; and (6) central pumping station. The protection of faces 
(element | 1) against inrushes -s can b be ef effected in several ways, such as grouting, = — 


> $F) ody 7 ort eager 


FIG. 1. —Elements of Underground Flood Control 


the “INSTANTAN” protection method, which provokes inrushes before full- 
mining scale starts (7), or else, preliminary or posterior sealing. Water collected © 
_or drained in the faces is conveyed through water cuts (element 2) to the block 
water cut (element 3). In other systems there may be special sediment settlers 
and booster pumps in every block, but gravitational conveyance of water and <7 
sediment is predominant in recently designed facilities. Water and sediment 
from block water cuts travel to the mine water cut (element 4) which leads 
to the central sediment settler (element 5), and the central pumping station 
(element 6). Sediment is removed from the ‘settler by special equipment and — 
_ pumped to the surface by hydraulic means. The number of faces and blocks 
increases as the exploitation of the “mine proceeds; the amount of inrushes 
and sediment increases as new underground spaces are opened. There 
_ are several uncertainties in the estimation of the load, i.e., the amount of water z 
- and sediment to be removed by the control system. A stochastic model developed 
earlier (1) is used to estimate the statistical characteristics « inflow quantities a 
: 


water 
ms 6Fséhazard has been selected. The structure of this system, which is shown in 
_ Fig. 1, is composed of six principal elements: (1) Protection of faces; (2) waters 


HY7 FLOOD CONTROL SYSTEM 
dhich constitute the gut for | the reliability This stochastic 
is summarized in the next section forease of reference. = 
4 System Failure Events.—The failure of an underground flood control system 
can be caused by a complex set of natural and engineering factors. It is necessary o 


Operation is concerned. More precisely, ‘‘disturbance of operation’’ and “Flood-— 
_ing’’ are distinguished. The disturbance of operation corresponds to a failure 
event which disrupts o or decreases | mining production but. does not | stop ‘it. On 


Depending on ‘the location. of ‘the failure, the following top events are defined: 

(1) Disturbance of operation because of flooding in faces; (2) disturbance of 

_ Operation because of flooding in blocks; (3) disturbance of operation in the 
mine with simultaneous disturbances in several blocks; and (4) flooding of mines 5 
In the : next section, the stochastic model of loading d due to inrush events is 
System Load. —System_ loading ; as a result of inrushes can be. characterized 
by the following three quantities: () q= eye, ay 9 or yield of inrush — 


wer area A 


—(Q2) Ymax = MA maximum inrush event ‘yield | over area A; and (3) Q(A) = Saas 


‘yield of inrush events over area / A. The ‘pdf « of these variates can be estimated 
1. Yield of inrush events. —In a reasonable hypothesis based on ‘a 


~~ q follows a lognormal distribution. Observation data tend to support > 


oa Maximum inrush event yield over area A. pe N(A), , the number of 4 
‘inrush events occurring over an area A, follow a Poisson distribution with mean 
AA, which is a second hypothesis based on phenomenological reasoning and 
supported by observation data. Then the distribution of Ymax i is derived from 


distributions of )a and as follows: 


Total yield of —The total yield for 

sum of inrush event yields, and is thus calculated as the sum of a Poisson : 


The DF of Q must be determined from the | DF of q and since direct 
data « on Q are ‘rarely available. There are three possible « approaches to estimate 
the stochastic properties of the analytic ‘method, d, the 


value and variance of Q, but leads to large errors of estimation of small 


ical difficulties. 
‘The first-order analysis consists essentially in calculating directly the 


| 
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| probabilities. Thus, as stated in Ref. ‘AL, a simulation ones ‘such as the one 
; used in Ref. 11 appears to provide the only practical method to estimate the 
é pdt of Q. The statistical dependence i in space of inrush Levents should be examined — 
- _ The yield of inrush events can be controlled by various methods such as — 
cuulien, grouting, and freezing (10). In fact, element | of the system, the protection c 
of faces, includes such yield control. The impact of the control method on 
-inrush yield should be considered in the reliability analysis. A decision rule 
impact expressing | the effect of control on a given inr inrush 
{ 
in which qi = controlled yield of inrush in face (i,j); f,(-) = impact function — 
in face (i,j); q, = natural yield of of ev event; and ind 4,,, 43), 43, = parameters - 4 
an example, consider the method of sealing with a decision 
If gy < ayy do not seal; then cos 


‘fay < qy < ay, seal a portion of the inrush yield: q, = =ayq ait 
3. > aw seal the total inrush: qy = =0. 


‘The ini function is considered in the following reliability model although — 
- notation q’ is not used. Also, uncertainty of the control can be accounted for 
by introducing a random error with a known distribution into the impact function. 
- To sum up, the calculation of system loading, based on a stochastic model, 
leads to the calculation of the DF of inrush yield, the DF of the maximum 
_ yield of inrushes, and that of the total yield of inrushes. The effect of the 
decision rule of control on the inrush yield can be accounted for by an impact — 
function having three parameters, 
_ Simulation which appears to be a practical method for implementing the model, 
calculates the block yield of ‘inrushes, given block i with number n of faces, 


in which Q(A,) = the total yield of inrushes for face (i,j), “characterized by 


Similarly, the total al yield of inn inrushes f for al mine having m m blocks | _ 


pans 


3 

Reliability Calculation.—Reliability models for ‘each event defined in the section 

“Systems Failure Events” is given in Ref. 3. As an example, the reliability 


of mine flooding is now calculated. _ Flooding of the mine occurs when any ~ 
the following six events occurs: bs brave Dae 


OLE Event E. —There is flooding in in every block o of ‘mine: 


i 
i 
| 
7 
q 
| | | | 
| 
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Lae 
“4 
in which D, = (Q,; > a,) = the event of flooding in block i; Cy, is is the le capacity — 


es Event F. .—The e sediment removal capacity CH of the mine is smaller than — 


sed 
the maximum ‘amount of sediment inrush s. Empirical evidence shows ——- 
_ a linear statistic relationship between this sediment volume and the yield ae iF ; 


in "which ky o=a specific sediment yield (in tons of per | ‘cubic meter 
of water); and ¢ = an error term assumed to be distributed normally — 


_3. Event ER. =” failure of the mine sediment removal equipment occurs: 
in 1 which | B = 


of mine 


/ 
Flooding in blocks 


(of mine 


settler water cut 


FIG. 2.— —Failure Tree fe for Mine Flooding 
and t = the time considered. The variable B is an exponential variate which — ; 
is characterized by \,, the average number of failures within unit time interval. nas 

4, Event G. —The total yield of mine inrushes, ‘is greater than the capacity 
of sediment settler (clement 


> 


aa Event H.—The total mine water yield, Q*, is larger than the rae 
cv of the mir mine water cut (element nt 4): 


| 
a 
0) 
nt; 
| bd 
7 ation — 
; 


~ 
This event may be caused by an excessive flooding in the mine, or failure 


of some of the pumps; in either case the real capacity, y of the pumping <a 
decreases. 


tm parameter A; “(the average failure rate e of one pump). The number of pumps — 


remaining in operation is a binomially distributed variate and y is the 
‘a this variate and the nominal of one 


1 —Data Estimation of System Reliability Mine Manyi 


q 


4a 

we 

INRUSH EVENTSS 


Statistics of 

Lognormal 


Number Distribution 


re per | in in 
blocks” ck ilometers® | ki r°| minu meters| seconds 
(8) 


“Data type |. 
{Data type 
Data type 3. 
Note: Data type 4 Gy 0) = 
= 5 m’/min; and cf, = water conveyance capacity from faces = a m*/min. | 
— type 5 includes Cv, = capacity of block water cuts = 550 m*/min. ame 
Data type 6 includes Onn = maximum amount of block water (Eq. 5) = = 180 m aden 
ky = specific yield of sediment = 0.1 ton/m’; S*(¢) = variance of sediment cman 
= 0; CH = sediment removal capacity = 3 tons/min; CQ = = water capacity of the = 


control sediment settler = 170 m’*/min; CV = = water conveyance capacity of the mine 
water cut = 350 m’/min; nm = number of pumps = 16; C = nominal capacity | ofa 


pump = 11 m’*/min; A, = specific number of failures of a pump: 0.11/yr; Ay = specific 
number of failures of the sediment removal equipment = = 0.03/yr; and 2Q,,, = ‘maximum — 
total amount of mine water 300 


flooding can be written as: 


The failure-tree of mine flooding i is in Fig. 2. 


‘The probabilities of failure events are estimated by a Monte Carlo ‘ciutatation ial 
i method. data for this simulation method should be given 


period, 
— 
Fr q 
q 


ar {8} vididsiiss .o% etimil sonebil delidates sd bt sow 
Let the reliability model be applied to the Manyi mine to be located in the 


_ level. Two alternative mine water strategies are considered: the first alternative 
is the application of “INSTANTAN’’ protection on combined with artificial sealing, 
while the second alternative consists in a passive control method. Passive control _ 7 
_ means that inrushes may enter the facilities in which a drainage system has | 
: been provided, but added costs emerge as a result of production slowdown. 5% 
‘Fig. 3 shows the effect of the artificial sealing ‘Strategy on inrush yield 
_ probability. The upper distribution, Fig. . 3, (a) reflects original inrushes; distribu- 
tions, Figs. 3 (b) and (c), refer to two different sealing strengths. Input data 
mobam ow? Yo si! to asd? ei to yiilidedorq 
cient 


Fm q 


6. 3 3. —Effect of Sealing on Probability Distribution o of Inrush Yield: (a) ) No Sealing: 


are given in Table 1. Appendix | I gives “sample simulation Tesults- using the 


_ Generally y speaking, : system reliability of engineering structures is very sensitive 
a the accuracy of loading statistics because the numerical value of reliability 


is generally found in the upper tail of the distribution, where probability estimates 
are most uncertain. In the present case, the mean and variance of inrush yields | 


- numerical value of reliability. At the same time, these statistics are themselves — 
uncertain, since they are based on regional data, experience or at best small # 
3 samples. A Bayes approach can be used to account for this parameter uncertainty _ 


(2,3,4,12). The use of Bayertan distributions for the pdf of different loading 


= the specific number of inrushes are parameters that strongly affect the © 


ransGanuDian region Ol ungary al a Gepih Of 5UU m Delow the Karstic water 
| 
| 
. 
E | 
a 
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doniianditi accounts forth the effect of this parameter uncertainty. Another aia a 
would be to establish confidence limits for reliability (6), $= 
‘The reliability of the systems considered can be augmented in a number of 


ways, such as redundancy, maintenance, and standby. It is east to include 


v-*. % Toi investigate the reliability of underground works such ai as mines, es, statistical | 4 


, 2. Both load and resistance should be considered as random variables: - the 
oe probability of failure is then a function of the ratio of two random variables. _ 
A failure tree analysis followed by Monte Carlo simulation is the eppeoach 
; recommended for studying complex failure events of underground structures 
subject to hydraulic hazard. The reliability analysis is an indispensible part 


the optimal design of the protection scheme. 
2. 
of this research was the of a 
research project between the University of Arizona and the Institute of Cultural f 
Relations in Budapest under National Science Foundation Grant INT 78-12184, 
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Appennix —Reuasiury Ana ysis For 6- 16 


04103 0.04103 


0.04103 0.06103 0.04103 


“For th the flooding of faces, the probability, p = 0.0 everywhere. poe 2 
For disturbance of operation in blocks, p is as follows: 
Block Block 2 lock 3 Block 4 Block 5 
@ For the disturbance of operation in at least one. block, p = 0.00504 Ker ane ain 
flooding in at least one block, r= lia? a si 
For disturbance of operation in the mine, p = 0.00001. 
_ The probabilities for events leading to mine flooding i is as follows: Ie ult bow 
= 


| 
| iV 
cial —The disturbance of operation 
00 
| 
| 
| 
| 


Using Passive Protection in —The d of ‘operation in faces is” 


Block = Block 2 Block 4 


‘Block 
0. 99809 0.99809 30146 30146 


For block ‘flooding, p= 0.0 everywhere. = HS 
flooding i in at least one block, p= = 0. 0. voltae bad 


0.24270 1.00000 0.0 volt =< 
or mine flooding, pe 00000. 
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the disturbance of operation in blocks; 
7 
_ The probability for events leading to mine flooding is as follo ~~ rn 
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The following symbols are in this paper: 


= sample standard deviation; pe ad 


D, flooding event in blocki; 
= (cumulative) distribution 
E = flooding event in every block; 


= impact function in face (i,j); 
i sample element ind index; 


4 specific sediment yield; tt 3 
L = event element of mine flooding; — hy vile 


= evertotmine flooding; 
7 = number of sample elements or number of mining blocks; G + 
number of cross or subreaches or mining faces; 
= probability density function; ‘nahin 
yield of controlled inrush event; I 


maximum yield of f inrush e events; 


R=1-P, reliability; * Do shi 
failure event; a 


mie 


| 
| 
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J 
4 
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= failure time; 
error term; ac 


average failure rate of sediment removal equipment; 


= av average failure rate of one p pump; and 


cea 


‘om 
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classified ; as potentially dangerous to life and property in ‘the event of a hve een 
_ (23). The majority of these dams are earth embankments, and failure of earth — 
_ embankments is usually associated with some type of breaching. Despite these 

impressive statistics, very little is known about the breach itself—its inception, ; 


The few documented earth dam failures have shown one striking common 
denominator: the failure is anything but sudden. For instance, in the arr o 
case a breach will form and grow gradually under the erosive action of the — 

waters. This gradual failure of an earth dam is of particular interest to disaster 
relief planners, because the rate of growth of the breach strongly influences ; 


be the assumed ator instantaneous) failure, 
4 although the former probably amounts to only a few hours in most cases. _ ar 
The objective of this paper is to present a mathematical model of the gradual — 
- failure of an earth dam. _ The main body of the paper contains a review of | 
=, the pertinent literature at and da description of the model f formulation. The mathe- 
matical model is tested using data from the Huaccoto dam failure, which occurred | 
At this stage of the model development, t the overall objective is to achieve 
a physically realistic simulation of the gradual failure of an earth embankment 
caused by an n overtopping flood ¢ event. t. To the writers’ knowledge, the approach — 


£ an implicit numerical solution "technique. The implicit arate solves s for 

_ water elevations simultaneously at all cross sections over the full channel length — 

: by successive time steps. Principles of unsteady open. channel flow, sediment 
transport mechanics, and channel morphology were successfully combined 

simulate the embankment failure. Further model improvements, , specifically in 


_ the description of the breach 1: morphology will be the subject of additional research. — ; 
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on . Reference is made to papers by Ritter (16), Dressler (5), and Whitham 

7 (25), to name a few. Instantaneous failure causes a positive wave in the 

downstream direction 4 wave in n the upstream. direction. . As evidenced 
by 


far from reality in n the case © of : a , gradual breach. In their investigation of the 
_ Teton Dam failure, Brown and Rogers documented a control section that formed _ 
“just upstream of the breach, and described how the reservoir level dropped - 
uniformly upstream of this section, preventing the negative wave fi from Propagating — 
upstream. The duration of the Teton breach was approx 3 

~ Cristofano’s work is perhaps the first attempt to simulate the growth of Z 
‘breach i in an earth dam. Using geotechnical principles, Cristofano ‘equated force 
of water flowing through the breach to resistive shear strength acting on the 
Be ex surface of the overflow channel. Thus, he was able to relate rate of 
change of erosion to rate of change of water flowing through the ae 
channel. This analysis le led to an algebraic equation amount (of eroded 
‘material to flow of water ‘through the breach. ow’ 

_ Cristofano assumed that the breach top width would | remain constant over 


ono 
hes the late 50's s, the United States Army Corps of Basincers Waterways _ 


Experiment Station (WES) used Physical models to conduct an extensive 
correlation was as found between peak ‘outflow from a suddenly breached ion ; 
and a shape factor describing the geometry of the breach. The WES sched dam 
support the conclusion that the Froude number based on peak flow would reach © 

a value of 0.29, which verifies the Schoklitsch ‘equation | (19) for wes: outflow | 
from a a sudden dam failure. vor 
_ The WES findings do not apply to the case of a gradual ities of a an _ 
dam. The experiments were carried out in a laboratory flume of specified bed 
slope, cross section, and roughness characteristics, and the failure was simulated _ 

by an almost instantaneous removal of part or all of the dam. Therefore, while © 
the tests were ‘representative: of sudden failure case, the results cannot be 


associated with gradual failures because the hydraulics of the two cases are 


Prince, et al. (15) of the TVA. have reported on models that use the 
_ for sudden dam breaches developed by WES. They do not apply their models | 

to earth dams, but rather to sudden failure of large gravity dams. Su and Barnes 


w 
| 
| 
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| a 
failure process. In addition, he fixed side slopes of the breach equal to angle 
i: repose of the bank material, and bottom slope of overflow channel equal 
a to angle of friction of bed material. However, the use of an arbitrary constant 7 ._ 
g in Cristofano’s formula renders it, in effect, empirical. In spite of this limitation, — q 
_ the formula has been used in preliminary estimates of dam breach potential, . i 
| 
| 
; : that both resistance and cross-sectional shape were significant in determining | = 
I Brown and Rogers (1) developed a com 
— 
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by I Harris and Wagner (9), in which the was used 
4g compute suspended sediment. They considered the failure of an earth dam 
immediately upon overtopping, degradation of the breach, and erosion to datum a 
_ Brownand Rogers make several rather lucid statements regarding the mechanics 
of breach development. They point out the need for incorporating lateral erosion | 
in the model simulation. They also address some of the differences in modes 
@ failure for exceptionally high earth dams as opposed to long, low embankments. P 
| . addition, they point out that the bulk of the material | eroded from the breach 


Fread (7) has substantially contributed = of 
in recent years. His doctoral dissertation dealt with a dam breach model which 
used the method of characteristics as its numerical solution scheme. The most 
current version of Fread’s model uses the four-point implicit finite oe 
_ Fread assumes the rate of ‘growth of the breach to be time- me-dependent, with 
either rectangular, triangular, or trapezoidal shape. He accomplishes this by 
considering vertical erosion to take place at a constant, predetermined rate. 
This scsumptice is convenient because it allows the time scale of the phenomena 
to be fixed a priori. However, it renders the: model | incapable of Predicting 
the breach-induced flood wave properties. It can n produce a range of flood events 
for | a given range of vertical erosion rates, but which flood event is likely 
to occur is not discernible. Fread does indicate that the outflow hydrograph 
is extremely sensitive to the chosen rate of vertical erosion, but assumes any | 
os errors in prediction to be dampened as the sharp wave moves downstream. wi 
 Fread’ $ main concern is the downstream valley routing of the flood wave 
- from a ‘a breached earth dam, which is the ultimate goal of any | investigation 
; _ of potential dam breaches. However, it should be emphasized that rate of erosion 
and mode of failure of dam determine to a large extent the shape and duration — 
of the flood wave. Until this mechanism is better understood and properly - 
described by mathematical modeling, Fread’s approach can be considered only 
as an approximation giving a range of probable events. ae ve a 
ol For a useful description of dam incidents i in the United Stat es, reference 
is made to an ASCE-USCOLD publication (24). This report contains a well 
- documented chronological listing of dam incidents involving major dams dating ¥ 
back to the late 1800’s. In striking contrast to the general concern for dam , 
safety, this book actually documents a superb safety record for United States 


The described herein | has five major components: (1) Numerical 
(2) channel resistance; (3) sediment routing; (4) breach morphology; and (5) 


Numerical Scheme. —Extensive literature on ‘unsteady flow y and numerical — 

“models justifies using an implicit solution of the dynamic wave to simulate 

a gradual breach. The rate of rise of the outflow hydrograph is fast compared 


the usual stream channel flood waves. Therefore, its 


| 
| 
| 
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should consider the in the equation of motion, i.e 
a _ The four-point implicit scheme, also referred to as Preissmann’s scheme (1! 1), 
is used in the simulation. In theory, this scheme is unconditionally s stable. Ih a 
practice, however, it is often necessary to introduce a weighting factor to damp — £ 
es out disturbances of the same scale as the grid size. In addition, proper specification “ 
“ of initial and boundary conditions is crucial to the stability and convergence Ys 
of the numerical solution. success of the -Preissmann scheme is well 
Channel Resistance. —No closed form expression for channel resistance is _ 
_ available, and empirical equations continue to be used in practice. These equations © ate 
have been derived for equilibrium flow conditions, and their application to ed 
_ The model uses the Manning equation to describe t the resistance of the channel 4 
bottom. For simplicity, the Manning coefficient is kept invariant in time and 
space. An improvement which may hold some promise is to relate the resistance 
an coefficient to the hydraulic and sediment characteristics of the overflow channel. as 
ig Sediment Routing.—The sediment routing component has two major features: 
(1) Numerical solution of the sediment continuity equation (the Exner equation); 
at and (2) computation of bed material transport using a bedload formula suitable 
- for application to high Froude number flows. In this particular application, — 
g the time scale of water and sediment waves is likely to be the same because 
. Mg the channel bed changes are taking place at a rate comparable to the changes - 
i in water surface elevation. Therefore, the same time step is used for water 
wa The choice of a suitable bed load function applicable to high Froude number wt 
flows 1 remains primarily o one of computational convenience. Phenomenological 
, considering the physical detail 
7 ‘that can be incorporated into the simulation. However, the complexity and 
7 “associated expense of these models is a definite limitation. On practical grounds, i 
a simpler and yet reliable equation such as that of ete and Miller a 
EE which q, = bed material transport rate by weight, per unit of channel width; 
7 = bottom ‘shear stress; t. = critical shear stress (shear stress at incipient 
motion); and a and b = coefficient and exponent, respectively. A value of 
7 7 b = 1.5 is widely recognized in connection with high Froude number flows, 
the value of ais s based on experience. 


which y = density of water; d = flow depth; and S, = ‘The 


maa 


critical s shear stress expressed as fo as follows: 


in which ‘density of and D,, = grain size of solid 
For steep slopes and high sediment transport rates such as in the dam a 


Breach Morphology.—The model a ‘component to tie the cross- 
 % | queatiotey of the breach to the outflow from the failing dam. This relation must 
_ function as an additional degree of freedom to represent the vertical and lateral 
development of the breach in a physically realistic way. At its present stage by 
_ of development, the breach morphology component consists of a relation between — 
top. width 2 and flow rate. The existence of such a relation is S rooted in if 
work on regime channels. It is used here for nonequilibrium « channels for lack 
of a rigorous theory of breach morphology. This aspect of the modeling can 
be readily improved as additional research helps clarify the governing morpho- _ 
_ _ logical principles. The recent works of Parker (14) and Chang (3) are tt ae 
~ Initial and Boundary Conditions.—Correct specification of initial and boundary 
ad conditions i is a crucial aspect of modeling. eee acme convenience dictates 


small but finite flow is specified as initial condition « on the downstream side a 
of the dam. Hydraulic conditions are determined using steady flow principles. __ ' 
The upstream boundary condition is the Teservoir water Surface level. In 


0 computation isa function of the specified s size of the rivulet. Progressive erosion | 

a widens and deepens the rivulet, increasing outflow and erosion rate in a 
self-generating manner. The upper cross section on the sloping downstream 

face ‘‘creeps’’ upstream across the dam top until it reaches the upstream face, 
whereby rate of fic flow and erosion increase at a faster rate. 
If outflow increases enough to lower the reservoir level faster than the channel - 

_ bed erodes, both outflow and erosion gradually diminish. Of course, outflow 

_ will eventually decrease even if the breach bed erodes all the way down to_ 
the stream bed. This mode of failure creates the outflow hhydrograph in the 

_ shape of a sharp but nevertheless gradual flood “a oe 7 

5 ‘The mathematical model assumes an overtopping failure of an earth embank- 

ment, whereby a breach begins to grow at some low or weak point of the 

crest and downstream face. This is a likely form of small embankment failure. 

- The actual task of modeling by computational means can be quite demanding. { 

Many calculations must be assembled in a logical and stable framework to 

achieve a workable numerical solution, 
_ The model has the following major parts: (1) Input data and initial conditions; - 

— (2) cross-sectional characteristics; (3) matrix coefficients and double sweep 
solution; (4) | breach geometry descriptor; (5) sediment transport; and i (6) reservoir 
‘- routing. The model simulates spatial and temporal variations of hydraulic and © 
eross-sectional characteristics along the breach channel, from inception until 

‘either the reservoir has drained completely or the breach has stopped eroding. — 

Additional details on | the various parts of the mathematical | model given 

“Input Data and Initial -Conditions.— —Subroutines. are used to data such 


ad 
a than the critical shear stress. Therefore, Eq. 1 can be simplified to the following: _ 
| 
| 
| 
| 
| 
| 
| 
| 
| 
fi 
| 
= 
we 


ve 


scribed the following ‘algebraic relations: 

7 in which y = stage; A = cross-sectional area; P = wetted perimeter; and a 


and a,, B and 5, are coefficients and exponents, respectively. A least square 
gy is used t to determine the coefficients and exponents of Eqs. 3 and 6 as 
a function of distance (along the breach channel) and time (achieved by ‘apdating is 
the cross-sectional parameters afterevery time step), 
_ Matrix Coefficients and Double Sweep Solution.—These model components com 
calculate the coefficients of the discretized Saint Venant equations : and solve 
te resulting matrix by the double sweep solution technique. For a sate 
treatment of the discretization following the Preissmann scheme, reference is 
_ Breach Geometry Descriptor.—To properly describe the gradual breaching of 
a an earth dam, the geometry of the breach must be related to the ang 


i wave. They state that the typical breach it is between one and three times» 
as wide as it is deep. However, no documented research exists to relate breach 
"geometry to breach hydraulics. For the sake of simplicity, the writers initially - 
_ assumed the breach top width to remain constant. After a stable and physically va 
realistic simulation was achieved, the constant width assumption was relaxed, 
4 and instead, a relation between the breach width and flow rate was specified. 
i relation was applied from inception to peak flow, after which the breach 
width was kept constant and equal to the peak flow width. ep +t shneees Bit 
_ This descriptor of breach geometry provided satisfactory results and allowed 
"greater flexibility and an altogether better simulation. Future model upgrading 
should focus on identifying the underlying erosional mechanisms and ee 


them into the simulation scheme. 


_ Sediment Transport.—The model uses the equation of sediment continuity 
(the Exner equation) to » donate the erosive processes at the channel bed. The 


in which Q, “= sediment transport rate, by weight; P= porosity of material _ . 
§ forming the bed; B = average channel width; and z = bed elevation. Ea. 7 
enables the calculation of change | in bed elevation as a function of prevailing 
‘sediment transport rate and bed material p properties. 
To describe sediment transport rate, the Meyer- Peter and Miiller equation . 
_ is used. On this basis, a simple yet reliable relation between sediment transport 


rate and mean velocity can be postulated. i aig breach erosion rate _ 
o- be related to breach hydraulics, ” 


Reservoir Routing.. This model componer nt the u; upstream boundary 


as initial area, width, depth and Tate, dam and breach geometry; soil J 
characteristics; roughness coefficients; spatial and temporal resolution; etc. 
/ 
| 
I 


in which AH = rise / drop in reservoir level from time level n to time level 
nt 1; Q.u. = flow rate at upstream section of breach at time level 7; i. 
= inflow to the reservoir at time level n; At = time step; and A, = reservoir 

Surface area at time level n. Eq. 8 amounts to al local of the upstream 
boundary condition; its “accuracy is ‘contingent upon its use with reasonably 


The aforementioned upstream ‘boundary condition computed as indicated 
on “provides the mechanism that “are the development of the breach. As the 
. upstream flow depth increases, the flow rate and velocity also increase, leading — 
; to the breach h enlargement. The Process is halted when the Teservoir is depleted 

to the point \ where the upstream flow depth and velocity | can no longer sustain 


Case Stuov Huaccoro ‘Narurat Dam, Peru 


7 - The Huaccoto natural dam was formed as a result of a landslide ‘ai occurred © 
Cochacay Creek, tributary ¢ of the e Mantaro River in Central Peru, in April, 
: 1974. The embankment created behind it a huge reservoir which took approx : 
60 days to fill. The failure occurred by overtopping, whereby a breach developed — | 
— and grew progressively under the erosive action of the waters. The large volume 
~ of material encompassing the embankment made this particular failure somewhat 
-unusual—Rodriguez (17) has stated that there was enough material to ‘build bes 
an engineered dam about three times as high in ‘the same location. First, the 
breach did not erode to datum, but rather, erosion stopped when the available _ 
; water in the reservoir was insufficient to sustain further breach growth. Second, — : : 
the duration of the failure was approx 48 h, including breach development and > 
reservoir depletion. This is in sharp contrast to documented failure times for — 


engineered dams, which usually take at most a few hours to erode to ¢ datum. | > 


Of course, the dam thickness, and thus, the amount of material to be e eroded 


TABLE 1 —Comperieen | between Estimated Simulated Flood Characteris- 


ties, Huaccoto Dam Failure, June, 1974, Mantaro River, Peru 


‘Maximum crest erosion, in meters 
Maximum flow depths, i in meters 
from 0 s-100 m */s, in hours 
- Time from 100 m wi] s to peak discharge, __ 
Hydrograph duration, 100 m’*/s 
rising-400 m*/s receding, in 


F 
| 
| 
| 
| 
| 
q 


in | an engineered tain would be much less than that of the Huaccoto dam. am. sie 
“a _ The overall characteristics of the Huaccoto embankment and reservoir are 
(1) Reservoir = 665 millions of m° @) maximum height of em- 
bankment = 
crest elevation = 
size = il mm; and (6) percentage e of ' material finer than 200 s sieve size = = 15%. 
These characteristics and other auxiliary data were used to perform a hindcast 
simulation of the Huaccoto embankment failure. Eleven cross sections — 
- sed, each 300 m apart; the simulation was carried on for 72 h using a time 
. , step of 5 min. ‘The channel friction, , sediment transport, and | channel geometry 
parameters were adjusted within physically realistic bounds. 
Table 1 shows a comparison of estimated actual versus flood 
characteristics. It should be emphasized that the estimated actual values la 
are by no means exact, considering the strenuous circumstances under which i 
_ they were made. . The simulated results agree reasonably well with the estimated 
Recast values, further demonstrating the potential of this modeling tool in the 


assessment of postulated embankment failures. | 


A simulation model of f the gradual failure of 3 an earth embankment has been 


‘the model i is its “ability toa account t for the ; pl of the breach and the eventual 
draining of the reservoir behind the embankment. Concepts of water and sediment 
; - Souting are used in conjunction with a channel geometry descriptor to arrive 
; at a self-contained mathematical model of the breach enlargement and the ensuing 
Zz wave. Unsteady flow elements of the simulation are an implicit nu numerical 
solution “of the complete Saint Venant equations s coupled with a sequential 
sediment routing technique. This approach provides an increased rational basis _ 
for determining outflow hydrographs from postulated earth dam breaches. _- 
The model is tested on the failure of the natural embankment which formed — 
7 in April, 1974 on the Mantaro River in central Peru, as a result of an earthquake- 
generated landslide. Agreement between estimated actual and simulated flood = 
"characteristics i is remarkable, considering the coarseness of the physical detail 
incorporated at this early stage of model development. 
- The need for this type of analysis is supported by means findings of a National 
Science Foundation-sponsored workshop which convened to delineate priorities 
ia hydraulic and hydrologic research. According to the workshop report (2), a 


"government regulations are now requiring that embankments be ar analyzed against 


based on scant data available from a few documented failures. With mathematical 

models becoming more widely available, there is a definite trend toward the 

use of these improved ways of analyzing gradual embankment failures. ie : 
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= coefficient, Eqs. | and 4; 


wetted perimeter; pot wae 
= = inflow to the reservoir; 
flow rate at upstream of breach; 
, = sediment transport rate, by weight; 


= sediment transport rate, by weight, per rr of channel w width; 


S, = friction slope; 


=z bottom level (ort bed elevation); 


A AH = rise /drop i in reservoir level; ~ 
+ = bottom shear stress; and BE 


wy = shear stress. we 
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BOunDARY SHEAR IN SMOOTH 


ai 
‘By Donald Ww. Knight 


The boundary shear stress is a particularly melee parameter in open channel | : 
flow. Tei is 3 also one of the ‘most difficult | to my from a ee. theoretical 


all combine to make the prediction of the boundary shear stress distribution a 

impossible, even for simple cases. Recourse thus has to be made to empirical 

methods. It is of interest to note that in striving to solve the ‘turbulence closure’ 

_ problem, new — has been | expressed in in ‘simple cases” such as rectangular _ 
ary 
_ Shear stress or velocity in resistance, sediment, and dispersion studies. This — 

_ paper therefore presents a method whereby me mean wall or bed shear stress : 

may be calculated for flows in a rectangular channel over a wide range of ¥ 

differential roughness and id breadth to d to depth ratios. The experimental data has 

been contained in two previous papers by the writer (see Refs. 3 and 4), but — 

7 is now presented in a more useful form that will be eae, helpful oi 
those engaged in flume studies. ir? . 


Experimental experiments were in a 15-m 

460 mm wide, set at a constant bed slope of 9.58 x 10~*. Artificial strip roughness 


was used on the bed of the flume at varying wavelengths, A, in order to 
_ differentially roughen the channel bed in comparison with the smooth flume 
- walls. Uniform depth flow was obtained for each roughness spacing at breadth - 
to depth ratios, B/h, of 1, 2, 3, 5, 7.5, 10, and 15. For each depth and 
- roughness intensity, the wall and bed shear stresses were measured by a Preston 
_ tube or by semilogarithmic plotting of the velocity profiles. The percentagt 
'Leet., Civ. Engrg. Dept., Univ. of Birmingham, O. Box 363, Birmingham, 15 
- Note.—Discussion open until December 1, 1981. To extend the closing date one month, 
a written request must be filed with the Manager of Technical and Professional Publications, _ 
ASCE. Manuscript was submitted for review for possible publication on May 27, 1980. a 
_ This paper is part of the Journal of the Hydraulics Division, Proceedings of the American - 
ag of Civil Engineers, ¢ ©ASCE, Vol. 107, No. HY7, July, 1981. ISSN 0044- 
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TABLE 1.—Summary of Experimental Data 


in Shear Forces 
New- 
tons, 


Experi- 
ment Witers 
per 


_ ber | second 


0.187 | 
0.266 | 
| oa | 
0.604 | 0.608 
0.945 | 0.816 
0.988 
o| | 108 
0.242 


702 | 0.810. 


~ 
illi- square square per | per per 2 
pters /ks,, meter meter | meter meter 
7 a 0.011 | 0.123 | 0.134] 83 
5.9 1 0.024 | 0.179 | 0.203 | 12.0 
3 0.04 | 0.202 | 0.248 
ee 0.110 | 0.280 | 0.390 28.2 
476 | 1523 | 0.208 | 0.375 | 0.663 | 434 
90 | 200 | | 0.560 | 0.454 | 1.014 
| 105 | 3004 | | 0.701 | 0.480 | 1.181 | 
| 322 0.022 | 0.145 | 0.167 13200 
200 | 91.6 (0.482 0.632 | 0.086 0.291 | 0.379 | 
222.1 | | 0.973 | 1.081 | 0.432 | 0.497 | 0.929] 465 
380 | 307 | | 0.290 0.007 | 0.133 0.140 | 
a 22 | «5.8 | 45.1 | 1,056 | 0.193 | 0.343 | 0.019 | 0.158 | 0.177 | 10.7 ; 
95 | 614 | 906 | 0.294 0.453 0.036 | 0.208 | 0.244 48 j 
ry 18.5 92.2 | 666 | 0.566 0.594 | 0.104 | 0.273 | 0.377 | 276 
ie | 38.5 | 154.4 | 296 | 0.676 | 0.864 | 0.209 | 0.397 | 0.606 45 _ 
(26 73.8 | 2249 | 162 1.083 | 0.968 0.487 0.445 | 0.932 | 523 
oh | 306.7 | 0.757 0.594 | 1.351 
32 45.0 2,280 0.206 0343 | 0.019 | 0.158 | 0.177 | 107 
3a | oes | 928 11762 | 0391 | 0.636 | 0.073 | 0.293 | 0.366 
| 68.0 | 227.8 1.360 | 0.369 | 0626 | 095 
| 0.343 | 0.024 | 0.158 | 0.182 | 3.2 
68 626 | 7,880 | 0.240 0.571 | 0.030 | 0.263 | 0.293 10.2 
“- 136 | 921 | 5,260 | 0.382 | 0.624 | 0.071 | 0.287 | 0.358 98 
(30.9 | 156.7 | 4,160 0.470 0.931 | 0.147 | 0.428 0.575 5.6 * 
=| $7.7 | 227.8 | 3,500 0.612 | 1.350 | 0.621 | 0.900 | 31.0 
| 303.9 | 3,080 | 0.819 | 1.658 0.763 1.261 395 
34 46.0 — | 0.183 0.375 | 0.017 | 0.172 | 0.189 | 90 
as 0.210 0.795 0.039 0.366 0.405 96 
6 | 235.2 | 9,300 | 0.548 | 1.657 | 0.258 | 0.762 | 1.020 | 25.3 
| 624 | — | 0.184 | 0549 | 0.023 | 0.253 | 0.276 | 8.3 
4 | 104 | 922 | — | 0.236 | 0.757 | 0.046 | 0.348 | 0.392 | 11.2 7 a 
| 1524 15,420 | 0.410 | 1.097 | 0.125 | 0.505 | 0.630 | 19.8 
48.5 | 230.1 | 11,620 | 0.520 | 1.484 | 0.239 | 0.683 | 0.922 | 25.9 
80.0 | 308.4 | 7,380 | 0.778 1.904 | 0.480 (0.876 1,356 | 35.4 
mm | wo} 91.7 | — | 0238 | 0.725 | 0.044 | 0.334 | 0378] 116 | 
76 48.0 | 229.0 17,120 | 0.478 | 1.350 | 0.219 (0.840 | 26.10 
4 $2] — | 0.498 | 0.017 | 0.229 | 0.246 | 69 
10.0 | | 0.028 | 0.357 | 0.385 | 7.3 4 
| 155.6 14,840 | 0.335 0615 (0.719 | 14.5 
| 2314 13,00 | 0.503 1691 | 0.233 0.779 1.012 | 23.0 
846 | 310.7 | 5,220 | 0.708 | 2.065 | 0.440 | 0.950 1390} 31.7 
0177 | 0.418 | 0.016 | 0.174 | 0.190] 8.4 
9a | | ono | 3,280 | 0.239 0877 | 0.08 | 0.403 | 0.44] 96 
* 96 | 48.7 | 223.7 | 3,320 | 0.600 | 1.700 | 0.268 | 0.782 1.050 | 25.5 a 
| | | 1.296 | | 0.396 | | 0.182 | 0.201 a 
$7.9 | 1184 0.231 | 0.027 | 0.240 101 
13.0 | 89.0 | 1,522 | 0381 | 0587 | 0.068 | 0.270 | 0.338 | 20.1 


tet 


30.5 | 149. 1.085 | 0.499 | 0.647 
58. 223.6 | 0. 0.683 | 1.018 
284 | 0.995 0.788 | 1.396 


~_ 


of the shear force carried by either ‘the walls or ‘bed was then calculated ons 0 
¢ correlated with | roughness intensity and aspect ratio. The Nikuradse equivalent | 
sand roughness size, k, , for the bed was also determined. Further details of 
the experimental technique and results may be found in the earlier papers (3,4). _ 
_ Experimental Results.—A summary of the shear stress data is given in hana : 
‘1. The individual experiment numbers are the same as those used previously. — 
_ Values of the depth and discharge are given for "each experiment to together with 2 
‘the mean values of the measured wall and bed shear stresses, 7, and 7,. The 
_ corresponding values of the shear force acting on the walls, SF, (= 2h7,,), 
, the bed SF, (= B7,), and the wetted perimeter, SF, (= 7,P = pgAS,) are 
me also tabulated, in which 7, = mean boundary shear stress; P = wetted perimeter; 
_ p = density; g = gravitational acceleration; A = cross-sectional area; and S, 
. = energy gradient. The percentage of the shear force “carried by the walls, a 
RSE (= SF, x 100/SF,) and bed, %SF,(= SF, x 100/SF, ) are also g given. 
_ Values of the Nikuradse roughness size for the bed, k,,, were determined from — 
values of the bed friction factor and Reynolds number and are shown in 
dimensionless form by dividing by the wall roughness size, k, , taken to be ~ 
- equal to the ‘smooth’ value of 0.0015 mm. The ratio k,. y ist then a measure 
of the differential roughness between the bed and the walls. 


=a results are presented graphically in Fig. 1 in which the %SF,, values 
shown against k, for different B/h values. Although there 


| 
2.9 
106 32.9 
—_ 8.6 
a 


1981 


is a certain amount of scatter, due to the difficulty it in determining k, accurately, 


_ the overall trend is systematic. This « ‘diagram ‘essentially re represents an entirely ; 


mew) way of looking at Figs. 10 and 7 in the writer’s earlier papers (Refs. 
7 3 and 4, respectively). It should be noted that although the results were obtained 
by using artificial strip roughness, the use of the Nikuradse — size 
& 
ANaLysis oF oF Resuits 
Equations for 7,,, and 7, .—The reduction in the percentage 
5 of the shear force carried by the walls, %SF,,, with increasing B/h and k, ,/k, 
; ratios indicated that some analytical function might be found to represent this 
variation. The shape of the curves in Fig. 1, together with wal 


FIG. 2.—Values of %SF, Predicted by Eq.1 } 


channel flow, an of the form ones 


i 


oe" [tanh - -0. 5 (tanh np - 


in -3. log” (B/h + 3) 621; =1- and = 
Ea. i resulted from a curve fitting exercise using the basic function ae 

(7 x) with a correction term -0. 5 {tanh (7 x)- - (1 = 
added to it to fit the mean of the data points for all B/A values. Having obtained - 
' rs relationship y = tanh (w — x) — 0.5 {tanh (7 — x) — (1 — x/7)}’, the. 
variation in the parameter y» for different B/h values was accounted for by 
Som a relationship between y and B/h for the smooth wall/bed data (i.e., 


= 0) in the form given by the equation for a , and the coefficients obtained 


by linear regression. Eq. | results if the substitutions y= 7 %SF,, and x = evs 


are made. The choice of Y= 5 as an arbitrary end point was made on the 


basis of the shape of the majority of the 
q 
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The corresponding for the mean ‘wall and bed shear stresses are 


pghS, 100. t ts Gierefore 


= -= 


and — 


7, by the equation 

The analytical relationship eupesened by Eq. | is shown plotted in Fig. 2 


for the same B/h values as used i in Fig. 1. ‘Tt will will be observed that the equation 


Smooth smooth channel data 


3.—Smooth Channel Experimental Results "Showing Variation of ‘ with 


predicts the same kind tien systematic reduction in %SF,, with increasing B/ h- 
and k,, /k s as shown by the experimental data in Fig. 1. However the precise 
form of Eq. 1 has been chosen in the light o of additional data relating to smooth 7 
>) in smooth 04 open ‘channel flow in which | ra “$5 The writer’s data is taken id 
from Table 1, the Myers and Cruff data from Figs. 7 and 8 in Ref. 7, and 
the Ghosh and Roy data from Fig. 13 Ref. 2. Where precise numerical results ail 
Bs are not available the writer has estimated the numerical values and computed 
the *SF, values accordingly. The main in im concentrating ‘upon 
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5.—Comparison of Experimental and and Analytical | Results od) 


| 
ne: experimental Results Showing Variation of 7, with 
= 


— 
SHEAR 
smooth is because the form of Eq. 1 between the 


influence of roughness (values of ®), and the influence of B/h (values of . 
_ For k,,/k,, = 1, y = 0, B = 1, and the term in square brackets becomes 
; 0. 9963. The equation for a thus governs the numerical values given by Eq. 
4 for different values of y. It is therefore particularly important that any equation 
for a should fit the available smooth channel data well. Figs. 3 and 4 indicate 
that Eqs. 1 and 2 fit t the data reasonably well although there appears | to be 


ss some difference between the writer’ s data and that of other pete The — 


vas 


th Channel Mean Wall ond Bed Sh 
by Eqs. 2 and 3 


in. 


source material, and upon an n equation a in the form 

a=-C, log (B/h + C, 

in which C,, C,, and C, are constants. If the writer’s data is considered on 
its own, then the best fit results are obtained with C, = 3. 775, 4, and 

—C, = 6.909. In the interests of trying to fit the expression to as many results _ 
as possible, the coefficients in the equation for a ae by Eq. | are to be 


is a composite of Figs. | and 2. The reduction in the %SF,, ‘values appears 


to be well | a wide Jk, values. 
of 2, 3, and | 5 may 


q 
om ear Stress Variation as Predicted =f 
| 


be attributed to the choice of coefficients in the equation for a, since it was 
felt that all available data and not just the writer’s should be incorporated 
in @ he analysis. If the writer’s data is considered one its own, then the correlation oi 


t t 10 


maa 


FIG. 7.—Smooth Channel Predict Predicted Mean Wall and Bed Shear as Proportion 
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_ FIG. 8.—Predicted Mean Wall Shear Stress | Variation for Smooth and Rough Channels, i a 


is improved provided that the coefficients in Eq. 5 are changed to those values 
Tew 


a... Mean Boundary Shear Stresses in Smooth and Rough Channels. —Having = 


weeth 
05 
Pe 
J 
q 


a obtained 4 an equation for the percentage of the shear force carried by the well, 
and demonstrated its general limitations and validity, it is instructive to use zs 
a in Eqs. 2 and 3, to obtain the mean wall and bed shear stress ‘variation — 


Then w 


ja Cary wih The 
FIG. 9.—Predicted Mean Bed Shear Stress Variation for Smooth and Rough Channels, 


icy 
af 


as desi 
q FIG. 10.—Predicted SF, | Smooth and Rough Channels ‘Using Eq. 


Figs. 6 and 7 show how the mean wall and bed shear stresses vary with © 
2 B/h for the smooth channel case. The shear stress is nondimensionalized by 
the two-dimensional shear (= pghs, ) in Fig. 6 6, and by the overall mean shea . 


of 

a 

sine ar’ 
OA 

| 
| 


ogRS,) in . Fig. 1, in which R hydraulic ra Figs. 8 9 the 
variation of mean stresses for channels with roughened beds in the same — 


4 could also be produced for the a channel case. 
-. Fig. 10 shows the variation with B/h in the ienateinaes of ‘the shear force 


- carried by the walls for givenk , /ks, _ values. This diagram represents an extension - 
of the smooth channel results given in Fig. 3. 
a ng 11 shows the limiting B/h values that will give 90% or 95% two-dimensional 


11. 11.—Limiting Values of B/ht for and 95% T Two- 


respectively. This Tie equ nt to 


‘equal to 0. 9 or 0. 95. 


| 


OF 


4 is ; somewhat scattered and locking in definition. The analytical "expression 
for a fits the available data with moderate success, although there is clearly — 
room for refinement as more data becomes available. The general shape of 
both curves is however, in good agreement with the data over a wide range — 
The analytical relationships for the mean wall and bed shear stresses shown — 
in Figs. 6 and 7, shows the complexity of boundary shear s stresses in 0] open 
SS The mean wall shear stress, 7,, is higher than the mean bed shear 
stress, t,, for low B/h values, is equal to 7, at B/h = 2.4, and lower than 

-§ at higher B/h values. When nondimensionalized by pghS,, as in Fig. 6, — 

_the mean wall shear stress reaches a maximum at mum at B/h = 7.2. Different analytical . 


F 
ghS, 
Rey 
month hannel Doeulte The sy nerimental data ecente in igs 2 
] 


BOUNDARY SHE. 


- expressions | for a however, » will ll yield maxima at widely differing B/h \ values — 
owing to ‘the flatness of the curve. All the available data suggests that the 
maximum stress is likely to be reached in the region of 5 < B/h or = pee 
a _ The ratio between the mean bed shear stress and the overall mean shear - 
stress shown in Fig. 7 is remarkably constant at around 1.04 for B/h > 6. 
_ The ratio between the mean wall shear stress and the overall mean — stress 
oa Rough Channel Results.—The data presented in Fig. | 1 was derived from flume : 
studies in which the walls were assumed to be hydraulically ‘smooth (k,, = 
0.0015 mm), and the bed was progressively roughened. The ratio k,. Ik, _ therefore 
_ refers only to this condition. Cases in which the walls have a higher "roughness — 
will not necessarily correspond to the same relationship for the same value 
of k,,/k., _. Further studies are needed to elaborate on this point, and in the 
absence of any data, caution is needed in applying Eg. | to any situation other 
é that for which it was derived. Its usefulness however is not greatly impaired bd 
since it is primarily directed towards flume studies which play such an important | 
- "The analytical expressions developed initially from the data shown in Fig. ; 
and given in Egs. are shown to be of considerable practical significance 
7 ‘through Figs. 8-11. Assuming that Eq. | is reasonably correct then Figs. 8 
a 9 allow an investigator to quickly determine the mean wall and bed shear 
stresses for given values of B/h and k, /k, . Fig. 8 indicates that the 


_ wall shear stress is very considerably reduced with increasing k,,/k, _. The 


= values of 7 w /(pghS,) are approximately halved throughout the range "LS + an 


on B/h< 15 when # Jk, is changed from 1-10*. Fig. 9 indicates that the mean = 
® bed shear stresses are considerably higher once the flume bed is roughened, 
and allows a bed shear velocity to be quickly estimated 7” given romghees 
and aspect ratios. 
_ Two Dimensional Flow. —The concept of a ‘wide’ iieaias channel is a 
popular one with hydraulicians on account © of the additional ee gt st that 


meaning of ‘wide’ is of course one of definition. Taking the mean wall — 
force as a possible criterion, then the limiting values of B/h which give only 
5% or 10% %SF,, values are shown in Fig. 11. As would be expected the 
limiting value decreases with increasing bed roughness, e.g., B/h = 22.6-12. 6 
ask, ,/k, varies from 1-104 “for 5%SF,,, and B/h=12.7-6. .6 for the corresponding 
Bjk.< 5, 10*, 2.1 x < 4aR/v < 3.7 x 
: < i/ (gh) < 0.57. The smooth channel results were restricted to the range — 
of 2.9 x 10° < 4uR/v < 3. 6 x 10°, and ( 0.46 < a/-V (gh) < 0.57. The symbols — y 
7 their usual meaning in that i= = mean velocity; and kinematic 


circumstances should the k, 


function in Eq. | has been formulated, i.e., ., for = 10°, 


| 
| 
4 
@ _ It should be emphasized again that the results were obtained in a flume with = fj 
smooth walls and smooth or rough bed. It does not necessarily follow therefore 2 
he rough wall/rough bed case. Under §f 
it 


‘The B/h range is sufficiently wide to define any major changes in the variation 
_of mean boundary shear stress, e. Bs Figs. (3, 4, 6, and 7. It may b be possible im 
therefore to extrapolate the re results “outside the given B/h range with some = 
degree of confidence, although once. B/h = 0.5 it is likely that other effects — 
Because boundary resistance and shear is Reynolds number dependent, at, it 
:a likely that Eq. | will need some further modification before ‘it can be used — 


in any design problem involving prototype flows. = a. 

_ General Comments. —The equations put forward in this paper are considered 

to be reasonably representative of flow conditions in small rectangular open 

channels, although it must be said that the available data is not extensive. 

_ Eq. | is therefore recognized to be somewhat tentative, and likely to be improved | 
upon as further data becomes available. However, the general trends predicted — 
* by Eq. | are sufficiently encouraging to believe that it is a reasonable first — 
_ attempt at analysing this particular problem. Indeed one of the writer’s objectives _ 
in putting it forward is to encourage other investigators to analyse and publish | 

- their data in the same way. It should not be too difficult to investigate the 
smooth channel case more ¢ thoroughly | and for the equation for a to be modified 

accordingly. Yo od cwode on spd baa , 
In terms of application, the proposed equations are simple and easy to use 
on account of the separation of the two dependent variables. There are many 
practical problems involving flows in flumes or rectangular culverts for which 
it is important to know the bed shear stress or bed shear velocity accurately. 
The majority of sediment or Gispersion f formulae utilize parameters exten- 


and ‘roughness. 1-3 should enable such estimates to be n made mor 


bas 


give the percentage of the total shear force ‘carried by the walls 
a flows in smooth and rough channels of rectangular cross section. __ a 

2. The range within which Eq. | is valid needs to be strictly observed on 

= account of the limited amount of data, and the curve fitting procedures employed. _ 
eae The reduction is %SF,, with inc increasing B/h andk,/k, is shown by both 

f the experimental results a and the proposed equation. 

a. . The mean wall and bed shear stresses may be calculated for different a 

aspect ratios and roughness conditions using Eqs.2and3. 

 §, The experimental data for the smooth channel case, which defines the 

_ equation for a, is shown to be limited and not very satisfactory. — ‘alse Rt 

Eqs. 1-3 provide a novel way of estimating the mean boundary 


stresses in flumes and small rectangular culverts. to 
a Further enpetnpenal de data is nome in order to check the general vali vy 
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Aprenoix Il. —Notarion 


The following ‘symbols are used in this Paper: i 
Az "cross-sectional area; 
B channel width; 
gravitational acceleration; 
channel depth; 
Nikuradse equivalent sand Toughness; 
wetted perimeter; 
slope; 
shearforce; 
= =. percentage shear force; 
Section 1 mean 


= coefficient used in 1; 
shear stress. 
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Several stochastic models have been used for modeling hydrologic ithe series 
4 in general and streamflow time series in particular. They are: (1) Ee 
models (15,28 34); (2) fractional Gaussian noise e models (14, 31); (3): autoregressive _ 
moving average models 3, 6,18); (4) broken-line models (16,23); shot- 
3 models (32); (5) disaggregation models (12,30); (6) Markov mixture models (8); 
me ARMA-Markov models (13); and (8) general mixture models (1). Supposedly 
i: justification for using these models follows from their ability to reproduce the 
characteristics which are in actual hydrologic 
Unfortunately, the exact model of a hydrologic time series 
never known. The inferred population model is only an approximation. ‘The af 
exact model parameters are also never known in hydrology; they must be estimated t a” 
from limited data. Identification of models and estimation of their parameters — 


_ from available data are often referred to in the literature as time series modeling 


4 
- time series modeling may be composed of s six main phases (24): (1) ldeatieaton 


of model composition; (2) identification of aon type; (3) identification of © 
model form; estimation of model parameters; (5) testing goodness of fit 


model type and the ‘identification of the model form. 
= The identification of the type of model is aimed at deciding on one among 
a various alternative models, say AR (autoregressive), ARMA (autoregressive | 
moving average), FGN (fractional Gaussian noise), BL (broken line), SL (shifting 
lev el), or any other model that is available in stochastic hydrology. | Once the or . 
type of model is ; identified it is _ necessary to ‘identify the form or order of 
Assoc. Prof., Dept. of Civ. Engrg., Colorado State Univ., Fort Collins, Colo. 80523. 
hy °Grad. Student, Dept. of Civ. Engrg., Colorado State Univ., Fort Collins, Colo. 80523. 
Grad. Student, Dept. of Civ. Engrg., Colorado State Univ., Fort Collins, Colo. 80523. et 
ve ee Note.—Discussion open until December 1, 1981. To extend the closing date one month, 
ar written request must be filed with the Manager of Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication on July 25, 1980. — 
_ This paper is part of the Journal of the Hydraulics Division, Proceedings of the American | 
Society of Civil Engineers, 107, No. HY7, 0044-— 
/81 /0007-0853/$01.00. 
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the model. For instance, if ARMA is is the we type — model selected for a particular 
case, ¢, finding the values Of p and a (the number of autoregressive and moving 


The selection or identification of the type of hydrologic gieciantte model 

- generally depends on: (1) The modelers’ input such as judgment, experience, 
_and personal preference; (2) the physical basis of the process under study; 
and (3) the statistical characteristics of a given time series. The approach developed 

in this paper leans on the second factor, namely, - delineating the physical 


Klemes (11) covers most of the work dane in the “Thomas 
‘Fiering (4) consider a watershed system where the annual precipitation is 
decomposed into evaporation, ‘infiltration, and surface | runoff. Then, using the 
mass balance equation for the ground-water storage, the autocorrelation function | 
of the annual streamflow was found by induction as a function of the first 
serial correlation of the precipitation process which was assumed to follow 
a first order autoregressive model. The streamflow autocorrelation function (ACF) 
given by Fiering (4, p. 74) is a fairly complex expression. ‘However, for an 
assumed independent precipitation process such an autocorrelation function 
7 (ACF) simplifies to p, = p,(1 — c)*~" in which k = the time lag; and c = _ 
a cone, This function being slightly different from the well- known function, : 


process is in fact an ARMA(I,1) Bow Thomas and ‘Fiering’ s 
model of a watershed is used herein to identify the type of model for the — : 
- ground- -water storage and streamflow processes given tl that the precipitation input r 
be represented i in general bys an ARMA process. 
The identification of the form or the order of the hydrologic stochastic model 
- been usually made by applying statistical techniques such as spectral analysis 
and autocorrelation analysis. The books by Jenkins and Watts (9), and Yevjevich — 


5) describe such techniques in detail. More recently Hipel, et al. (6) examine 


FIG. 1 on Representation of Precipitation-Storage- Streamflow Processes 


0 

> resnectiveivi) is the nrocess Of identification of the ora 
a P; » (i) ICU Pieris to deheve = 
‘that it may correspond to a multilag AR model (see also Jackson (7), p. 60). , ; 
However, Salas and Smith (25) note that the ACF p, = p,(1 —c)*~' has the 
* 
— 
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IDENTIFICATION OF STREAMFLOW 
identification | based on the autocorrelation function, partial autocor- 
- relation function, inverse ACF, and inverse partial ACF. However, often such 
~ techniques are subjective and do not uniquely identify the order of the model | 
_ when the model is ARMA. An alternate identification technique is used herein 
is based on the R-functions and S-functions proposed by Gray, et al. (5). 
In the following section of this paper the problem of identification of the 
type of stochastic model for annual streamflows, using a conceptual model 
representation of a watershed is analyzed. Once the type of model is defined, — 
ae identification of the order of the model is presented in the third section. 
An example of model identification is given in the re | section. The paper bl 
ends with a section 1 of summary and conclusions. 
IpenTiFICATION OF Type oF Move. 


4 


which = the precipitation in the year, ¢; = the ground- water si storage 
at the end of year, t; and z, = the streamflow in the year f. Assume that 
the amount ax, infiltrates, percolates, and reaches the ground-water storage 
_ and the amount bx, evaporates from the soil, plants, and surface storage. Thus, 


sm amount (1 — a — b)x, = dx, Tepresents the surface runoff reaching the 


; stream. Assume further that c ‘Ss, , is the ground-water contribution to the ~ 
stream in which S,_ = the ground- water storage at the beginning of year, 
| 
conditions are required: 0 = a, b,c = 1 and 


wa The streamflow z, is made up of the ground-water contribution c S,_, and © 


“the surface runoff dx,. ‘Therefore ak 6 


The mass balance equation for the ground-water is written as 


Eqs. land 2 gives be 
Let »,, w,, and p, represent the n means of x,, z,, and S,, From 
Eq. 2, the expected values and are related asp, =a Likewise, 


from Eq. 3 the expected and are related as p, = (1 — by. 


write x, = p, + z,=p, + and S,=p, + S’, in which 


Xp and S; = = the correspondi ing ‘Tesiduals with mean zero. Substituting these 
expressions into Eqs. 2 and 3, it can be shown that —— — 


—c)S'_,+ax 
and =(1—c)z)_,+ ax’ [d(1 — — ac] 
In the Box-Jenkins notation Eqs. 4 and 5 can be written aS sit? mio? 


in whieh = 1 =( 


| 
a 
| 
6 


in which @ = [dal c) — ac] /d; B an opener ‘such = 
Assume that the precipitation input, x/ = x, By is represented by 
in which 4'(B) = ... — 6) B’; 0(B) = 1 - 0; 


bd 2 and 1; and @/ = the autoregressive and moving average coefficients, respectively 
is an ARMA | (p + 1,q) process. ‘Eq. 9 can be also as 
in which ) = — ,B); and = ae! 


1.0 


which is an an ARMA + + 1) process. ll can be | 
= (1 6,B)0"(B); and n, = de 


be readily derived. For instance, if the precipitati on > ae is an independent 

process, i.e., p = 0, and g = 0 in Eq. 8, or x; = €’, then $(B) = 1 —- 
B, = and 6"(B) = 1 — 6, B, so Eq. 10 simplifies to 


(1 — 6, B)S; =&, (IS) 


(10) = 


| 
/ 
§ 
| 
a 2.—Parameter Space for ARMA (1,1) Streamflow Model of Eq. 14 
| 
3 | 
Prom the several expressions given in ihe nreceding 2 


IDENTIFICATION OF STREAMFLOW 


process 3s with =l-e. Likewise, Eq. | simplifies 
an ARMA al, 1) with = —c;and 6, = [d(1 - —c)—ac]/d. TA 
4 Note that the ARMA processes derived for both the ‘ground-water adie 
and the streamflow processes are what are called here restricted ARMA processes 
in the sense that their parameter space is a subspace of that corresponding 
to the general models. For instance, a, 1) Process of 
Puodila 4 tar 
*-Function for an Model 


ce 


that 0 <= >, = 1, and , = 8, so that the conditions 0 < a, b, c = 1, and 
0 < a + b < 1 are satisfied. Furthermore, it can be demonstrated that ‘the 


parameter space of such ARMA (1,1) process is as shown in Fig. wr a 
autocorrelation function of an ARMA A (1, 1) process is given by (2) 


16) 


+6, 4 — 26, bes 


| 
Undefined 
4 
4 | 


Note that the s sign of Py is determined by t the sign of , — 8, and 1 since , 
= 6, for any combination of a, b, and c (within the parameter space as above 

indicated) p, will always be positive. Furthermore, since 6, = 1 — c = 0, 

the ACF, will also be positive for 2 = Thus, p, decays 


conceptual ‘physical model of a “watershed with» an n independent 
4 input leads to an ARMA (1,1) streamflow process whose autocorrelation a 

is always positive. Note that it is not unusual to obtain negative first serial 

correlation estimates from historical data of annual streamflows. The foregoing — 

analysis indicates that such negative "serial correlations , although statistically 
possible, are not mathematically justified, given the above physical 2 assumptions. — 
_ The main conclusion of this section is that, based on a simple conceptual. 

physical representation of a natural watershed it is demonstrated that the 

ground-water storage and streamflow processes belong to the general class of 

This means that under the linear reservoir sumption the 


eamfl flow -_ 


IDENTIFICATION of Formor Mover 

Once the type of model is identified, it remains to identify the form aA th 
model. That had it is necessary to determine the order of the ARMA nor . 


“identified i in on previous section. This is done following an ain developed 

b> Definition of R-Functions and S-Functions.—Consider an ARMA &, q) model 
- 


= the MA dependent variable assumed with mean 
a variable with zero mean and variance ..., 


the set of autoregressive parameters; {6,, ..., 6 v= the set of moving average 


parameters; and p and q are the number of eutoregressive and moving average 

The R- R- -functions | and S- ‘functions are defined by 6) | + 


/ 
q 
_ time series is justified based on physical grounds; | 
| 
— 1.0000 
a 
Donte: 
7 


for 


TABLE 3.—S*° -Function for ‘ARMA 1) Model = 0.6, 


us 


0.3530 


= 0000 
Of 


0.5000 


0.9167 0.5000 


- The functions H(-) are in tur 


Pr+n— 


| 
f 
— 
2 


in which bh H,(p,) = 1;k = an integer; andn=a positive integer. Similar expressions 
“are: also defined when the p,’s in 20 and 21 are replaced by (-p* Px: in eel 


case the H functions are denoted by H, [(—1)‘p,], andH,,,, [1,(- -1) “px 


_ The R-functions and S-functions can be obtained recursively by - 


k+l 


ees = R,(,.:) 


i 


provided ‘there are no zero 9 divisors. The initial values in Eqs. 22 and 23 are 
S,(p,) = 1; and R,(p,) = p,. Note that the form of the foregoing functions 

- remains unchanged when the argument is th the alternating autocorrelation function, 


(9) p,, instead of p P,- It has been shown (33) i: the shifted S-function given 


y 


TABLE 4. Function, j Pye 30 for Flows of River 


d —0.2007 


0.462900 — 


0.2012 
0.1624 


-1)' ‘Pal = C; (constant), ks k< 


which | C, = ig = p; and k= 
‘st -1)"p,] on 
ifn =p. 
q 
and S* "Pal 


(22) 
/ 
5 
For the ARMA (p,q) model of Eq. 17 with p > 0, the shifted S*-function ; 


IDENTIFICATION 0 OF STREAMFLO 
which 1=0 for the argument p, and / = 1, for the argument (-1) 


__ The foregoing properties can be better illustrated in a tabular form. Table — 
4d gives the S*-function of an ARMA (p,q) model. It has a certain pattern” 
. that is useful for identifying the order p and g of an ARMA model. For instance, 


Table 1 shows that precisely at column p there are two constants (C, up and - 
Cc , down) around 2g nonconstants. The column where this occurs identifies — 
_ the order, p, of the model. Furthermore, since in this column there are 2g > 
-nonconstants, this identifies the order, q of the model. Another feature of 
_ the S*-function is that for n > p, St, ) = too fork =-q- 1 (S* alternates 
in sign as n increases), Ss (p,) = (-1)'C, for k = p + i (S* alternates in 
sign as n increases), and Ss? (p,) is undefined for g + 1 = k = —q — 2. The 
- of this table is elaborated by an example. Note that the infinite values 
always occur above the centerline in Table], 
_ Example.—Consider a particular ARMA (p,q) model as in Eq. 17, with known 
-pand and q and known parameters. The autocorrelation function for such a model 


TABLE 5.—S*-Function with Argument or 1 Pe for A Annual | Flows of Niger | River ia 


oC 


pina la ga) cid) 


~2.2390 -2.5532 | 4.1371 
72.9873 3.6189 2 3996 
| 665.4361 | 0.4419 5305 
28699 | -4.6588 | -10.698 | 9.4689 
-1.5348 -1.2941 | 1.4990 
by Ob, 8576 5883 
| 292. 9447 | -1.6619 | 


is given in Table 2. It is often difficult to identify the order of a model by 7 


_ visual inspection ‘of the autocorrelation and partial autocorrelation functions, S 
- especially for higher order mixed ARMA models. Gray’s identification method — 
follows. The S*-function is obtained from Eqs. 22, 23, and 24 and it is given | 

E in Table 3. Clearly, in column n = 2 two sets of constants appear, the constant 
C= = 0. 5 for k = | and the constant C,=1.0fork = = —2. Since such constants 


Furthermore, there are two nonconstants in column 2, that is, 2q = 
Table 1) or q = 1. In effect, the autocorrelation function of Table 2 and the — 
corresponding S*- function of Table 3 were derived from an ARMA (2,1) model — 
- with parameters , = 1.0, 6, = —0.5, and 6, = —0.5. It is important to note 
the other features of the S*-function. For instances, the infinity appears for 
> 2 and k = —2. The constant C, = 0.5 also appears in column 
1 but with a negative sign. | In general the the constant 0. 0.5 appears wi with acne 


| 
a 
| 
| | 
| 


had sign for n >2 and q= 1. 

oF Nicer River Annuat Fiows 


are used as an 1 example of model identification. Data is in the form of modular 7 
3 coefficients (flows divided by the mean) as appears in Yevjevich (34). y i Ee 
The autocorrelation function, p,, for the annual flows of the Niger River eh. 


a 


4 


the estimated S*- function for = = 3, and 

this table shows approximately constant values above the center line as —_. 
as below the line. There is no other column where such constant patterns appear 
as clearly as in column 1. Therefore, it can be concluded that the order of 
the autoregressive term is one or Jaw It remains to determine the order, Es 

Analyzing Table 5 it should be decided whether there exist nonconstants _ 

around the center line. If it is assumed that no nonconstants exist, that is, 

the constants, C, (down) and C, (up), would start right at the center line, then 

S*-value for n = 2 and k = —1 would tend to +o, However, 

—4.6588. On the other hand, the value S}(-2) = 


‘likely value to represent —o than any other in that second column above | the Se. 
center line. This observation would indicate that two nonconstants should appear ee 
in column |. The values S,(—1) = —2.8699, and S,(0) = —1.5348 would be ead 
_ the two nonconstants. This seems more likely because column 1 shows that 
Ss (- -2) ~ S#(-3) = -2. 158, and ~ StQ) = -1. 864. Therefore, 
~2.15, and the constant 
££ below the line would be about <1. 86. This is supported by the fact . 
S3(1) = 1.86 [about the same but with opposite sign as ST a0) Ee From the — 
: foregoing analysis it can be concluded that 2g = 2 or q = 1. ane 
model for the annual flows of the oer River i is identified as ARMA (1, 


the main subject of this paper. The identification of the type of stochastic | 7 . 
model has been made based on a conceptual physical representation of a natural 


based on the recently developed R-functions and S-functions. 
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CATION OF STREAMFLOW 


-water Storage is an ARMA 39) process and the streamflow = 


ans ARMA + l,g+ 1) process. In general such ground: -water and streamflow 
, processes belong to the class of restricted ARMA processes in the sense that 
«their parameter space is a subspace of that corresponding to the general ARMA 

_ _ The form or order of the ground-water and streamflow ARMA processes 
= for given historical time series can be uniquely identified t by using the R- R- functions 
and S- functions. It is "expected ‘that such a a technique will have a widespread — 
application in modeling of hydrologic time series. 
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IDENTIFICATION OF STREAMFLOW 


The fol pllowing symbols are used in this paper: 
= abbreviation for ‘autoregressive moving average’; 
fraction of the precipitation reaching ground-water storage; 
Be backward shift operator defined by Bz, = z,_,; 
= fraction of precipitation lost by evapotranspiration; 
constant (down) of pth column of table of S*-function; — 
At oe constant (up) of pth column of table of S*-function; — 
c = fraction of ground-water storage released as streamflow; a ~— 
d = _ l-a-b = fraction of precipitation reaching stream as direci runoff; . 
AY determinants used to define R-functions and S-functions; 


_N = sample size of annual flows; 
p= order of autoregressive component; he | 
R-function used in Gray’s method of identification; 
S- function used i in in Gray’s s method of identification; 
= storage at beginning of year, f al 
_ & == symbol used to denote undefined terms in table of S*-function; iy 
- 7 = input precipitation to conceptual model of watershed during year, : 
used as annual streamflow to define autocorrelation function 


deviation of x, from its mean, p wale 


= sample mean of annual steamflows; 
‘ = streamflow during year,t; 
deviation n of 2, from its mean, 
with zero mean and variance, ? (white noise), 

= white noise with zero mean and variancea/ ; at 

= white noise with zero mean and variance 

), =, jth moving average parameter of general ARMA model; _— 
jth moving average parameter of ARMA model for input; ‘a Prue 
polynomial in backward shift operator B with coefficients, 0’ ree 
‘polynomial in backward shift agi B of ARMA model for 
population mean of random variable, S,; 
population mean of random variable, x,; 
mean of random variable, z,; 


white noise with zero mean ne variance, a 2g’ : 


pong 
variance of white noise, 
jth autoregressive parameter of general ARMA ‘model; 


| 
| | 


= jth eutoregressive parameter of ARMA model for input; 


= polynomial in backward shift operator B with coefficients, /: and © 
‘= ‘polynomial backward shift operator of model for 


wes, 


| 
» ie @ 
a 


 gwobdru avon 


water resoure trends is cially 


(7). Flood plain and land-use management efforts around a large lake are eseverely 
handicapped by not having the estimates of future lake levels. =f, 
_ Statistical forecasting is a relatively new "application in resources 
management and it may provide a way to project and plan for the future. : 
_ Hydraulic engineers are currently using a variety of methods ranging all the - 
= wey from subjective intuitive methods to rigorous, cause-effect models such 
_ as estimating runoff from rainfall data. The cause- -effect models depend wholly 
on information pertaining to variables which cause, say, lake levels to rise 
and fall; and are relatively difficult and expensive to develop. Moreover ‘such 4 
_ models may not provide much lead time for the planners to develop preventive 
measures. There exist, however, mathematical techniques that make use Zs 
: historical data and project future values. Such techniques have proved to be 
«very useful and often ‘accurate’ in other disciplines such as economics and ] 
_ The purpose of the present study i is to utilize a forecasting technique ‘called f 
: : Exponentially Weighted Moving Average (EWMA) method and forecast water 


levels of a large lake in Alberta, Canada. The forecasting ability of aera 


_ causes for floods; (2) the flood-prone lake shore lands require advance planning» 
_ for rehabilitation measures; (3) the real estate around the lake shore is always 
in demand for recreational and other in demand for recreational and other development purposes; and (4) the average purposes; and (4) the average © 


__~ T§ection Head, Design and Construction Div., Alberta Environment, 16403-102 Street, Head, Design and Construction Div., Alberta Environment, 16403-102 Street, z 


lake levels due to. hydrometeorological and other nonrelated factors are ‘prime = 


Note.—Discussion open until December 1, 1981. To extend the closing date one month, : 
a written request must be filed with the Manager of Technical and Professional Publications, — 
ASCE. Manuscript was submitted for review for possible publication on September _ 
Ss This paper is part of the Journal of the Hydraulics Division, Proceedings of the 
American Society of Civil Engineers, ©ASCE, Vol. 107, No. HY7, 
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= data for the years 1913-1977 is readily available for Lesser Slave Lake, 
Tt has become conventional to to classify forecasts as being either ex-post or _ = 
ex- -ante. In the latter case, the true values of exogenous variables : are unkriown 

at the time the forecast is made, therefore, they must be estimated. Forecasting — 
errors may, therefore, be due to errors in the estimates of the predetermined — 
variables or incorrect specification of the forecasting model or both. In this 
study the forecasts presented are of an ex-post nature, which uses ectual 
observations of the predetermined v variables to predict values for the endogenous" 7 
variables; any errors are therefore attributable only to the model. | eRe 


_ The monthly data for the period 1913-1972 (see Fig. 1) is used in = 


building and estimation of the parameters of the model. the predictive ability 
_ of the forecasting model is tested by comparing the ex-post forecasts with = 
60 observations of the period 1973- 1977. The predictive accuracy is 


Xe? /60, a and Mincer- ‘Zarnowitz criterion Yow in which e e= = the « difference 


FIG. Average Lake Lesser Slave Lake, 191 


the actual and predicted lake levels. Mincer and Zarnowitz (6) enunciated a a 

definition of forecast accuracy basedon 

in which A(t) = observation at time, F(t) = the corresponding 

and a, and a, = regression coefficients. If a, = 0 and a, tends to unity, — 
the forecast is deemed accurate. the coefficients, a, and a,, are determined 7 
by the application of ordinary least squares. ct onal 


, x(t)} be a a | stationary stochastic se series in discrete 


1). The The tar By g a(t + 1), may be given as a weighted sum of the the past observations 


in which c,, c,, ... = the weights attached to each of the known observations. 


It seems sensible to attach more hy to recent observations and less weight 
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PREDICTING LAKE LEVELS 


those which in PR progression (2). Eq. 2a then be expressed 


“im which = a constant such that0<A<10 ai itasa 


2b implies an infinite number of past observations that are to 
- estimate X(t + 1), however in prastion there will only be a finite number. Eq. 


which i is called ai an exponentially weighted ‘moving average. Thus, for t = 


2c may be expressed as a series of equations: 


which time the observation x(1) is available. If we set %(1) = x(I), i.e., 

_ predicted value of the series at time ¢ = | equals the actual observation at 
that time, then £(2) = Ax(1) + (I — A) x(1) = x(1). This means that the —— 
value forecasted for time ¢ = 2 of the series is equal to the observation at 
time 1 t= = I. Eq. . 2e can now be used ‘Tecursively to compute forecasts at future 


be easily updated using the latest and the — 
Eq. 2c may be rewritten 


4 in x(t) - 20), the prediction error at time, ¢. The of 


the constant, \ depends on the characteristics of the time series. For an arbitrary 7 


chosen the may be calculated: 


com 


‘prose 
and so 


£(t) = he,. it 1) con 


and the expression Zina? is computed. This aeeiniun is repeated for various 


values of ) in the sange of 0-1. The value c corresponding to the minimum: of 


| 
§ 
| 
q 
3) j 
R(3) = re 


is chosen as the optimum dt to ma weed in computing - 
Winters (8) has generalized the foregoing method to deal | with h time series — 
that contains trend and seasonal variations. Let M (t) denote the ¢ estimated current 
mean of x(t) in period, t; T(t) denote the estimated trend, i.e., the expected | 
_ change in current mean; and S(t) denote the estimated seasonal factor in period, 
_t. As each new observation becomes available, the three terms M(t), TO, 
and S(t) 3 are updated. The ‘Seasonal variation in the time series 5 may possess 
pattern be proportional to ‘the level of ‘the observations a multiplicative, or 
ratio, seasonal effect is said to exist. However, if the amplitude is independent | 
. of the levels then an additive effect should be considered. A graph of the 
_ data should be examined to ascertain whether an additive or multiplicative seasonal 
: effect is present. The method may ‘blow up’ if the wrong seasonal effect is 7 
applied. The updating ‘equations for M(t) for an seasonal 
are 


= B {x(t) - (l 


~ 
= 


in which a and 8B are ancuthing constants such that 0 < a; B <1; and Ss 
= the seasonal span (e.g., s = 12 for monthly data). The current mean and 
“seasonal factors are thus updated by linear superposition n of past - values. If 


_ the seasonal variation is multiplicative the updating equations will be srlay 


in which the. current values for M(t) and S(t) are ‘derived ‘through a Process ~ 
of ratio, or division. The — — term 


Eqs. 4a-4e are of such a nature that if a state of a time series is specified 
at some initial time, ¢ = ¢,, a solution exists for ¢ > be and is uniquely determined 
by Eq. 5. Starting values for M(t), S(t), and T(t) may be approximated from 
- initial observations in the time series, e.g., for monthly data, the first 24 


MQ) 


| 
| 
L 


"PREDICTING L LAKE LEVELS” 


S(l)= 
M(2) 
+ 


The three smoothing odpatiie a, B, and ¥, are varied in the range of 0.0-1. 
and the quantity Le? is calculated. The set of values for (a, B,-y) corresponding 
to the minimum of the calculated Le? is the optimum for a, B, and y. y The - 
optimum (a,f,) is used in updating the equations for M(t), S(t), and T(t); 
and forecasts are made using Eq. 5 for the horizon of length, c. eo fe oe | 7 
As mentioned in the | the first. in| the | 
time series are used to develop the initial values of M(t), S(t), and T(t). The 
- next 60 observations are utilized to test the forecasting model—making a forecast; % 
a moving along one period (12 months for monthly lake level data); comparing 
_ the forecasts with observed data; absorbing the observed data into the forecasting © 
model; making forecast for the next period; and the cycle is repeated. er o 
A \ grid of values for the smoothing constants, a, B, and y, is used to calculate 
the predicted mean squared error, = (actual- Seiotest" /60. The grid was made 
= of all possible combinations of 0., 0.2, 0.4, 0.6, 0.8, and 1.0. The values 
of Se’/60, in which e = the error between the actual and forecast | values, 2 
are given in Table 1. Table 2 shows the areas around the minimum Le? /60 — 
fora finer grid values of a, B, cand y with i increments equal to 0.01; as — . 
4 be seen the value of Ze’ /60 is rather flat 1 near its minimum. Perhaps” lower a 
: values of Ze”/60 could be found, however the benefit gained may not be worth 
4 the effort. The lowest calculated value for Le’ /60 is 0.35 at a = 0.80, B -— 
0.64, and y = 0.29 for a multiplicative seasonal model. 
Few general comments about the optimal values associated with the shiities 
constants | can be made. It is possible to to imagine a series in which some 
the constants are subject to little or ‘no drift; and also a series s where there 
is little random effect. If there is no change i in a constant over the series then | 
if the random effect is small, the value associated with that constant would 
be small, or even zero, because there is no use in changing the original and 
still accurate estimate of the | constant. If the drift in a constant is large over 
the series, even to the point of abrupt - changes | in its value from time to time, 
there are two possibilities: (1) Little random effect would lead to large values 
weighting current estimates heavily; and (2) a large random effect would yield 
substantially smaller values depending on the relative importance of the changes — 
_in the constant as compared to the random element. These ‘intuitive deductions 


ized in 
In the case ase of Lesser Slave Lake level dota the etnies constant, 
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TABLE 1.—Value of 27/60: Grid Results” 
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TABLE 2.—Value of Finer Grid Results 
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“ments are quite pronounced | in Lesser Slave Lake levels. Finally the trend 
parameter makes only slight revisions in the smoothing process as is evidenced — 

by a low value for y. This shows that no abrupt shifts in lake level sala 
occur. of the EWMA the actual lake 


7 large revealing that the effect of earlier observations is quickly attenuated. This = 
: _ generally occurs when the mean of the series changes quickly, and this would — 
seem to be the 
| 
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—— ACTUAL OBSERVATIONS 


«1876 
a =| 


is presented in Fig. 2 the 60 from 1973-1977. 
A visual inspection of Fig. 2 shows that the forecasts in general follow the 
trend of the actual observations and the | a deviation appears to be of 


—" 


Comparison with OTHER FORECASTING Movets 
wry 
Two other forecasting methods were to the lake level by ‘the: 
wrleer (5) and the results are used to make a comparison with the results of 
EWMA. The models were eased on Box-Jenkins technique (1) and Harmonic» 


> The results of these forecasting models together with those from EWMA 
are given in Table 3. A summary of the comparison between =e of the forecasts 

and the actual observations is presented in Tables 4 and 5. The predicted mean 

‘squared error, 2e’/60, for the Harmonic | analysis is substantially higher than _ 
‘that of the other three models. The Box-Jenkins and Composite models however 

are more serious contenders; which in some ways have results similar to those _ 
of EWMA for the period being forecast. The Box-Jenkins approach identifies 
stochastic components in the time series, i.e., autoregressive (AR) and moving 
average (MA) components. The Composite model embraces not only the AR * 
and MA stochastic components but also information inherent in harmonic cycles _ 
of periodicity greater than 12. On the other hand, EWMA readily incorporates — Pa 
any drifts over time and can filter out substantial random effects in the — 
‘Table 4 shows that EWMA is consistently superior in its predictive ability 

as evidenced by the lowest values of the predicted mean squared error, Xe? /60, 
the mean absolute e error, =| e|/60. Both these criteria are important because 


{ 
| 
| 
| 
| 


PREDICTING LAKE LEVELS levers 
TABLE 3.—Actual Observations and Forecasts Made over Horizon of Time Period pussies 
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1,895.30 1,896.30 1894.49 q 1,892.68 
1,894.80 895.80 | 1,894. 1894.39 | 1,892.28 = 


TABLE 4.—Comparison of Forecasting Models 


Method Zle|/60 /60 


Box-Jenkins 0,712 0.687 


‘TABLE 5.—Mincer-; Zarnowitz Forecast Comparison 


A(t) = 376.55" + 0. 8012 E(t) 
Box-Jenkins = —254.83" + 1.1345 B(r) 
Harmonic ‘bh A(t) = 255.27” + 0.8660 H(t) 
Composite A(t) = 0.004" + 1.0 1.0 C(t) 
Composite 
a Not significantly different from zero because of the nibatnae a? of standard deviation 


“Equation made to pass through the origin. 


1 


is s difficult to to determine the consequences of the forecast errors in this particular 
case study. Whenever the consequences of one large error are more serious 
than that of several small errors the predicted mean squared error would be 
‘amore appropriate criterion for r comparing the forecasts made by various methods. 
.. Table 5 contains the regression results obtained utilizing Mincer-Zarnowitz 
- criterion. As is evident only EWMA and Composite forecasts satisfy the criterion. | 


| 


PREDICTING LAKE LEVELS 


To further evaluate these two forecasts, the respective regression equations — 
were made to pass through the origin; and the coefficients associated with 


: E(t) and C(t) are found to equal unity and are significant at 0.95 level. The 


= 


grails wad the (ublowieg 


_ multiple correlation coefficient, R’, which explains the percentage variation 
_ inherent in the actual observations, however is different for the two methods. 
The EWMA forecasting model explains 78.5% while the Composite model explains ; 
only 53.6% of the total variations in the actual observations. Furthermore, the 
q 
are the closest to the actual observations. 


_ The forecasting method presented in this paper provi easy, convenient, _ 


and inexpensive procedure to estimate future lake levels. It must be strongly 
emphasized that the object is to provide reasonably good estimates for hydraulic 
and planning engineers who would otherwise may have no flood hazard informa- 
tion at all. Another important application of this paper has been to demonstrate 
the advantages of the EWMA forecasting model over other methods that use 
_ the historical data as the ‘sole input: (1) EWMA requires less information and 
5 ‘storage space; (2) EWMA responds more rapidly to sudden shifts in the time | 
4 series; and (3) EWMA provides routine forecasts without human intervention. a 
‘The EWMaA method is generally applicable to data which contains a minimum 
_ of 50-75 observations taken at regular intervals of time. The nn the number 
. of observations the ‘more’ accurate will be the forecasts, i haa closer will | 


EWMaA model possesses a smaller standard error of regression. Thus a complete z 
statistical evaluation reveals that the results from EWMA are. significant alll 


be the forecasts to the actual observations. 


_ By fitting a mathematical model to the historic lake levels it is demonstrated 
that possible future occurrences of the lake levels can be generated. Furthermore, 
it is shown that the forecast values are statistically indistinguishable from the | 

data. The future lake level(s) is(are) on one of th the crucial cements in 
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symbols are int this paper: 


& 


= Box-Jenkins forecast at time, 


waerdo 
integer; 
Corl constant t weights; 


= seasonal span, e.g., s = 4 for quarterly data; 3S) 

To) = trend factor in period, ¢; be 

x(t) = lake level observation at time, 
 F® = estimate of lake level at time, t; 


= smoothing constant for 


smoothing constant for TO; and 


— 
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ont By Chao-Lin Chiu," M. ASCE and David E. Hsiung’ ban 


mostly deal with one- or two- dimensional analysis and modeling. To treat the 
flow three dimensionally, ‘‘secondary flow’’ consisting of the vertical and 
transverse flow components, must be included, along with the ‘“‘primary flow”’ 
4 (the flow component in the longitudinal direction), in hydraulic studies. It AJ 
also necessary to deal with the “complex reality that channel cross sections, — 
shear stress (including the boundary shear), velocities of primary and secondary 
flows, sediment concentration, etc., all vary three dimensionally, and that these oi 
hydraulic variables interact each other. A change in one will affect the others. __ 
Earle papers (2,3) dealt with various other aspects of three- ~~ 
d mathematical modeling of open channel flow; this paper is concerned with relations 
_ and interactions among the secondary flow, shear stress, 2 and sediment tconcentra- _ 
“g tion in alluvial channels. Like those earlier papers, in hydraulic analysis this — 
paper uses the framework of a curvilinear coordinate system consisting of isovels— - 
_ (curves along which the velocity is equal) of primary flow, marked as & ouves 
in Fig. | and their orthogonal trajectories mar marked as asmcurves. 


Coorpinate System ano Basic callings 
4 


In the &n coordinate system the distribution of primary flow velocity al 


be represented quite well by the following logarithmic ‘equation (Ref. 2, 3,4): - 
in which a suitable le equation for the € coordinate, representing a family of primary 
‘Prof. of Civ. Engrg., Univ. of Pittsburgh, Pittsburgh, Pa 
?Engr., General Electric Company, Rockville, Md.; formerty, Res arch / Asst., Dept. 
Of Civ. Engrg., Univ. of Pittsburgh, Pittsburgh, Pa. 
= Note.—Discussion open until December 1, 1981. To extend the ‘eslias date one month, 
a written request must be filed with the Manager of Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication on September 25, 
_ 1980. This paper is part of the Journal of the Hydraulics Division, Proceedings of the 
American Society of Civil Engineers, ©ASCE, Vol. 107, No. HY7, » July, 1981. ISSN 
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direction); V mear aie 
in which g = the gravitational acceleration, R = the hydraulic radius, and S__ d 
= energy or channel slope; D = the water depth « at the x,-axis; B, for i equal - 
to either l or 2 = the transverse distance on the water surface between the ¢ 


= the mean shear velocity, evaluated as 5 aks at a cross section, 7 ‘ 


approximate ™ 
channel cross-section) ~ 


The mean elevation 


the x5-x_ coordinat 


system at a cross- 


FIG. 1.—x,-é-n Coordinate Sys 


coordinate is in - transverse direction; x, = the coordinate i in the vertical disection 

(the x,-axis is selected such that it passes through the point of maximum primary 
: flow velocity’; ce and B, = empirical coefficients which can be determined 
- according to the variance reducing technique from the linear relation between 
‘ V, and € ona 1 semi-log plot. Eq. | has the same form as the Prandtl-von Karman 

eniversal velocity distribution equation. In addition, « is a coefficient, actually — 


evaluated with observed data on V,, which also tends to be equal to about | 
0.4; thus, if'o no data are peer 0. 4 can be used | as sa A good approximation 
of 


/ 
- 
“yy 
| 


br Orthaganal trajectories of € curves (isovels of of primary flow) can be n be obtained 


_ The momentum equation in the x, direction (longitudinal Gisection) in the 
x ,-€-n coordinate system directly gives the € component of flow, a, , the 7 
component perpendicular to the isovels of primary flow, as: 


-p V7, H 
p — (pg )+ ax, 


+— + 


while the continuity equation gives the n component (the f nis eee tangent a 

+ V, 


In Eqs. 4 and 5, and = =] the: and ‘components of mean ean secondary 
flow velocity; h, ‘ae h = the “scale factors’’ or the “‘metric coefficients” 
on the — and y curves: p = the fluid density; ¢ = time; g= - the gravitational “i 
acceleration; H= the 1 mean elevation of the bottom of a transverse cross section / 
of open channel; 4 A . = the value of V,, at a boundary point in in which a = 
is the “shear stress in the "direction in the plane perpendicular to 
‘the € direction; and o, = the normal stress in the x, direction. iy enewer pnd 
_ The €- and - components of secondary flow “velocity given by Eqs. 4 and 
5 are related to the x, and x, components in the (Cartesian) X>-X, system according 


to the following transformation rule: 


_ To compute secondary currents, Eq. 4 can be used first at every on -_ 
of the €-n coordinate network to obtain V,. The other component, V. 

__ then be obtained by integrating Eq. 5 along each € curve starting from a feast - 
-- point. The boundary point can be selected on the water surface where the | 
>. component of flow, V,, is equal to zero. This boundary condition, V, = | 

(0 along with V, from Eq. 4 can be substituted into Eqs. 6 and 7 to determine 
V, and the x, component, V,. With Eqs. 4 and 5, along with these boundary — 
- conditions, the secondary flow velocities, V, and V,, can be computed at every 
grid point in the €-» coordinate system which can then be transformed into 


we 


_V, and V, in the Cartesian coordinate system by Eqs. 6 and yf a 2 Sindee alt 


(4) | 
| 
| 


The scale on the and n curves can be. detived ‘from Eqs 9% an 


8) 


G 


As indicated by Eq. 4 1, the secondary flow interacts with the 
By definition of & coordinate, t,,, Should be zero as the gradient of primary 
velocity in the y direction (tangent to a curve) is zero; then, + Should 
: represent the total shear at a point in a &y plane (x,-x, plane, or a transverse 


cross section of a channel). Eq. 4 can also be expressed as: 


Seconpary FLow AND Snear Sraess 


One can the solution « of f Eq. 10 10 as od: ome 


in which E,(n) = ‘the value of ¢ € on the water surface along an 7 curve, i.e., 


on each curve the value of varies from &, to which can be ver 


and the last term on the right ‘side is the s shee stress at the water surface ; 
on an 7» curve. Eqs. 4 and 10 show the interaction among the secondary flow ti 
- stress, gravity, hydrostatic pressure, and the variability with time and 
‘space of primary flow velocity. For &, = €, Eq. 12 represents the boundary 
_ Shear stress. If the wind effect is — and oe shear on the water surface 
e x, direction, Eq. - 


pgS 


the primary flow | are horizontal and parallel so that n 
- parallel, vertical lines, e.g., in a “‘wide channel,” h, and h, ty an a 
are bok constant, of Then, h = h, (Eon), h 


88 7 
3 
x; 
| 
j 
| 
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= 


= Dt = l, = x,/D, R = D, = V, = Vz, so 


_If the secondary flow velocity, V,, is 


‘The first integral in 1 Eq. 13 is the. x, ‘component of ‘the weight of enn in 


‘in Fig. 1, multiplied by the unit length of channel in the x, direction. If "E, ; 

= = &, the perimeter, ab, of length equal toh Eo.) will be on the channel = 

bed. So, this integral represents the portion of shear stress due to gravity. 

The second integral represents another portion of shear stress owing to the 
‘net momentum transfer above ab, in the € direction (perpendicular to ab and 
other primary flow isovels above it), by & (the secondary flow velocity in 

the € direction). The direction of Ve may change, positive or negative, i.e., 

- upward or downward. If the net momentum transfer i is in the negative éd direction, . 
ie, , downward and perpendicular to the perimenter, _ab, so that the integral 
9 is a negative quantity, this portion of shear stress contributed by the secondary 

_ flow is to increase the total shear, t,,(&,n). The algebraic sum of the two 

integrals divided by A , (€,n), which represents the length of ab, gives the resultant 
a shear stress under a steady, uniform flow condition. If the flow is unsteady — 
or r nonuniform, or both, the remaining terms in F (defined by Eq. 12, other 
than the gr gravity and secondary flow terms, will also affect the total shear stress. 
_ To use Eq. 12 or 13 to evaluate the shear stress, one needs the secondary _ 
oy flow velocity V,; however, one must know the shear stress to estimate V, . 
Both the shear stress and secondary flow are difficult to measure, therefore, _ 
i it is desirable to develop formulas to estimate them. It seems to be much — 
easier to derive a formula for the shear stress than to derive one for the secondary = 
flow, since the shear stress has a close relationship with the distribution of 
primary flow velocity, as indicated by existing formulas for the viscous and 
turbulent shear stresses in fluid mechanics. A shear stress estimated by aformula _ 
can then be used in Eqs. 4 and 5 to calculate the secondary flow. 
Based on Eq. 12 or 13, tal (é,9) is a function of & and », which varies with 
the difference, - — on an curve. It is, therefore, ‘reasonable to expand 
the solution of Eq. 10 into the form of a polynomial in [&,(n) — 4, such ol 


in which 0 < 
a of two, rather than a higher degree, is 


— 

ie 


1981 


malay based on the number of available boundary ne of the ee 
and secondary flow needed to determine the coefficients, aed the desirability 


ofasimple formula, 
At the water surface on an curve, € = so that = 
? ae wind effect is neglected the shear st: dives on the water surface can be = 
assumed to be equal to zero everywhere, even at the banks, for a trapezoidal 
channel; a rectangular channel tends to have nonzero boundary shear at intersec- 
tions of two banks with the water “surface (6). Alluvial channels have | cross -_ 
: _ sections which in general are approximately trapezoidal and, thus, the assumption 
: of zero shear stress everywhere on the water surface mane be reasonable 
if the wind effect is negligible. The assumption implies a, = 0 everywhere. 
If the wind effect is not negligible, A o Will ill represent ‘the nonzero shear stress 
— 
With the boundary conditions 2 = 0 0, along with the derivative 
of t,, with respect to € on the water pas which can be obtained by 
differentiating Eq. a and letting € = &,(n), Eqs. 6 and 10 (or 4) give V, on 


(19) 


and F, defined by Eq. 
At the point (x, = D,x, = 0) where the water surface intersects with the 
X-axis, V, = V, = 0. With this boundary condition, Eq. 19 gives 
_ Eq. 19 also gives V, and, thus, explains the mechanism that generates he 7 
transverse flow on the water surface of an open c channel. The transverse flow _ 
V, is zero at the x x,-axis in which A: = 0; it also vanishes at the banks where - 
€ = & = 0. Between the x,-axis and the two banks the transverse flow on 
the water surface occurs whens magnitude depends on all factors shown by 
Eq. 19 while the direction depends on the signs, positive or negative, of A 


and In addition, has the positive sign on the left side of x axis, 


_is directed towards the x x-axis (the center) from the banks ¢ on ‘the both | sides, 
_ as often seen on the water surface of open channel flow. 
On the channel bed (bottom and walls) in which = 0, i.e. 
flow bw velocity to | the is zero, , so that Eq. 10 becomes 


With + Te ‘from Eq. and with respect to Eq. 22 


a 
| | 
| 
7 


SECONDARY FLOW 

in which Fy, h., h, and its derivative are evaluated at each pi point (€,,) on 
- the channel bed, so that a, is a function only y. In computing a, and a, 
by Eqs. 21 and 23, the x, gradient of normal stress (d0,)/(dx,) will be needed 
_ which may be approximated as equal to the gradient of hydrostatic pressure 
unit area , —(aP)/(ax, ), which arises owing to to the variation of channel 
cross section in the x, direction, a as other terms are - usually ‘relatively small. 
; It may even be approximated to be zero for a faisty uniform, open — 
flow. For a steady uniform flow, Eq.21 becomes 2] ree 


that A= = pgS(h, D)=0 (which steadily d to from the 


bank towards the x,-axis), 


[Eof) —€ol 


hap 


_— 


_ all the coefficients determined, Eq. 17 can be used to | ania the — 

i shear stress in the flow : and along the channel bed, that includes among other 
things the effect of secondary flow. If the flow is steady and uniform, the 4 
P difference between the shear stress computed by Eq. 17 and the first term © 
on the right side of Eq. 13 should give the portion of shear stress owing to 
the momentum transfer by secondary flow. 
7 7 The mean boundary shear along t the wetted perimeter of a channel cross” 


Fig. 2 ‘shows the shear stress lincluding that in 
- flow and along the boundary, along with corresponding primary flow velocities 7 
(measured)] and the computed secondary flow directions in a cross 
Ss (Section No. 255) of Rio Grande Conveyance Channel, near Bernardo, — 
'N.M. (5). The shear distribution is expressed in terms of the ratio of local 
_ Shear to the mean boundary shear. The cross section generally has a trapezoidal — , 
shape, although the side slopes appear very steep because of two different 
scales used for the vertical and horizontal dimensions in the Fig. 2. The boundary | 
a shear distribution on the channel bed has two peaks near the two corners of 
: ‘the cross section, which look much sharper than actual because of the exaggerated | 
vertical ‘scale. There are two peaks of boundary shear on 
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is anual quite well by results of direct measurements in trapezoidal — 
(6) as shown in Fig. 3. . Fig. 3 also shows results from using two indirect methods 
of estimating the boundary shear (one by Preston tube, and the o other from 


FF 
--@--@-- Measured Chennel Bed 


No.=2SS (x=25S00 ft) 
= 0.2579 (1. 2618 


ft; Sec. No.=25S 


x, 

FIG. 2. iti Stress, Primary Flow Velocity, and Secondary Flow as ime. 

= Rio Grande Conveyance Channel (5); 


= 1,280 cfs; 1 cfs = 0.3 mm/s; 1 
= 0.3 m]: (a) Computed Distribution of Shear Stress in Flow and Along Boundary; 


wy Distribution of Primary Flow Velocity (Measured and Computed); and (c) Computed 


the distribution of primary flow velocity near the heaieeay. 3 It shows that 
without using the direct method of measurement, the peak value of boundary 
shear that tends to occur near the corner on the channel bed cannot be detected. 


method, using the of primary flow near the 


| 4 

q 


Bs shear decreasing ‘continuously from the center to the side wall. In 
alluvial channels, or even in laboratory flumes, as long : as the channel is not a 
smooth the boundary shear is not distributed as indicated (misled) by ‘such 
an indirect method. Direct measurements, even in laboratory flumes, is very 
- 7 difficult especially near the side walls; thus, enhancing the value of a mathematical 

: model of shear stress such as Eq. 17 and the associated computational technique. & 

_ Fig. 2(@) also shows the shear stress: distribution including only the first — 
of Eq. 13 which represents only the portion of shear stress owing, to the gravity 
Ss A comparison of shear stress distributions with and without secondary 

flow indicates that the effect of secondary currents" is very p pronounced near — 

the secondary flow velocity, V,, is not very great in in magnitude, generally 
med than 10% of cross- -sectional average of p primary flow v velocity, but its effect 


on shear stress near we we walls | is not small (not secondary). as reason 0 


section 


i} 

el 


Cross 


Direct Measurement 
Indirect estimate from 
velocity distribution 
Preston tube estimate 


3.—Distribution of Boundery Shear in Trapezoidal Channel (Ref. 6) 
is that, as shown by the or the right side of Eq. “2, 
ss contribution of secondary flow to the shear stress at a point is cumulative i . 
or net momentum transfer in the & direction, along an y curve between the 
point and the water surface, and that near the side walls the secondary flow 
is is mostly directed toward the channel t bed (including th the bottom and side walls), 
as shown in Fig. 2(c), so that. V, carries the negative sign in Eq. 13° and, ‘thus, 
its effect is to increase the chest, The magnitude of h,, one of the metric — 
coefficients, also has effect on the contribution of secondary flow to the shear 
ss shown by Eq. 13. It is relatively large near the side walls and approaches — 
Zero near the x,-axis. Near the x,-axis, both v, and h,, are very small so that 
the effect of secondary flow to the shear stress on or near the x-axis is relatively 
small. The shear stress shown in Fig. 2(a) also includes the effect of flow 
| being nonuniform, i.e., aV,/ax, 4 0 as indicated by Eq. 12 which shows that 
th _ shear stress could be affected by every term of F defined by Eq. 11 a 
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a. aes a study of bank erosion problem, an equation such as Eq. 17, rather 7 


than Eq. 16, should be used to determine the peak» values | of boundary : "shear 
(drag force) that tend to occur on ‘on the side walls and on the channel bottom — 
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near the side walls. Peak boundary shear along with the convection by secondary | 

_ flow near the side walls should be a major mechanism responsible for bank 

and bed erosion of alluvial channels. Figs. 4(a), and 46), sl show / computed 
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—_ total velocity, i.e. , the vector sum of primary y and secondary flow | velocities, 
; which is also expressed in percentage of cross-sectional average of primary 


ees shows the distribution of computed shear stress at three consecutive _ 

sections (Section Nos. 240, 250, and 260), facing downstream and separated : 
of These | Sar. | thon om tie 


. atio iGrevity 
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= 
= 0. lbs/ft? (0.4087 kg/m? the boundary 


= 


Sec. Ne. =250 (x#25000 ft; =1. 104) 
= 0.0886 ibs/ft? (0. 4335 


200 (xF26000 Ofte 
= = 0.0 0.0821 (0.4017 kg/m?) 
FIG. 5. 5.—Three- Dimensional Shear Stress Distribution [Rio | Grande | Conveyance Chan- Ks 


(5), ender a steady flow of 630 cfs (18 m*/s). In Fig. 5 the values of 7, and 


(ratio of mean boundary shear stress to pgRS) are shown: and the shear 


distributions are again expressed in terms of the ratio of local shear to the 
mean boundary shear. It is not surprising to see that the shear distribution, 4 
= like the primary flow velocity distribution, varies from cross : section to cross 
section. The peak values | of boundary ; shear on the channel bottom and side — 


secondary flow velociti ages of 
& 
- 


Woe 

1 MILE 


. 6.—Relative Locations of Cross Sections Studied [Rio Grande ——— 
near Bernardo, N.M. (5)] 


Approximeted Channel Crose- 
d Chennel Bed 
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mae = 630.0 cfe: = 1.070 
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FIG. 7.—Variation of Shear Stress with Discharge Rate [Section No. 240, Rio o Grande ; 
Channel (5); 1 cfs = 0.3 mm/s; 1ft= 0.3m) 
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‘é the geometrical shape of actually measured channel bottom is very similar to 
_ the distribution of boundary shear on the channel bottom. The two deepest | =] 
locations of channel bottom that occur near the side walls coincide with those 
of peak boundary shear. The two peak values of boundary shear (one on the 
fee bottom and one on the side wall) on the left side of x,-axis in each 
of these cross sections are e greater than corresponding values on on the on 
The shear stress and primary flow velocity distributions have an irregular and 
_ alternating pattern of asymmetry. These phenomena seem to be owing to: (Ql) 
The asymmetry of cross-sectional geometry; (2) a fundamental instability of © ee 7 
flow in a straight channel; and (3) the effect of an upstream bend as shown 7 
Fig. 6 which also’ shows and relative locations 


similar. Probably, the most ‘striking finding from observing Fig. 5 ‘is that 


"generally increases although the general pattern of shear distribution does not 
change appreciably, and the ratio of local (point) shear stress to the mean boundary 
shear may even decrease, as shown in Fig. 7. For example, the peak boundary 

Shear on the left side increased from 0.417 lb/ft (2. 04 kg/m’) to | 0. 499 lb/ft | 
(2.44 kg/m *) as the discharge rate increased (about doubled) from 630 cfs (18 
to 1,270 cfs (35.6 m*/s) although the ratio of boundary shear 


about doubled, as the i increase (doubling) i in discharge rate occurred. 


At A general differential equation that ‘that go governs srs the se sediment eae is: 
which C = the concentration of suspended sediment; t = time; ¢, ‘denotes 
the diffusion the xy direction; u, = = the x; component of 


e@.g., dso, under the effect of gravity. Both C and U (and its components) | 
are time- raveraged quantities \ without including instantaneous random fluctuations. 


Then, u, is equal to the primary-flow velocity, V,; u, is equal to the sum 


of the x, component of the secondary-flow-velocity, V,, and the settling velocity, 
_ —V,, of the representative sediment particle, or u, = V, — V,; and u, is 


equal to the x, component of the secondary-flow velocity, , V;. Fora may : 
and nearly uniform flow in which fot) eat 


HY7  SECONDAR | 
hear increaced from 0 0825 Ih/ft (0 4 1569 Ih / ft (0) 7647 kg/m’), 
, sum of the fluid velocity and the settling velocity of the representative particle, sg ; 
| 
| 
| 
| 


> 

_ umerically with boundary values of C. If the secondary aah, 

Y, 


in which y = the proportionality constant. It hes been ‘Teported (1,7,8) that 
1 y = | for fine sediment particles, e.g., smaller than 0.1 mm in diameter; and 
< coarse sediments. = 2 and 3, €, can from 


62) The 


2 n 


in which €, and = the and ‘components of the diffusion coefficient 
for sediment transport. The y component of momentum transfer coefficient, 
» in n Eq. 31 should be z zero as the n direction i: is perpendicular to the ahead 


is zero. The E component of €,, can n be defined =, 
by Eqs. 1 and 17, respectively. _A similar expression ‘of and consideration 
(€,,),, being proportional to dV, /an, e.g., to the Prandtl’s turbulent 


0. In using yg and 17 together in Eq. 34 to “compute €,,)¢> OF in "Eq. 
4 to compute V,, Eq. | should be considered only as an empirical equation , 
which is often found to be quite good from statistical analysis of actual data 4 
of primary flow velocity. If other equations of different forms are found from _ 

7 data analysis to be better, they should be used instead of Eq. 1. In any case, = : 
a it is not even implied here that Eq. 1 can be derived analytically from the A 7 
shear stress equation such as Eq. 17, although the Prandtl-von Karman universal _ 

_ velocity distribution formula, which has the eee form, was ee 
j te The evaluation of the effect of secondary currents on sediment concentration a 
can be accomplished essentially by comparing the sediment concentration — 
¢. estimated by solving the aforementioned equations with or without secondary © Pz 
currents. effect point to within a channel cross section, 


j 
| 


as the direction and of “secondary flow tend to vary vary from point to 
2 of Rio Grande Conveyance Channel near M., but this 
_ shows the suspended sediment concentration in the section measured by the _ 
United States *s Geological Survey (5) along” with the boundaries s of a region in 
“which ‘computations for sediment concentrations using Eqs. 29 and 30 were 
performed to study the effect of secondary flow. The size » and location of 
he! the study region were restricted and determined by the measured data available. 
The data had to be used as references to determine the accuracy of computed 
‘Tesults; some of the data were also needed as boundary values in numerical — an 
“computations. The measured primary flow velocity distribution along with 
shear stress and secondary flow properties, to the 
4 
which were used i in computing sediment concentration. “ee Shen 


Measured Sediment Crncentrations 


8.—Distribution of Sediment Concentration, in grams per liter [Rio 
Conveyance Channel, Now | Mexico: in K = = 0.4 4451; Q = 1,280 cfs; 1 cfs = 0.03 mm/s; — 
- While properties. of flow and channel are indicated in Figs. 2 and 4, suspended 
sediment particles in the cross section have the median size, d.., of about 
0. 12 mm, for which the settling velocity, V,, at the water temperature of 20° e ; 
is 1.1 cm/s (0.04 fps). The value of y in Eq. 31 was made equal to unity — 
in computations of sediment concentration, in absence of a precise estimate 
. of it and also because of the sediment particles being fairly small. “giles ie 
In Figs. 9(a) and 9(b) the values of sediment concentration inside the a 
_ region, computed with Eqs. 29 and 30, respectively, are compared with the 
measured. The difference between computed concentration in Fig. 9a) and 
that in Fig. 9(b) represents the effect of secondary flow, as the former included __ 
the secondary flow in computation while the latter did not. The sediment 
- concentrations on the boundaries of the study region : are not shown in Figs. 
9(a) and 9(b) since they were the measured values" (with some interpolation, 
_ between measuring points) used as the boundary conditions, as shown in Fig. 


8. In with Eq. 30, which does not include the 


| 


393 
a 
4 
@ 
| 
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by secondary flow, tend to overestimate the sediment concentration at a = 


Lines of Equal Sediment Concentration, Measured 
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FIG. 9.— of Sediment gallons per liter [Rio Grande 
Conveyance Channel Section No. 255 (5); « = 0.4451; Q = 1,280 cfs; 1 cfs = 0.3 
mm/s; 1 ft = 0.3 m]: (a) Computation Including Soeenany © Currents (Eq. 29); ene 


gives very good results. The vertical component of secondary flow velocity, 

: in the ‘Study region, which is located in the central portion of the cross 4 

section, is weak but consistently downward in the same direction as. gravity, 
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SECONDARY OW 
as shown in Fig. aa). ‘It means that more downward sediment transport than 
that by gravity alone exists in the region. This is more so at upper levels of 
the region where the downward Vy is stronger than those at lower levels. AS a 
_ shown i in Fig. 4(a), V. roaches = 
“near x,-axis; it even Seoemees upward (carrying the positive ve sign) at lower levels” 
of the region, relatively far from the x,-axis. Therefore, errors in computed a 
_ sediment ccncentrations owing to neglecting the vertical (downward) sediment 
transport by secondary flow are relatively large at upper levels, and seemingly 
small or negligible at lower levels. The relatively large discrepancy between 
the computed s sediment distribution and the measured at lower levels, shown 
in Fig. 9(b), is mainly owing to neglecting the sediment transport by secondary | 
“4 flow in the transverse direction. The sediment transport in the transverse direction 
in the study region is mainly accomplished by secondary flow, as the gradient — 
a of sediment concentration in the transverse direction is very small as shown 
y by Fig. 8 or 9 and, thus, the sediment transport by diffusion should also be 
very small. As shown by Fig. 4(5), at upper levels the transverse component 
of secondary flow carries sediment into the study region from sediment sources __ 
at the two sides. Therefore, the neglect of this flow component should result __ 
in an underestimation of sediment concentration at a point, while the neglect 

__ of the vertical (downward) Componant of secondary flow results in an overestima- — 
tion of sediment concentration. The exclusion of the components 


upper levels of flow ‘should not be too bad because of ‘the mu mutually a 
effect of sediment transport by the two components of secondary flow (although 
they may not cancel each other completely). On the other hand, at lower levels - 
the transverse component of secondary flow is quite strong, and transports 
g sediment ¢ out of | the study Tegion, towards the two ) sides. Therefore, the neglect — 


of secondary flow at lower levels is generally weak (quite negligible); and, 

_ thus, the neglect of this component cannot offset the effect of neglecting the — 
transverse flow component. The main portion of sediment transport in the 7 
longitudinal direction through a transverse cr cross section of an alluvial channel» 
a achieved through lower levels of the cross section because of the large s sediment 7 

_ concentration. Therefore, the role that the secondary flow plays in transporting ‘ 

= especially in the transverse direction, should not be ignored. Also, 7 : 
_ wherever primary flow isovels have significant curvatures, e.g., under side-wall ‘ 1 


of sediment concentration at lower levels. As stated cartier, the vertical component : 


py the role of diffusion in the transverse (x,-) direction (a component of of 
_ diffusion in the € direction) should be included in sediment transport s studies. es. 


‘Summary ano Concusion oe 
oti 
‘Under the modeling technique used explicit links have 
been established among interacting variables in open channel flows, such as 
the distribution of primary flow velocity, channel cross section, discharge rate, _ 
secondary flow, shear stress distribution, and sediment concentration. The = 
_— of numerical computations and simulation depends on that of primary 


_ flow velocity distribution used as the basis of forming the | coordinate o system. 


— 
| 
| 
| 
| 
| 
| 
> 


rece 


Eqs. 1 and 2, used to represent isovels and veleahty distribution of primary 


_ flow, were developed from statistical analysis of United States Geological Survey 
data from Rio Grande Conveyance Channel (5), and appear quite adequate. | 
However, some adjustments will be necessary to cope with primary flow velocity 
distributions significantly different from that indicated by the Rio Grande data. 
_ For instance, in many open channels, isovels of primary flow near the water 
_ surface tend to curve toward x,-axis (center) and the maximum velocity occurs — 
_ below the water surface. An approach to such a situation is to use a separate Bo 
coordinate system near the water surface. Regardless of adjustments that may 
he 


r pen stress in the flow and along the bottom and side walls of alluvial channels. 


A Amy waren of degree two, the formula has coefficients determined from a 
quation (Eq. 10) and boundary conditions of both the 
secondary flow and the shear stress. No data on shear is required. It is capable ; 
of including the effects of gravity, primary and secondary flows, and nonunifor- | 
mity of flow, etc. Such a formula is essential to hydraulic studies in alluvial — 
channels , especially near the corners and side walls where no other existing © 
formulas work. The shear formula A values of boundary" shear on 


_ The geometrical shape of boundary shear distribution computed along the 
: channel bottom is found to be similar to that of actual (measured) channel 
bottom; the deepest locations of channel bottom that tend to occur near the 
7 corners and side walls almost coincide with the locations of peak boundary 
shear, , indicating a direct correlation between the boundary shear and the erosion 
or scour on the channel bed. The distribution of boundary shear, as well as 
the pattern of secondary flow that moves into the corner from above and out 
- of the corner along the channel bed as shown in Fig. 2(c), is clearly a main 
: mechanism responsible for the erosion and scour of channel bed and banks. 


Such knowledge and the computational technique presented herein for computer — 
; simulation of secondary flow and boundary shear distribution should, therefore, 
be an effective tool for studying and selecting better designs for protecting 
channel beds and banks from the scour anderosion, 
‘The secondary flow affects the shear stress which in turn affects the diffusion 
process. Therefore, the secondary flow affects sediment transport both directly _ 
(convection) ; and indirectly (diffusion, through its interaction with the shear 
stress). At upper levels of the main region of a channel cross section, not 
including the near-walls area, the secondary flow transports sediment into the 
region from the sources on the two sides, and helps the gravity force to a 
sediment downward. At lower levels, it penny ‘sediment out of the region = 
in the transverse 
sediment transport is sdone by diffusion. The downward secondary flow component 
_ is much weaker than the transverse component. The main parts of mechanism 
governing the sediment concentration is, therefore, the gravity, diffusion, and 
_ secondary flow. The traditional analysis considers only the first two; therefore, 2 
the inclusion of the third component should receive an emphasis in future studies 7 


= main. problem with using Eq. 29 in computing the: 


il 


convection with secondary flow, is the pa for data on all boundaries of a 
study region. An accurate estimation of y in Eq. 31 for computing the diffusion — 
coefficient for sediment to be used in Eq. 29 also represents a problem to 
be selved in the fetuses, = 
ie At present, there are no directly measured data on shear stress in alluvial: 
channels. There are also Shortages of ‘weaaeay shear data on on laboratory flumes, 


= in the central portion of channel bed. ‘This indicates, therefore, a great pod 
for data collections in both the field and laboratories, to check and, if necessary, 
improve a shear stress formula such as Eq. 17. An improved shear stress formula 
can improve the accuracy of secondary flow computations, and accurate second- 
ary flow measurements can in turn check and improve a shear stress formula. Ss 
A close interaction and coordination between mathematical modeling and data a 
collection i is clearly essential inthe future, ated 
im C. H. Lin and G. F. Lin, Research Assistants, Dep artment of Civil 
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ApPENDIXx 
a The following symbols are used in this paper: 


| 
q 


= cross-s -sectional area; 
A, cross- -sectional area on one side of x,-axis; on thew 
B, = top width of cross section on one side of x-axis; 
= water depth along x,-axis at cross section; 
= median diameter of sediment particles; in 
mean elevation of channel bottom at x,; ai 
index, denoting either left or right side of x,-axis; 
total velocity vector of sediment particle; 
component of velocity of sediment particle; 
x component of flow velocity; 


and y components of secondary flow respectively; 
x, Cartesian coordinates in longitudinal, vertical and transverse 
a, = coefficients i in shear distribution formula; 
vy = = constant in velocity distribution aan" 
= = constant in logarithmic velocity distribution equation; 
€,,€, = and y- component of diffusion 
‘= momentum transfer coefficient; 
s(€ we = and components 
‘SOV 
= value of for which V, = = 0; 


= value of on along y-curve on water surface; 
4 density of water; 


= shear stress in x, ‘direction in plane perpendicular to € direction; = 
= mean boundary ‘shear at cross section; 


boundary shear at point on on channel bed at cross ‘spetion. a 


eqalyale ox 


| 
| 
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Fiow RESISTANCE IN 


open n channels i is a a problem of central importance in hydraulic ai and va i aoe nos 
engineering. Flow resistance or hydraulic roughness coefficients are used _ 
principally to derive depth-discharge relations necessary, e.g., in waterway design _ 
and in the computation of sediment transport rates. 

aa Leopold, et al. (30) classified the hydraulic roughness of natural channels — 

‘into three elements: (1) Skin; (2) internal distortion; and (3) spill “itiies. 

Skin friction is a function of the roughness of surficial material; internal distortion — 

: ‘sadebiine is caused by particular boundary features such as discrete boulders, 
bends, bars and bed undulations that create eddies, and secondary circulations, 
spill resistance is associated with flow accelerations and decelerations. The 

_ Relatiogshipe between these various forms of roughness and flow variables, 


are complex and poorly understood. Spill resistance is probably negligible in 
the more-or-less clean, straight, symmetrical, coarse gravel bed river reaches 
under fully turbulent, near uniform, subcritical water flows considered herein 
and is therefore notexamined further, = 
Bi sediment transport regimes are studied herein: (1) The bedload transport 
ate is vanishingly small and the channel boundary i is rigid; and (2) the bedload © Z 
rate is finite and the channel boundary deforms. With the former a 
regime, possibly the chief contributors to hydraulic roughness (for the given 
_reach specifications) are grain roughness and internal distortion or form resistance p : 
generated by relict or remnant bedforms of previous floods. (With the latter _ 
regime, reactivated or developed bedforms probably are responsible for 
of the hydraulic roughness.) In contrast to sand bed rivers, which can have _ 
a variety of bedforms (35), apparently only dunes or bars develop under subcritical F. 
flows on the channel bottoms of gravel bed rivers (13,14). Indeed, Hill (23) 
has observed the persistence of dune forms, moving downstream under super- 
critical flows, in a laboratory channel with a gravel bed; whereas for a io 
~ Scientist, Water and Soil Division, Ministry of Works and Development, P.O. Box 
(1479 Christchurch, New Zealand. 
-Note.—Discussion open until December 1, 1981. To extend the one month, 
a written request must be filed with the Manager of Technical and Professional Publications, 
_ ASCE. Manuscript was submitted for review for possible publication on July 23, 1980. 
_ This paper is part of the Journal of the Hydraulics Division, Proceedings of the American 
Society of Civil Engineers, ©ASCE, Vol. 107, No. HY7, July, 1981. ISSN 0044-— 
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TABLE 1.—New Zealand Gravel Bed 


Reach cubic meters | in meters | width, B, 
per second in meters 
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FLOW 


perimeter, | radius, R, | Y,in | surface | rial, dso, | T, in degrees” 
P,in meters | in meter s | in meters 


0.00865 
0.00125 
0.00490 
0.0110 
0.0036 
0.00035 
0.00042 
0.00140 
0.00016 
0.00017 
0.00102 
0.00039 
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bed under these conditions, antidunes moving upstream re occur (28). a. 
» The purpose of this study was to derive gravel bed river flow resistance 
prediction equations suitable for use in engineering design and fluvial studies. 
_ Theoretical considerations are employed to assess functional relationships estab-— 

: lished by a dimensional analysis which is calibrated empirically by linear, multiple 

“regression. The large field data set listed is also useful for research and practical - 

studies of gravel bed river hydraulics and sediment transport. 
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0.00468 
0.35 | 0.00380 
0.76 =| 0.000458 
273 | 0.00200 
0.00121 | 
0.00530 
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0.00099 
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0.000399 

0.000338 
0.000268 | 0.038 
0. 000314 0.038 
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0.00038 0.076 
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= Data CoLtection AND SPECIFICATION ous 
_ The 136 field data sets used in this study were obtained from nN am . 
4% 46 New Zealand gravel bed rivers (Tables 1 and 2); and were collected 


‘mainly by the Ministry of Works and Development and Catchment Authorities 
‘as part of their routine water flow gaging programs. The following selection 


7 criteria were adopted in the compilation of the data set: (1) Straight channel 7 


Continued 
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River Dats for Mobile Bed Channels 
| depth, Water perature, 
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Bare 


Wairoa a 
_ alinement and a more-or-less constant width; (2) little vegetation on the channel - 
oo and banks; (3) stable banks and bed devoid of major irregularities; and 
relatively wide having a simple geometry that contains 


the “Weisbach friction factor, f defined 


in which g = the gr gravitational acceleration; R=t e hydraulic radius; S = 
the average longitudinal water surface slope, to average energy 
. slope; and V = the mean velocity of flow. Direct measurements and soundings 
were utilized to determine waters surface width, B, wetted perimeter, P, maximum — 


depth, y, and cross section area, A. From this information hydraulic radius, 


ond Benson (8); and\ water er temperature, T, was measured with a 1 standard n mercury 7 
thermometer. All the hydraulic variables were determined at a single cross section — 

7 during conditions of reasonably uniform flow, i.e., at or near flood crests - 
at low flows. The Froude Number (V/V gR) was always less than 7 0 and 
therefore flows were subcritical. Median size of surface bedmaterial, was 
determined using the ‘method of Wolman (39) or variants thereof (27), or by . 
sieving of surface material at reaches with fine gravels. __ 


Frow Resistance in Ricio Bounpary Channers i 
Lae 


— 


The criterion employed to distinguish rigid boundary (vanishingly small rates" 
of bedload transport) from mobile boundary conditions was Shields Entrainment 
Function, which for fully turbulent flows i by 


| 
K, was compute and in conjunction with current meter measurements oO point 
velocities mean velocity. Vand water discharce O were deduced _ a 
| 
- — 


6 | 17 | 238 | 0.00652 
Row 20, | «(1.65 | 0.00646 


in which S,= = specific gravity of | of bed particles (S, = 2.65 throughout t this paper); “ 
u,= = average shear velocity; an and v = . kinematic viscosity of water. em: 


Henderson 2) Shields’ original value of F, = 0.056 was ‘With 


1.60 x 10°°m’/sEqs.2and3 yield | 
50 


which is by data of Tables 1 and 2. Thus, boundary 


conditions vasa — if, from Eq. 2 and sac aforementioned _ 


Data satisfying this criterion arelistedin Table], | 


a Dimensional Analysis. —Since only water discharges in clean, ait uniform = 
channels are being considered it can be assumed that friction factors for subcritical — 
flows, satisfying Eq. 5, are a function of V, R, P, y, dso, g, p = density 


or equivalently 


The independent variables, P, y are included because, with the exception 
a of mirror images, the set R, P, y uniquely defines cross section shape (21). 
Because of a lack of other bedmaterial data, it is assumed that the surface - 
: 7 bedmaterial o or rbed pavementis a adequately described by the median grain d diameter, a 


| [based on intermediate particle axis (31)]. Choice of median, as opposed ; 


‘to other values in the bed pavement particle size distribution, is supported _ 
_by Bray (4) who showed that no significant difference obtains, with best fit 
Pa _ equations for roughness estimation, using several other measures of particle 


| 038 
7 
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; 
Th heoretical Considerations. —Keulegan (29) the Prandtl -von Karman 


flow i in | rough, , rigid bed channels: 4 


_ in which E = a coefficient; k = von i alate s constant of turbulent = 
and k, is the + et — of the bed surface. Eq. 8 may b be written (22 


) 

& 


In 


in which a = 
K 0.400, then c = 2.03 and, a, varies wit "The parameter, 

TABLE 3.—Correlation Matrix for Logarithms a of f Dimensionless Variables pares. 


Perimeters Vas /v | V/VgR V/V 

j 

Vdso/v 
VR/v_ 
V/VgR 88 | -0.027° 

0. 3 0.631 | 0.025 


, is on the shape of the channe an 
Koulegen’ s analysis (29) that for c = 2.03 


"provided flow width/depth ratios are not very small. For the data of Table 

1, in which the cross sections are reasonably trapezoidal and mostly very wide, — 


H 22) plotted funct Eq. 10) yields 
ey’s plotted function ( hon (Eq. ) an 


cross section yn shape o on friction factors is minor. The pate of suds ame 


- Apelt (25) also indicates similar independence for the data of Table 1. 
_ Moreover the Boundary Reynolds Number (VR/v) and the Froude Maniies. 
( riven )may be neglected because sdh are re likely tc to be unimportant: 


; ai can be ‘discounted (21). It may be ‘inferred from the preceding ‘that Eq 


et 
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= RESISTANCE 


for friction Quetere of reaches and hydraulic conditions specified i in this study. | 


computed for the data of Table 1. Where water temperature ‘had ‘not been 
measured, é an average value, determined ed during « other flow gaugings ; at the relevant 
reach, was used. As the dimensionless \ variables were found to be approximately 
log-normally distributed, statistical analyses were made with their logarithms. 
In order to examine intercorrelations a correlation matrix of their logarithms _ 

_ (Table 3) was computed; and then a linear stepwise regression technique was . 
applied in which the intercorrelations evident it in Table 3 were eliminated as a) 


TABLE 4. Friction Factor Relationships Based on Power Form of Equation 
nel “ 
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state 


=1. 125 (R/dy)” 


corp. 


1/Vf = 0.886 (P/ vied 

1/Vf = 1.33 

=: 2.22(VR/v)°"* 

V/V 2.21 (V/V 9d) 0.384 | 

1/Vf = 2.17 pee 37 | 0.216] 0.074] 13. 
= 0.321 q 


far as (9). A level of significance corresponding to a Students t value 


"explained 71% of friction factor variance inherent in n Table 
relative roughness (R/d,,) explained 59%; Froude Number (V/V gR), 9%; = 
nwod factor (P/R), 3%. Eq. 12, and indeed all equations in Table 4 involve | 
a degree of spurious correlation (3) in that all dimensionless variables contain — 


5 
eter- | error 
Equation ber | tion | units | F-ratio 
q Rigid 12 | 0.710] 0.0933] 656 
Dad 
Mobile 
- (n=52) 
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"Preceding co considerations, the fact that the a 


ersu 


is involved in ‘the regression, ‘and ‘the sa ai’ degree of variance explanation, the 
Froude Number (V/V gR ) was discounted from Eq. 12. The e regression equation 
for the two remaining dimensionless variables (Eq. 14, Table 4) shows that 
only a marginal improvement in logarithmic standard error results between Eq. 
which is solely a function of relative Toughness (R/ds0). Also 
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Friction factor. f 
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FIG. 1.—Friction Factor vs Relative Roughness for Rigid Boundary Conditions in 


Coarse Gravel Bed ‘Rivers 
‘the s spurious correlation between friction factor parameter, a7) al ‘shape 


factor (P/R) may contribute significantly to the 3% of variance explained by ; 
wg On these grounds the accepted result of the statistical analysis is. (Fig. 


in with wif 


Using data from 67 Canadian bed river Bray (4) obtained 


significant (r? 0.472, = = 60) statistical 


- = 0.0519 - 


| 
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. But since 


geen: (1) The statistics of Eq. (Table 4); the high correlation between 
Boundary Reynolds Number (VR/v); and relative roughness (R/ds.) (Table 
4); and (3) the low « "correlations between Particle Reynolds Number (Vd,, /v) 
ta both friction factor parameter (1/V /) and relative roughness (R/d..); a 
conclusion for the data of Table | is that the statistical ee of Eq. 


Resistance Equation For Ricio Graver Beo CHANNELS oH 
“ie The results of the foregoing analysis indicate | that, forthe given n reach ‘conditions 
or selection of descriptive parameters, the friction factor is given by (Eq. 11) 4 


in agreement with enstier work with gravel rivers by Leopold, etal., (30), Kellerhals 
(26), Limerinos (31), and Griffiths (19) amongst others. Contributions to the 
scatter of Fig. | are thought to result from: (1) Inadequate deecription of relev aoa 


in data collection ‘procedures; 3) in channel alinements and 
sections; and (4) rugged bed topography. Of these the last is qualitatively — 
: considered as the most pronounced for the New Zealand data (Table 1). At 
- most reaches, residual or relict bars were present to provide significant iatemnal 
distortion resistance. 
: In the estimation of friction factors for river _ works, th ne 
‘theoretically based, _ empirically calibrated equations (26), (30), and (31) of 
_ Keulegan form (Eq. 9) give reasonable predictions for ee regular gravel © 
ek : recommended, in ‘addition to the equations previously referenced, for use in i] 
‘ 2 natural channel design; instead the use of a more comprehensive — 
calibrated by 186 data sets including data from Table | and elsewhere (1, 
‘11, 24, and 40), is advocated. The equation is (19) 
whic Is R/ dg : < 200. For low Values of slative 


= 0.760 + 1.98 log io 


f. which should be compared with Eq. 16 (Table 4). Bray states that his finding " 
‘= at variance with the expected result (/vVf a function of R/d,, only) based : 3 
on widely accepted ideas related to the simplified case of two-dimensional flow _ 
‘ 


Lal 


Mobile boundary conditions, occurring when the — transport | rate is” 

dx. < 


4 4, 


y Data satisfying Eq. 20 are listed a Table 2. ‘Most of the Gate | sets correspond ‘* 


i, to vigorous bed conditions as as the n mean value of the Froude. Number (V/V gR) 


is 0.606 with a standard deviation of 0.154. Except at Shotover River reaches — 
Table 2), dune bedforms were not observed directly as suspended fos 
almost invariably clouds the water at high flow rates; but residual or relict — 
a bedforms are ‘commonly apparent subsequent to flood conditions. A | qualitative on 

- description of various bedforms | in a New Zealand gravel bed river is given 


om Veiga da Cunha ® analyzed mobile bed data —_ the fine gravel bed Mondego 


cand divided the friction factor, into two 


in which hf’ = friction factor due to grain roughness; end | f’ = friction factor : 


~ due to form roughness produced by bedforms and other channel irregularities. — 


This separation (Eq. 21) may be of conceptual use; nevertheless it is arbitrary. 


7 Moreover there is doubt (38) concerning the ability | of the plane bed resistance — 
Py relations, used to evaluate f’, to predict the grain roughness of a gravel bed hen 
a covered with undulations. Laboratory studies with coarse bedmaterial (32,34) _ Ay 
4 have generally conformed Cunha’s results. Furthermore, Gladki (15) has shown 7 
that the flow resistance relations derived by the previously mentioned workers ; 
(7, 32, and 34) fit data from the coarse gravel bed Raba River, moderately well. 

& 2 Theoretical Considerations.—Einstein and Barbarossa (10) argued that form 

, resistance contribution to the friction factor must be dependent on the ey | 


of the bed configuration which, in | turn, is a function of the sediment transport : 


terms it can be argued that for a mobile gravel bed with developed bedforms, jf 
“flow resistance is mainly due to form drag with “minor contributions , given * 


‘the high Boundary Reynolds Numbers (VR/v) ‘observed i in natural rivers, from 
‘grain skin the profile drag, F,,, of a bed undulation 


i. which C, = coefficient of profile drag. _ 
_ The average resistive stress, t, is from Eq. ‘ad 


| 
| J 
— 
| 
| 


now RESISTANCE 


But for the conditions in this p paper we have (20) 


In experiments with sand- coated two dimensional triangular roughness elements, 
simulating ripples and dunes, Vittal, et al., (38) deduced an empirical cnet 


be foregoing, C,~ y? Ce and Eqs. 26 and 27 see 


Interestingly Engelund (12) found, using a pipe expansion 


R 


= 


in which 1, = average bed shear stress at critical conditions of bedload Cees 


31 “applies to sand dunes, bu but the range of (38) ‘commensurate. 
with that of Table 2 (i7< R/ dsp = < 247). wrt 
Assuming that 30 and 31 apply to gravel dunes or bars, then 
30, and 31 give ‘ 
=fni 
with J F, constant (Eq. 2). eat 


: 
, 
| lags | 


has demonstrated however, for ‘sediment of specific 


oo in this study as average | shear velocity, u u,, may be rather insensitive. 7 

and mean velocity, V, is more accurately determined than | u, in ‘the field. 
These several arguments, together with Eqs. 33 and 34 and the neglect of grain 

roughness effects (Eq. 23), that is, R/d.. is assumed to be eannpastent, — 


for friction factors of wide reaches with mobile _— -s Eq. 35 takes implicit | 
account of grain roughness ef! effects, given Eq. 34 ar and the relative roughness 
term (R/d;s) in E in 22. 


TABLE Matrix for Logarithms of Dimensionless Variables 28 Computed 
from Mobile Boundary Data of Table 2 _—s 

om 


Perimeters 


Statistical Analysis. —The statistical analysis detailed previously was repeated _ 
for the data of Table 2. A correlation matrix for the seven dimensionless variables 
_ (Eqs. 6 and 7) is given in Table 5 which should be compared with the markedly 
different Table 3. Linear stepwise regression indicated that only the mobility 


parameter, vi the result was (Eq. 36, Table 4, 
4 
“which is in accord Eq. 35 and supports earlier wor work 15, 32, ah 
34) based on the Einstein and Barbarossa (10) parameter, | (Eq. 22). Note that a 
_ the relative roughness relation (Eq. 37, Table 4) is different from that for rigid 
_ boundary channels (Eq. 15, Table 4). The low value of the exponent of Eq. 
37 together with its statistics (Table 4) indicate that the friction factor parameter 


(a /Vf), is only slightly dependent upon relative roughness (R/d.,.). Comparison 
‘Statistical data 4 for Eqs. and 38 supports this contention. 


| 
0.273 | -0.095 | -0.466| 0.063 | 1.00 4 
VR/v_ 0.497 0.558 | 0.786 | -0.057 0.704 
| 0.242 | -0.304 0.025] -0.014 | 0.544 | 0.197 | 1000 
v/V gdw | 0.620 | 0.714 | -0.505| -0.126 | 0.306 | 0.660 | 0374 | 1.00 
| 


FLOW 


in 
; ‘the friction factor is given by (Eq. 35): f= jn (VIN Bd), and the ‘regression = 


result (Eq. 36) is i/Vf = = 2.21 (V/V gd,.)°™, will be less useful, for flow — q 
a resistance ‘Prediction it in | gravel bed | rivers, than it its rigid ren counterparts = 


= 35). In particular, hysteresis effects caused mainly by dune development 


persisting on the falling stage of 2 a flood, together with the ‘dependence of dune 


size on flood wave fe steepness (16 and 1) are likely to Wioay important scatter 


the channel boundary is rigid, the resistance to fully | turbulent flow 
D.. in ‘straight, ‘regular reaches of coarse gravel bed rivers is dependent, largely, z 
on the relative roughness, provided that bed and bank roughness are not 
distinguished and surface bedmaterial is described ibyt the median size. A statistical 


Resistance Equation Tor MoO avel De annels.— The foregoing analysis 
Ba 
| 
m*eq 
| 
AG. 2.—Friction 
Gravel Bed Rivers = 
q 
1 


3 (= for sdf tad? 


peeing nearly 60% of the friction factor variance inherent in a large field 

data set from New Zealand gravel bed rivers, 
cc _A theoretically based equation, calibrated by a more comprehensive data — 


 - set (including that of Conclusion 1), applicable to rigid beds and defined Sante “a 
— = 0.760 0 + 1.98 108 


_ where the mean velocity. of fl flow is is determined f from 


* consequently recommended for use in design, gn, particularly in ‘Teaches with 


arugged bed topography, 
3. For channels with active bedload transport, when the is 


"specification « of ‘maialy ona ‘a mobility parameter. Statistical model of this 


4. With the New Zealand data given ‘herein, the flow resistance matte of 
Conclusions | and 3 are useful for reconnaissance work, that is for coarse 


prediction. This concineten probably also holds for all models of the form comel an 
5. Other apart from measurement errors, are thought to 
for the residual variance in the statistical models. The main contributor in the 


rigid boundary state is probably form roughness provided by relict or remnant 


bedforms; and for mobile beds, an inadequate account of bedform euany 


Acxwowteooments 


The writer recognises and appreciates the to 


‘difficult and hazardous conditions. For assistance with the compilation of this 
paper the writer thanks A. Alderidge, Freestone, ‘Dz D. Hewson, M. ‘Hicks, ith 


Barnes, H.H., H. H., “Roughness of Natural Channels,’ Water Supply Pa, 
a I 1849, United States Geological Survey, Washington, D.C., 1967. bow gn 


2 C.. “Flow Resistance of Large-Scale Roughness,”’ Journal of the 


&g 
| 
| 
= 
pricy, ror mis puldanice 


be Division, ASCE Vol. 104, No. HY12, Proc. Paper 14239, ee, 1978 . PP. 1587- 1603. re 
Benson, M. A., ~Spetions Correlation in Hydraulics and Hydrology,”’ Journal 
r _ the Hydraulics Division, ASCE, Vol. 91, No. HY4, Proc. Paper 4393, July, 1965, 
4. Bray, D. ‘Estimating “Average Velocity in Gravel-Bed Rivers,” Journal of the 
Hydraulics Division, Vol. 105, No. HY9, Proc. 14810, Sept., 1979, 
5. Buchanan, T. J., and Ww. P., Discharge Measurements at Gauging Stations,” 
_ Techniques of Water Resources Investigations, Book 3, Chapter All, United States — 
Geological Survey, Washington, D.C., 1969. 
6. Burkham, D. E., and Dawdy, D. R., “‘Resistance Equation for Alluvial Channel — 
_ Flow,” Journal of the Hydraulics Division, ASCE, Vol. 102, No. HY10, Proc. Paper 
. Cunha, Veiga da L., ‘‘About the Roughness i in Alluvial Channels with Comparatively _ 
Coarse Bed Materials,”’ Proceedings, Twelfth Congress of the International Association 
for Hydraulic Research, Vol. 1, Paper No. A10, Sept., 1967, pp. 76-84. 
8. Dalrymple, T., and Benson, M. A., , ‘‘Measurement of Peak Discharge by the Slope 


_ Area Method,” Techniques of Water Resources Investigations, Book 3, Chapter A2, 
United States Geological Survey, Washington, D.C., 1967. 
9. Draper, N. R., and Smith, H., Applied — Analysis, John Wiley and Sons, 
10. Einstein, H. A., and Barbarossa, N. L., ‘ aie Channel Roughness,”” Transactions, — 
ASCE, Vol. 117, Paper No. 2528, 1952, pp. 1121-1146. 
‘ll. Emmett, W. W., and Seitz, H. R., a 


in the Snake and Clearwater Rivers in the Vicinity of Lewiston, Idaho,” Basic- Data 
Report, United States Geological Survey, Boise, Idaho, 1974. 
12. Engelund, F., “Hydraulic Resistance of Alluvial Streams,’ Journal of the Hydraulics 
ASCE, Vol. 92, No. HY2, Proc. Paper 4739, Mar., 1966, pp. 315-327. 
a Fahnestock, R. K., and Bradley, W. C., ‘“‘Knik and Matanuska Rivers, Alaska: a 
_ Contrast in Braiding,”’ Fluvial Geomorphology, M. Morisawa, ed., State University on 
of New York, New York, 1973, pp. 221-250. 
14. Galay, V. J., “‘Observed Forms of Bed Roughness in an Unstable Gravel River,” 
_ Proceedings, Twelfth Congress of the International Association for Hydraulic Research, ‘ 
Gladki, H., ‘“‘Resistance to Flow in Alluvial Channels with Coarse Bed Materials, ta 
Journal of Hydraulic Research, Vol. 17, No. 2, 1979, pp. 121-128. ae 
16. Griffiths, G. A., ‘‘Transport of Bedload by Translation Waves in an Alluvial Channe ; 
_ Research Report 76-1, University of Canterbury, Christchurch, New Zealand, 1976. 
wv. Griffiths, G. A., and Sutherland, A. J., “ Bedload Transport by Translation Waves,’ 
Journal of the Hydraulics Division, ASCE, Vol. 103, No. HY11, Proc. Paper 13363, 
ay 18. Griffiths, G. A., ‘“‘Recent Sedimentation History of the Waimakariri ines New € 
_ Zealand,’’ Journal of Hydrology, New Zealand, Vol. 18, No. 1, 1979, PP. 6-28. 
19. Griffiths, G. A., ‘“‘Rigid Boundary Flow Resistance of Gravel Rivers,’ Report WS 


ae 127, Ministry of Works and Development, Christchurch, New Zealand, 1979. ice o 
, F. M., Open Channel Flow, The MacMillan Company, Inc., New York, _ 

. Hey, R. D., ‘“‘Determinate Hydraulic Geometry of River Channels,”’ Journal of the a 

_ Bydreniics Division, ASCE, Vol. 104, HY6, Proc. Paper 13830, June, 1978, pp. 869-885. : 

. Hey, R. D., “Flow Resistance in Gravel-Bed Rivers,”’ Journal of the Hydraulics — 
Division, ASCE, Vol. 105, No. HY9, Apr., 1979, pp. 365-379. 
23. Hill, I. K., ‘‘Fluvial Sediment Transport at a Large Bed Shear Stress,” ’ thesis presented 

to the University of Canterbury, at Christchurch, New Zealand, in 1967, ba partial 

fulfillment of the requirements for the degree of Doctor of Philosophy. eee ee 

. Judd, H. E., and Peterson, D. F .» “Hydraulics of Large Bed Element Channels, 

Report No. PRWG 17-6, Utah Water Research Laboratory, Utah State University, 
A. K., and Apelt, C. J., “Shape Effects on Resistance to Uniform ¥ 
Flow in Open -Channels,”” Journal Hydraulic 17, No. 2, 1979, pp. 


| 
4 


Kellerhals, R., Channels with Gravel- Paved Bed, Journal of the Waterways 
and Harbors Division, ASCE, Vol. 93, No. WWI, Proc. Paper 5091, Feb., 
= n. Kellerhals, R., and Bray, D. I., “Sampling Procedures for Coarse Fluvial Sediments,’ ; 
28. Kennedy, J. F., ‘‘Hydraulic Relations for Alluvial Streams,’’ Sedimentation Engineering, 
an Vanoni, ed., ASCE, New York, N.Y., 1975, pp. oo" Seo 
_ 29. Keulegan, G. H., ‘‘Laws of Turbulent Flow in Open Channels,”’ Journal of Research ‘ah 
of the National Bureau of Standards, 21, Research Paper 1151, ‘Dec., 1938, 
30. Leopold, L. B., Wolman, M. re Miller, J. P., Fluvial Processes in Geomorphology, 
W.H. Freeman and Co., Inc., San Francisco, Calif, 1964. 
31. Limerinos, J. T., 
in Natural Channels,” Water ‘ater Supply Paper | 1898-B, | ‘United States Geological = 
; 32. Raudkivi, A. J., ‘Analyses of Resistance in Fluvial Channels,”’ Journal of the Hydraulics 
Division, ASCE, Vol. 93, No. HY5, Proc. Paper 5426, Sept., 1967, pp. 73-84. > 
_ 33. Rouse, H., “*Critical Analysis of Open Channel Resistance,” Journal of the Hydraulics _ 
Division, ASCE, Vol. 91, No. HY4, Proc. Paper 4387, July, 1965, pp. 1-25. 
ss Shen, H. W., ‘Development of Bed Roughness in Alluvial Channels,’’ Journal of 
the Hydraulics Division, ASCE, Vol. 88, No. HY3, Proc. Paper 3113, May, 1962, 
4 35. Simons, D. B., and Richardson, E. V., “Resistance to Flow in Alluvial Channels,” 
Professional Paper 422-J, United States Geological Survey, Washington, D.C., 1966. 
7 36. Simons, D. B., Khalid, S., Al-Shaik-Ali, and Ruh-Ming Li., ‘‘Flow Resistance in 
w Cobble and Boulder Riverbeds, ** Journal of the Hydraulics Division, ASCE, Vol. 
105, No. HY5, Proc. Paper 14576, June, 1979, pp.477-488. 
37. Thompson, S. M., and Campbell, P. . “‘Hydraulics of a Large Channel Paved 
with Boulders,” | Journal of Hydraulic Research, Vol. 17, No. 4, 1979, pp. 341 -354. 
38. Vittal, N., Ranga Raju, K. G., and Garde, R. J., “‘Resistance of Two Dimension! 
“J Triangular Roughness,”” Journal of Hydraulic Research, Vol. 15, No. 1, 1977, pp. 
39. Wolman, M. G., “‘A Method of Sampling Coarse River-Bed Material,” Transactions, — 
"American Geophysical Union, Vol. 35, No. 6, Part 1, Dec., 1954, pp. 951-956. | eal 
40. Wolman, M. G., “The Natural Channel of Brandywine Creek,”’ Professional Paper 
7 Yalin, M. S., and Emin, K. ., “Steepness of Sedimentary Dunes,”’ Journal of the 


_ Hydraulics Division, ASCE, Vol. 105, » No. HY4, Proc. Paper 14502, Apr., 1979, pp. 


vA 


hee 
| 


_ STANDARDS FOR COMPUTER-BASED DESIGN Srupies 
William James,’ M. ASCE and | Mark A. Robinson? 
ott 24 Yo ytilideilss od) lo yideid 
_ This paper examines some special requirements of computer-based studies. a 
These requirements are supplemental t to those Felating to the specifics of a 
particular study. The purpose of these a additional ‘requirements is to guarantee 
_ credibility and confidence in the computer modeling. It is suggested that these 
- special computer-modeling requirements could be set out in the initial study 
specifications or request for proposals. If these requirements are excluded 
therefrom, it is suggested that consultants submitting study proposals recommend 
, inclusion of these items to their clients. The suggested additional specifications © 
summarized at the end of this paper. gt | seals 
= In acquiring design services from a reputable engineer, the client expects — 
that the study will be carried out properly, as indeed is almost always the 
+ case. The activities listed surtenggli then, apply more to those cases in which 


level of internal system detail. An example is the Stormwater Management _ 
Model (SWMM) (7) used for urban hydrology and drainage hydraulics. Its — 
lies in the ability to take into account the interaction of a reasonable number 
of processes and to report on the status of certain important variables within 
these Processes throughout the simulation. Of course, there would be little point © 


therefore a appears | to > be ‘that many, ‘if not all, of the eleven , processes. and 
i interactions can be sufficiently accurately investigated, at a large number 
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American Society of Civil Engineers, ©ASCE, Vol. N No. HY7, July, 1981. ISSN_ 
-796X / 81 /0007-0919 / $01.00. 


iy 
f 
| 
| 
| 
| 
| 
| _ with modeling will normally carry out these activities automatically. The ideas — 
of are not entirely new; ina sense they have been applied to undergraduate sical ; 
are canable of describing civil engineering svstems to ga high | 
iil USINE Such il the syste Wad CMUUEM 
_modeled with the use of less expensive calculating machines than the main — 


by the engineers responsible for the design. _ Thus, programs rec requiring n main net 
support should provide a thorough evaluation of reasonable alternative rel 


a A difficulty that arises is that the person responsible for devising and coding © 
_ the original program is usually not available for interpreting the performance 
of the program in every application. Also, ‘the engineers “responsible for the 
computer simulation, often on a considerably modified version of the o: original 
program, are frequently remote from the interpretation of the design and ultimate 
presentation of the results to the public. Whereas the original author of the 7 
program and the engineer responsible for submitting the computer runs may — 
be highly confident of the reliability of his particular computer program, the 7 
client, e.g., local authority or municipal engineers, and the people responsible — _ 
for budgeting and ‘supervising the ultimate project, may ‘not have the same 
confidence in the computer methods. Indeed, they may be wise to be sceptical 
until the results are provento be sound. 
_ This paper attempts to establish controls that ensure that the best model 
- ‘ne program are used correctly. Many of these control measures coincide wih 
guidelines for using computer programs — by the Association of Profes- 


that “credibility and confidence in ‘the “results are created and 
~ sequentially from the design office to the client, e. g., from the junior engineers 
to the project engineer, perhaps in the consultant’s office, and thereafter to 
the municipal engineers responsible for supervising the study, and ultimately yp 
SS the engineering committees and other political bodies to the beneficiaries, 
frequently the public. The p process is depicted schematically in Fig. 1. In other 
D words, it is not sufficient that the computer model represents all relevant physical = ‘ 
processes of the study accurately enough, or that the results are correct or 
- sufficiently accurate. It is important to assure that the study be carried out — 
in such a way that there is little chance of using a wrong model, or using 
wrong data, or wrongly interpreting the results, or simply of not understanding 
model or its relationship to the design objectives. 
a The writers argue herein that certain precautions be taken at the outset of 
It is recommended that the study specifications include, among other 
things, a minimum number of standard tests for model verification, sensitivity 
More ‘Comptex Computer-Baseo Desicn esicn ENVIRONMENT 4 


Several factors have increased the complexity of computer- based urban 
_— drainage sates over the years: ( More and more programs are becoming 


e 
"hardware is increasing: (4) the cost of computing - is decreasing: (5) software 
_ capabilities, especially of the smaller computers, are becoming more sophisticated; 
(6) computer communication between so offices and remote main frames 


920 
| 
| 
| } 
| 
| 
» improve their understanding of computer modeling, despite the longer learning 
7 q times associated with this methodology, and higher salaries and related costs 
L (and a tendency for some experts not to divulge the mechanics of key packages). = 


STANDARDS. 
‘These developments increase the pressure on regulatory agencies to approve 


_ sophisticated design technology which they may not have had time to study — 
_ fully themselves. Some of these factors are further described subsequently. _ 
he Computer Models = * is evident that the next decade will see an explosion — 
4 in software development, especially applications programs. Recent studies (3,6, 8), = 
. = , list several hundred computer programs currently available for solving 
problems associated with water resources development. In fact, it seems that 
every university and consulting engineering office is bent on developing their _ 
a own computer programs for the solution of civil engineering problems. Moreover, = 
while most of the popular programs are being enhanced to ‘better: describe 


constituent processes, or to include more process, other programs describe some 


‘COMPUTER CENTER: | COMMITTEE: an 


| 


| ENGINEERS: progress reports 
repeated 


«of these constituent processes in greater details. Thus, the problem not now 
Gg arises is the selection of the correct sequence of systems and process models" 
to be used in a study. Similarly, if a program adopted is very complex, it 
_ may not be clear which processes may be safely ignored, or that the data — 
used is the best of many Possibilities. At some point i in the ‘study i it may 


ot prefer s studies based on the same model. Problems such as as these which occur 


| 
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Detter tO USE SEQUENCE OF Process MOUEIS tO One OF The 
models taking into account similar processes, but doing so inaccurately. Evidently 
| 


in the mur area were reviewed in a recent paper (12). 
4 Most tried-and-tested programs are considerably enhanced each year. The 
latest versions of SWMM, the Illinois Urban Drainage Area Simulation (14), | 
and the Hydrological Simulation Program—-FORTRAN (11) may operate in . 
continuous mode, including water quality Processes. Since these programs have 
been under development for 10 yr or more, the precise stage of the evolution 
of the program must be correctly identified when the program is used. Program — 
documentation seldom matches the capabilities of the program at any point 
in time; consequently it becomes important that the actual program —— 
are accurately and carefully recorded. In their guidelines for standards of use 
for computer programs, the APEO recommends appropriate d documentation - of 
versions, theory, and software cs: These should be met, as a minimum, > 
organization distributing a program for general use. penal 
_ Many consultants and network computer centers maintain their own versions — 
of large program packages such as SWMM. Often these versions incorporate 
4 substantial modifications, usually aimed at core reduction, improved turnaround, 
_ simplified data preparation, or more relevant output interpretation. Again proper 
-Hardware.—Readers will be aware of recent advances in programmable 
calculators (5), and the rapid evolution of microcomputers with their bewildering 
range of configurations and capabilities. Similar situations exist within the range 


of minioomputers and main frames and ‘the machines that lie between | all these 


Increasingly, in order t to offer a large amount nt of computing for available 
“computing dollars, it is becoming necessary to distribute the computational effort 
- involved i in any one study between local microcomputers or minicomputers and 
folio main frames. Already, many programs have been modified so that their 
{ 


= 


solution methods are better suited to the large, cheap memories and small, 
bit words available on the small systems (13,15). Again, there are few gi guidelines * 
for this distribution of computational effort. 
Systems Ware.—New systems software products include local disc- -operating. 
"systems, local graphics, and sophisticated word processors. These products 
P considerably aid both output interpretation and input data editing. These, in 
turn, have benefitted from enhanced communication on speeds (e. dial- “up speeds 
to 1,200 BAUD). In addition, 
hee and international computer networks “available, so that ‘the smaller design 
offices may easily access the largest computers. These trends help smaller . 
- consulting engineering offices to gain remote job entry into the large programs ~ 
. maintained by various vendors. With the skills of a recent graduate or post- 
graduate e engineer ar and an inexpensive terminal, any office has has it immediate access 


without the burden of maintaining or updating the packages. 


anda ources. 
~ Consultants usually 1 review the complete study problem and available = 
‘resources at the outset in order to establish a mutually agreed upon scope > 
and basis for the model —" The —— review should not only be done 


| | 
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all those characteristics of the study problem that will aid in the selection of — . 
computer programs. | . It should attempt to list specific problems that ought tc to” 


be considered in any model study undertaken. 1. If hydraulic jumps or energy 


dissipation at manholes are important in a drainage s study, e.g., or if diversion 


structures are present, these should be identified. Few, if any, programs — 


_ adequately describe these phenomena. 
‘The study objectives and relevant objective functions should be listed and 
Benny All questions to be answered must be written down in order to explain : 
how the numbers generated by the « computer model will relate to the general 
3 objectives of the study. Urban runoff computer models, e.g., often produce 
_ only hydrographs and pollutographs resulting from a certain limited hydrologic | 
- time series. Simple parameters based on these response functions (hydrographs) 
include: peak values, time-to-peak, and total runoff or loadings. On the other 
_ hand, the study objectives often require a least-cost alternative to a specific 
flooding or drainage problem, such as the production of toxic or other harmful 
pollutants as a result of a certain hydrometeorological time series. This concise | 
explanation of the logical relationship between the model study 1 results, and = _ 
_ the objectives of the design is important for the rest of the study. ore — 
_ To many observers, the model produces only a sequence of numbers. It 
od is advisable to describe early in the study how these numbers may best be — 
compared b between several computer runs. It is difficult, e. Ss , to compare whole — 
_ response surfaces and curves such as hydrographs. In other words, the perfor- 
mance criteria for assessing the various design proposals must be carefully — 
_ explained. In some cases it may be sufficient to estimate the total annual loadings, 
whereas in another case the peak response is more important; perhaps a correlation 
eS whole response surfaces is necessary. The criteria for the oniemene 


be settled initially, a point established and agreed on by both the ‘client nt and 

_ Accuracy is of overriding importance. If the required accuracy can be | 
established early in the study, it would make the selection of the models and © 
_the level of discretization for the model study an easier task. In a sense, , the 

accuracy level predetermines the | type e of model to ‘be u used. There seems to 

be a mismatch between the nny of the results demanded from the computer 


ok one design against another ; should not be vague or appear to \ vary. It should 


systematic and the random sources of error should be investigated for both ~ 
the model (input data, concepts, solution methods, - computational accuracy, 
ete.) and the field equipment. There is little point in producing simulation r results” 


_ for validation that are many times more accurate than the observations which 
are used to validate the model. Consider, e.g., the accuracy of rain-gage sampling 
__ of moving storms, and the whole question of design storms. In any case, high a 
accuracy may not be necessary in view of the often tenuous relationship between " 
: the study objectives and the performance criteria. ‘It would be useful t to » 
automatically plot the spread of results rather than a thin’ response line such — 
as a hydrograph (or a smudge instead of a pollutograph)! tot 
_ Model Selection Criteria.—There is much published literature and data on — 
4 the structure and performance of various models. In water resources, see, e.g., 
Refs. 4, 6, and 8. a review of the do documentation w will indicate 


HYD STANDARDS 
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= 
| 
— 
| 
1 
qa 


the accuracy with which selected, required important processes are modeled. 
; In their published guidelines (2), the APEO has stressed the fact that program ; 


; - documentation should provide this information to any prospective users. The oi 
engineers conducting the study should attempt to show why selected models __ 
& y P 


Teview and special hydrologic and hydraulic processes required for their particular — 
study. The purpose here is to show how models may be used to support one 

another, or to establish a reasonable sequence of model use in the study. No 
- data and available study resources including the models, then establish the criteria 

for selecting the models and thereafter, select the relevant models. Concurrence 


ent 

for the design study is a central point of this paper. Adherence to this fundamental 
7 “approach to model selection should naturally bring to light all the points considered 

~ herein as guidelines. Assuming that the design engineer has the specific skills 
 (e.g., in hydraulics, hydrology, water quality) these familiar, basic, and funda- 


It is much simpler to seek out and review the available data than to collect — 
_ mew data from scratch. In two Canadian studies currently in progress, e.g., 


Srp years of sewer water level records, as well as rainfall records were _ 


- surprising finds. Such data could form part of the model selection criteria; 7 
- to make maximum use of available data the study should not be constrained 
to use a model unsuited to these (unknown) existing data. Water levels, e.g., 
_ could be directly computed and plotted rather than tediously converted > 
The study resources are limited to the available manpower, time, and money; 
the management of a study is an attempt to maximize the level of detail subject - 
to these manpower, time, and money constraints. Collection of field data is 
a primary drain on these resources. The field program must be dovetailed to +=§ 
complement the modeling effort. The key is sensitivity analysis. ==” 
Model selection criteria should be clearly stated in the study report. It is 
easy to criticize a study at a later date for not using a model which would © 
 deactibe processes to a much higher level of system detail if it is not made | 
_ Clear that the simpler model was used because of too little expertise, time, 
_ Or money to meet the effort for higher accuracy or detail. Other criteria include 
; d availability of the model, availability of supporting advice, experience with similar — 
or alternative models, availability of data, special costs related to proprietary 
versions of the models, and the professional interests and jealousies of the 
_ engineers involved in the mechanics of the program. This could mitigate against 
the use of proprietary programs, forinstance, 


produce sufficient information to answer a wide range of 


; are deficient or sufficiently accurate for their purposes in the light of problem 2 
| 
between client and designer should be sought stage. Finally specify 
© Pertinence listed inthis paper. 
= A certain amount of time should be invested in searching out existing data. _ 
| 
Sters to Estasusn Moon 
- series of verification, validation, and sensitivity tests can be designed to -— 
questions about the 


The tests should also to detect certain types o of error. 
_ Verification Tests. —Verification tests use some specific conditions for t which © 2 
* model response can be a to chank if indeed the model is 


7 rather comparisons between model | Fesponses and theoretically anticipated results 


Of course it will not normally be necessary for the design group to verify” 
_the same model in the same way for every design study. Summary information| 
will normally be sufficient for subsequent studies involving essentially the same 
model and the same processes. But it must be done at least once on receipt 
of the model ahead of the first design application. 
_ It would be useful to create a standard input file for verification tests. For . 
stormwater models, e.g., a hypothetical system comprising two simple, square ra 
subcatchments, of e.g., 1 acre each would be suitable. The subcatchments could 
be joined by a pipe of standard diameter and simple form. The hydrograph 7 
_ from the first subcatchment is attenuated in the pipe. At the downstream end 
this p pipe the hydrographs from both subcatchments are re superimposed. At 
the outlet of the pipe, the combined hydrograph is routed through a simple, — 3 
standard, perhaps rectangular, storage tank. The purpose of such a verification 
_ data set would be simply to test the algorithms for: (1) The generation of the 7 7 


overland flow hydrograph; (2) routing in the pipes; (3) the addition of two 


hydrographs; and (4) storage routing in storage tanks. Moft 
3 As a further example, we have used the following tests in our stormwater i) 
- runoff programs: (1) Zero rain to ensure that no runoff is generated; (2) very 
_ steep catchments to ensure that the hydrographs generated are very similar : 
to the input hyetograph; (3) light rain and high infiltration rates on low percent : 
impervious areas to ensure that no runoff hydrographs are generated; (4) 
completely impervious ‘catchments to ensure that» the total volume of runoff 
is equal to the total volume of rainfall; (5) low continuous infiltration rates 
and high continuous rainfall to ensure that the correct amount of infiltration — 
is subtracted; (6) very flat pipe gradients to check the surcharge calculation; 
(7) similar tests with small-diameter pipes and high pipe roughness; (8) high _ 
and low initial abstractions; and (9) tests on the storage- routing parameters, 
namely, the outlet rating curve and the storage curve. 
In all verification tests particular attention is paid to the quantitative summary | f 
output: e.g., total precipitation, total gutter flow, total snowmelt, total infiltration, 
ete . This is an essential check that the numbers generated by the computer 


carries out a series of tests. In our instrectional applications, e.g., 
undergraduate classes, the following verification data are normally used er 
s any stormwater drainage analysis: (1) Low imperviousness—5%; (2) high 

(3) high detention storage— —0. 84 4 in. (15.7 mm—46.7 


= _ model’s performance. Evidently, such tests will not materially increase the _ 
| 
7 
a 
| 
7 . are sensible. Ihe user must check the numbers by hand and ea = 
a sppeecciiall algorithm nnd data ffile should be built into the program. 
| 


impervious area with zero detention (High), a as a a percentage 5%; (6) om antiga 
—10 
in. (254 mm/h); (8) low « continous: ‘rain—0. in./h (2.54 mm/h); high 
continuous infiltration—2 in. /h (50.8 mm/h); and (10) low continuous infiltra- 
7 _ Calibration and Validation Tests.—Calibration implies the comparison of model 
> simulation results to field measurements, to another model known to be accurate, , 7 
or to some other adequate criteria to ensure | that the model of the system 
is producing accurate data. If these co comparisons indicate that the model results - 
are not sufficiently accurate, then the model of the system is altered, usually — 
by adjusting one or more program parameters, and the procedure is repented. 
_ This process generally involves several iterations before a satisfactory confidence 
, level is achieved. Techniques used in calibration include: (1) Comparison of 


“results against field observations; (2) cross-correlation of continuous model results 
with those of another proved, usually discrete event, or process model equivalently 
7 initialized, and (3) some combination of field observations and modeling. Wee? 
_ Validation implies testing the model of the system with a data set not t used 
in the calibration procedure. The most accurate method of validation is the 
_ comparison of output from the calibrated model against corresponding sets of — 
a The calibration and validation process should be limited to the most important * 
or sensitive parameters, and to the range of parameter values applicable to — 
“the normal operating conditions of the system. First, acceptable tolerances must 7 _ 
be established. They should be related to achievable field observations and 


the accuracy of ‘the field equipment. ‘Then, _ emphasis should be placed ‘upon 


‘the system. Reasonable assumptions may be satisfactory for less critical parame- 
ters. Thus, calibration tests are closely related to sensitivity analysis, examined 


4 Once verified, calibrated, and validated, the | model can applied 


input quantities that cannot be analytically measured may be considered to 
be empirical parameters. Following this approach, important empirical parameters 
_- discrete event urban rainfall-runoff models are: (1) Width of subcatchments; 
@ ground slope; (3) initial infiltration rate; (4) final infiltration rate; (5) Manning’s- 
roughness for overland flow; : (6) depression storage for both pervious and 
a impervious areas; and (7) infiltration decay rate. 
} Whether or not an input quantity should be considered an empirical parameter, 
initial condition, or a field-determined NSE, depends on the user’s 
opinions of the basic modeling concepts. For example, the SWMM models 
subspaces asa sasingle rectangular plane, forthe purpose of convertingahyetograph 
into a hydrograph of overland flow. Lumping the hydrological characteristics _ 
of all surfaces in the real subspace into an equivalent rectangle is widely known " 
as “parameterization. * The equivalent width of the rectangle cannot be deter- 
- mined analytically in the field. It could be considered a quasi-empirical eubopecs 
"parameter. Again, another user may not support the basic Horton- infiltration — 
approach | used in the SWMM and consider such such infiltration cid at toast for 


for combined quantity and quality “modeling. (10). Spatially- -averaged 


- 


Level of Discretization.- —The p procedure. for systematic 
_ described i in an earlier paper (9). Disaggregation implies modeling more subspaces 
of smaller size using a finer time step. For coupling the time increment io 
the size of the ‘Subspatial elements, the concept of an impulse response function 
| useful, e.g., an instantaneous unit hydrograph. It may be accurate enough 
for our purposes to represent this response function by at time © vector of 20 
i elements. _ This effectively sets the time step. to 1/20 of the time base of the 
in impulse response functions. Sensitivity tests in which the time step is systemati- _ 4 
_= cally changed may also be appropriate. In practice, disaggregation and sensitivity 
analysis proceed simultaneously (9). Often, careful disaggregation and appropriate 
selection of the time step will produce better results than extensive — 
of empirical factors such as Manning’s roughness for pipes and overland flow. _ 


A recent paper | provides some guidelines(1), 
_ There is an obvious conflict in the selection of the level of discretization. e 
Consultants’ costs increase rather rapidly with an increasing number of subspaces. - 
‘ More data have to be abstracted and prepared, and more expensive computer 
runs will result. On the other hand, the client will usually prefer to have a 4 : 
= high level of discretization for the agreed upon fixed study price. It is useful 
to have, at the outset of a study, an indication of the desirable level. As a . 
- guide, it is important to identify all significant elements in the system. These — 
_ represent locations for which it is necessary to generate a response surface, 
such as hydrographs, e.g., at all diversion structures and outfalls in hydrologic 
_ Studies. If these elements can be identified at the time the: terms of reference 7 
are set out, consultants will have a a better idea of t the scope of the work. It 
‘makes sense to face this possible conflict squarely z at t the beginning | of the study. 
aw Sensitivity Analysis.—Sensitivity analysis proceeds by holding all parameters 
but one constant at their expected values, and perturbing that parameter within 
_ reasonable expected limits such that the variation of the objective function 
ms can be examined. If apparently small perturbations of the parameter produce — 


’ changes in the objective function, ‘the system is said to be sensitive to 4 
_ that parameter. The user must obtain a measure of how accurate that parameter 7 
om must be represented in his model. If the objective function is not sensitive 
_to the pertubated parameter, then the parameter need not be accurately represent- — 
4 If the system is insensitive to the pertubated parameter, the parameter — 


: and its 3 associated process is redundant and the he process | should be deleted. The 


ta the actual values of the constant parameters may affect the sensitivity — 
analysis, and so, their values should be typical of the conditions being modeled. 
_ Here again, algorithms should be available that permit the user to easily conduct a 
4 sensitivity analysis. Another special button within the program is desirable. 
When this ‘‘auto-sensitivity’’ button is pushed, the computer requests the user — 
to identify the parameter whose sensitivity is to be tested, and the range of ro) 
es values. The data file is then automatically rebuilt and the tests 
carried out. All output response functions, such as hydrographs, should be 
on the same family of curves in order to the impact immediately 
tig se function are peak response and 


discrete-event models, a re or empirical parameter to be __ 

4 

— 
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“total mene wa'e as the pe peak | flow and n runoff unoff volume ina runoff 
Following the general approach to parameter determination 
earlier, the runoff parameters which affect each of these components are 4 
: Peak Flow: slopes (pipe and land), Manning’s n, subcatchment width, 
_ 2. B: Volume: imperviousness, as a percentage, infiltration, surface retention, — 


Itis critically i important to rank environmental parameters that affect the important 

components of your response function for your model in this way. *\ entiiogine 4 
Control of Errors.—Having carried out the required verification, calibration, — 
and analyses at \ various ‘levels of discretization, ‘the results 


a the tests have been achieved. The author of the report must explain or 
interpret the results carefully for readers who cannot submit the computer runs. f 
The verification tests must be shown to produce expected results. The calibration. 
and validation ‘Tesults must be shown | to accurate. 


evidence that the model is indeed performing properly: (1) The verification 
tests demonstrate that there are no serious errors in the coding of the model, ; 
(2) the calibration and validation tests at a simple level of discretization, say _ 
one or two subcatchments and pipes, demonstrate that the model is being used 
in a reasonable way. way, (3) the validation tests on the full data set, i.e., for all : 
of the subcatchments acting together as a complete system, indicate whether 
serious blunders have been made in preparing the input data set, and (4) the 
sensitivity analyses will indicate whether the level of effort put into estimating 
the individual parameters is appropriate based on their significance in affecting 
% The required output interpretation mn provides: the argument that proves that 
the results in fact make sense, and helps to ensure that the study authors have 
not misinterpreted the model results. All of these points are summarized in 
he 
= This paper suggests that the following topics should be included in any > 
computer- -based study. It is appropriate to specify them in the initial study — 
3 Problem Review. —The problem re review should identify all significant elements — 


processes, and to be disaggregated to the correct degree. 
_ Study Objectives.—The study objectives should be reviewed to show a 
the objective functions, such as pollutographs, or hydrographs, or both, relate | 
< to the design alternatives. New problems may | become ¢ apparent during th the study 


| 
| 


80 it must be possible to redefine ey objectives and all subsequent « deductions. ’ 
_ Performance Criteria.—The performance criteria for the | comparison 1 of one 
design alternative against another include peak response, total loadings, time 
_ to peak, cross-correlations, etc., and must be correctly identified so that the 
simplest modeling can be justifiably used to select the best design alternative. i 
Requisite Accuracy.—Accuracy of field measurements for validation should 
A be carefully reviewed in order to ensure that the model runs are not inordinately — 
expensive. Systematic and random errors should be defined. 
_ Review of Available Programs.—Several programs should be s suggested o 
“sect for review. The review should consider process models, as well as 
system models, and an appropriate sequence of models, 
Study Resources. —Study resources include time, manpower, money, and 
“information, and these, in turn, will determine which of the models may | be 
Model Selection Criteria.—All available data must be reviewed. Models should , 
be selected where expertise is available and data easily 
the limited study resourcesavailable, 
a Model Verification. Verification tests should be required on a specified simple — 


data set consisting of simplified elements « chosen to test each p process of interest. 
_ The verification tests should check the summary output and demonstrate that 


on one of the subspatial data sets specified parameters to demonstrate 
that the model is being correctly used. Validation tests should be carried out 
onan independent data set to demonstrate that the input data are reasonable. _ 
Level of Discretization.—The smallest number of subspaces required 
for modeling the system should be , selected commensurate with the objectives _ 
of achieving the best design at a reasonable cost. These minimum levels should 
correspond to the disaggregation necessary to identify response surfaces at all 
significant elements in in the system, which should be specified at the “outset, 
Sensitivity Analysis. .—Sensitivity analyses should be carried out on a pn 
number of parameters, e.g., infiltration parameters, roughness values, widths 
_ of subcatchments, etc., to identify which are of most significance, and to justify 
the effort put into their estimation. malt 2) 
x Data Preparation and Output Interpretation.. —All input and output s should be — 
re _ interpreted to demonstrate that the model is performing in a logical way. The - 
magnitude and direction of errors should be estimated. ‘WAS: 
_ Documentation.—The version of the program actually used in the study should 
be identified and appropriate documentation sources listed in the report. In 


addition, the machine readable agers and output files should be listed and archived — 


The delivery of computer-bas based engineering. be improved 
by enlightened study specifications and attentive control by the responsible — 

municipal or local authority engineers. The suggestions in this paper should = 

reduce the chances of inexperienced engineers using an inappropriate model, 


inappropriately u using an suitable model, through erroneous data 
P 


wien the are acceptable. Credibility would thus be 
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Journal or on any paper presented at any Specialty Conference or other meeting, the Proceedings irs 

of which have been published by ASCE. Discussion of a paper/technical note is open to 
anyone who has significant comments or questions regarding the content of the paper /technical a 

note. Discussions are accepted for a period of 4 months following the date of publication “Taw | 
of a paper/technical note and they should be sent to the Manager of Technical and Professional i 
Publications, ASCE, 345 East 47th Street, New York, N.Y. 10017. The discussion period may 

extended by a written request from adiscusser, = 4 
_ The original and three copies of the Discussion should be submitted on 8-1 /2-in. (220-mm) a 

ia by I1-in. (280-mm) white bond paper, typed double-spaced with wide margins. The length of a 

a Discussion is restricted to two Journal pages (about four typewritten double-spaced pages 
x of manuscript including figures and tables); the editors will delete matter extraneous to the =. 


- _— subject under discussion. If a Discussion is over two pages long it will be returned for shortening. i. ? 
 s Discussions will be reviewed by the editors and the Division’s or Council’s Publications —__ 
oo _ Committees. In some cases, Discussions will be returned to discussers for rewriting, or they 
_ may be encouraged to submit a paper or technical note rather thana Discussion. = = = 
. Standards for Discussions are the same as those for Proceedings Papers. A Discussion is LN 
‘subject to rejection if it contains matter readily found elsewhere, advocates special interests, 
is carelessly prepared, controverts established fact, is purely speculative, introduces personalities, 
or is foreign to the purposes of the Society. All Discussions should be written in the third 
person, and the discusser should use the term “‘the writer’’ when referring to himself. The — 
author of the original paper/technical note is referred to as “the author.” 
Discussions have a specific format. The title of the original paper/technical note appears 
at the top of the first page with a superscript that corresponds to a footnote indicating the ais” 
‘month, year, author(s), and number of the original paper/technical note. The discusser’s full KA ; 
name should be indicated below the title (see Discussions herein as an example) together with | 
his ASCE membership grade (if applicable). 
_ The discusser’s title, company affiliation, and business address should appear on the first — 
page of the manuscript, along with the Proceedings paper number of the original paper / technical 
_ note, the date and name of the Journal in which it appeared, and the original author’s name. - ae 
_ “4 Note that the discusser’s identification footnote should follow consecutively from the original cai 
- paper /technical note. If the paper/technical note under discussion contained footnote numbers 35 
1 and 2, the first Discussion would begin with footnote 3, and subsequent Discussions would 


Figures supplied by the discusser should be designated by letters, starting with A. This also 
applies separately to tables and references. In referring to a figure, table, or reference that =a 


Du Discussions may be submitted on any Proceedings paper or technical note published in any _ 
_ appeared in the original paper/technical note use the same number used in the ny 
ms It is suggested that potential discussers request a copy of the ASCE Authors’ Guide to 2 
; - the Publications of ASCE for more detailed information on preparation and submission of - 


or Waves By SHoRE-ConNECTED BREAKWATER" 


Closure by Volker W. Harms,” ASCE 


the to Memos his comments. In his 
discussion, Memos focuses principally upon the problem and probable cause 
of the wave- -height nonuniformity encountered. As described in the original 
7 paper, it had been found that regular waves of very-nearly smnacitel form = 
_ could be generated in deep water (d/L = 0.39 and H/L = 0.43), but that 
_ the resulting three-dimensional wave field was certainly not of uniform height. 
Wave heights varied both in the direction of wave motion (but not monotonically, 
ie., no _ evidence of significant viscous dissipation) and laterally along the | 
wave-crest line, but at every location the wave record was regular and did 
not vary significantly with tiene. If this in itself was not surprising, the degree 
of wave- height variability was. . For some wave conditions, the ratio of extreme ' 
wave heights, « , exceeded a value of 2.5; here ¢ is the ratio of largest to smallest 
wave height recorded within the wave-tank test area, i.e., at 520 points, with 
regular waves of constant height at each location (as in Fig. 6). This finding . 
. raised serious questions about the results of diffraction studies in which the 
7 incident wave height was measured at only a single fixed position along the 
_ wave-crest line; depending upon the spatial wave-height distribution that happens 7 
- to prevail (see, e.g., Fig. | ), this measured incident wave height could vary 
by a factor of 2.5 depending upon the position adi) 
sf Upon recognizing the seriousness of this problem, it was first ascertained 
b that it was not unique to our wave tank or wave generator. Thereupon, a 
 trial- rand-¢ -error Search for wave conditions H/L, at) that would 1 minimize 


sional wave field was undertaken recognizing that some wave are 
_ more likely to be troublesome than others, and therefore are to be avoided © 
q possible. Among those to be shunned were: (1) Wave-generator cut-off 
frequencies—to prevent ‘“‘wave-sloshing modes’’ (transverse waves), as predicted 
_ by linear wave-generator theory; (2) lateral basin-resonance conditions—to avoid » 
~ the generation of standing waves between basin walls; (3) substantial Benjamin- 
Feir-type amplitude modulations—to ¢ ensure that wave heights remain reasonably 
constant with time (e.g., € < 1. 1) at each wave-gage position, ne matter how 
g from the generator; (4) cross waves—since they are not only standing waves 
but are also of lesser frequency than the primary waves being generated; and 
i, (6) highly-reflective basin boundaries—since these reflect energy back into the 
test area and produce standing- s-wave effects throughout the basin, i.e., 4 


7 "heights may be locally constant but vary spatially. Recognizing the —— 


“December, 1979, by Volker W. Harms (Proc. Paper 15075), 


— 


= 
i such potential problem areas was certainly important but did not provide a 


@ The following criteria were finally chosen as the most ail 
and important for the selection process: (1) The cneyes time history at “" 


testing area and well alined when upon the beach. re 

- tees discusser concludes that irregularities in the wave generator contributed 
to the spatial wave-height nonuniformity, (as, e.g., stroke variations along the 
length of the paddle), as did reflections from the wave-tank boundaries which, 

¥ though perhaps very low, can never be totally eliminated in the continuous- 
_ operating mode. The writer shares this | conclusion, but would like to clarify — 
a misunderstanding that emerges from the discusser’ . arguments supporting this 
‘position. From the near-center line wave-height distribution shown in Fig. 3, 


any wave in the basin, 
it must not be inferred that the lowering of € was caused ie the lanes 

_ of the model (and might therefore lead one to conclude that strong teflections 2 
“from the main beach were being blocked). Instead, the difference in « values — 
i em should be attributed principally to the difference in sampling spaces: « = 1.6 
is the largest value that was encountered upon surveying the whole basin (wave 


_ Fecords at 520 points), whereas the lower \ values of € shown in Fig. 3 are 


for a much smaller data subset acquired near the basin center line (see Fig. ad 
_ The discusser also feels that basin-resonance conditions contributed to the — 
wave-height nonuniformity, particularly since the wave length of L = 0.69 m 
_ (that was finally adopted for the test program) is close to one of the natural 
modes of oscillation of the wave tank, albeit the 56th mode, with L = 0. 700 
m. The writer was never able to detect any influence of this ‘ “closeness” to 
a basin-resonant mode or to one of the wave-generator cut-off frequencies, 7 
and wonders whether such high basin-oscillation modes can actually be excited, _ 
_ and whether linear wave-generator theory can be meaningfully applied to such 
; wide basins (the wave maker extended across a tank 28 L wide). It should 
be noted that resonant wavelengths could hardly be avoided during our exploratory 7 
5 on since, for our test condition, they are so very-closely spaced, e.g., lateral 
_ basin-oscillation modes “paper 55, 56, 57 correspond for resonant wavelengths 
of L = 0.709 m, 0.700 m, and 0. 685 small frequency a 


laboratory experiments have shown them to be detectable (H/L > 0. ee 
they were present, then at it such | a low level that, in the opinion of t the writer, - 


the question becomes somewhat academic; serious cases of Benjamin- Feir 
‘instabilities could easily be avoided by simply ensuring that the wave record © 


| 
: ; th time; (2) wave attenuation should be negligible over the length of the a 
7 basin; (3) the variation in wave height across the tank should be as small as _ 
i 
7 _ which were obtained with a variety of models in the basin (and correspond — 
= 
ES 
‘The discusser considers Benjamin-Feir modulation instabilities (1) to have = 
. heen present because the adonted test conditions fall within range that m Ss : 


+ 
4 


for each gage on the 
1.1 (as in Figs. 6, 8, 9). When such wave-height stability i is evident at each 
wave gage, in our case over a distance of 21 m (70 ft), then the influence — 


contributor to the problem. Recent findings from research in the area of deep-water 
wave instabilities (28,29) indicate that the Benjamin-Feir analysis is - valid only a 
“for the initial growth period and that, in the long-time evolution, unstable 
- modulations typically grow to a vet subside, and then repeat themselves 
in a series of modulation and demodulation cycles, without reaching a led 


of complete breakdown of the wave train. The presence of this so-called 
Fermi- Pasta- Ulam | recurrence could not be detected at any time 


_ of Benjamin-Feir-type modulation instabilities can be ruled out as an a 


of Putnam and Arthur’s 1948 ward (16) which indicated that “diffraction oll 
is conservative,’’ and the present investigation, in which diffraction theory was . 
not found to be conservative; in fact, at large distances in the shadow zone, — 
measured wave heights: consistently ap theoretical values. This is indeed 


limited laboratory data base available to Putnam and Arthur at the time (some 
30 yr earlier!); wave heights were generally measured at only four locations i 
in the shadow zone, as opposed to 50 in our study, the incident wave was 
obtained from only one fixed wave Tecorder, measurements were not performed y 
as far from the structure, etc. This in no way detracts from Putnam and Arthur’ Ss 
very substantial contribution: that their results are still being used is tribute | 
7 in itself. Instead, the present results are intended to augment and complement | - 
_ their findings and in so doing, establish an experimental data base adequate 
— for a comprehensive understanding of the diffraction process. 
‘The discusser suggests that errors in the measurement of wave height may 


be responsible for the relatively high » values of K' at large. distances in the - 
‘Shadow zone (compared to theory); a given measurement error would 
7 a greater impact since wave heights are relatively small in that region. This — 
argument is initially appealing but, upon inspection, demands that distant wave a 
gages (see Figs. 7f-7h) have larger measurement errors than other gages, 
otherwise the agreement between measurement and theory that is evident close 
to the breakwater (see Figs. 7b-7e) even for s small K', could not be explained. 
In fact, all gages were similar in construction, had very nearly the same calibration © 


factors, and were similarly operated, and would therefore be subject to wave- 


height error bands of the same magnitude. The initial supposition is ont 
most probably it incorrect. ave pay 
_ The discusser has a apparently conducted diffraction experiments on shore-con- 
nected breakwaters under very similar conditions. The results from these studies 
will hopefully become available to us in the near future and perhaps shed more 
= on some of the questions that have been raised. ts 
With regards to experimental procedure, it should be clarified that the *' ‘I5- “y 
waiting period between transverses”’ (p. (1514) of the instrument carriage, , refers 
to the continuous data-acquisition mode and not, as understood by the discusser, 
to the intermittent mode. Each traverse of the instrument carriage- -generated © q 
amalog wave traces such as those shown for wave- 2, 4, 


and 6 in Fig. 9. WES 


| 


28. Yuen, H. C., Lake, B B. ‘Instabilities of Waves on Deep Water,” 
Review of Fluid Mechanics, Vol. 12, 1980, pp. 303-334. 
29. Lake, B. M., Yuen, H. C., Rungaldier, H., and Ferguson, W. E., “Nonlinear 
Deep-Water Waves: Theory and Experiment: Part 2, Evolution of a Continuous Wave 
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The writer thanks Gogiis, for his comments and interest in the paper. 
3 discusser has expressed two concerns, namely how is a evaluated and why are a 
inconsistent results obtained upon application of the writer’s model to Fischer’s 
test data of series 2800, 2900, and 31000. 
To determine the value of a, one may oiiiadl as indicated by the discusser. 3 
A slightly simpler method has been used by the writer: “a starts with 
Eq. 25 and requires that f(t) pat and df/dt = 0 at t = t,. The coefficient im 
a can then be determined by trial-and-error for different’ values of B, as 
summarized in Table 2 (for B = 0, a limit — was performed, resulting ne 
in a = 8 In 2 = 5.545). Table 2 shows that a = 5.55 for B =< 0.1. For B 
» 0. l, a increases slightly, being equal to 5.763 at B = 0. 5. When preparing .. 
the paper, the writer did not ‘suspect that values of B ‘exceeding 0.1 could 
be encountered (see Fig. 7), and thus suggested that a be a constant equal 
to 5.55. Recently, however, the writer found values of B as large as 0.28 for 
dispersion in a pool-fall sequence (the information is contained in an unpublished 
report by the writer). Where B exceeds 0.1, Table 2 should be taken into account. 
To find B from Eq. ats plot af versus B, based | on Table 2; enter graph at 
ap = (VdAT/dx)’ and ‘determine corresponding values’ of B. Apply a similar 
procedure when evaluating from Eq. 31. 
_ The fact that B may exceed 0.1 introduces an additional consideration. It 
_ may be recalled that Eq. 15 is an approximate version of Eq. 14. The term 
which has been neglected from Eq. 14 may become significant when B > 0.1 
: x/L is small. It can be shown from Eq. 14 that, for small x/L, the right-hand a 
side of Eq. 16 should be multiplied by V 1 +8 . This implies that, small 
values of x/L, the definition of a should read: = 


a January, 1980, by Spyridon Beltaos (Proc. Paper 15118), sts a 


* Research Scientist, Environmental Hydr. Section, Hydr. Div., National Water Resea 


~ Inst., Canada Centre for Inland Waters, P.O. Box 5050, Burlington, Ontario, Canada — 
4A6; formerly, Research Water Research Council, Canada. 
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wr 
a for which numerical values of a are given in Table 2. For 0.5 < x/L < 3.0, 
the ‘multiplier of (AT/G,)’ on the right-hand side of Eq. 39 will be less than 
a 1 + B, becoming ~1.0 when x/L = 3.0 even if B is as large as 0.5. When 
- B = 0.1, Eq. 39 reduces approximately to the definition given in Eq. 30, regardless _ 


0.005 

0200 


5.763 
= 
: rn comparison of AT with | G, using experimental results, the © writer 
“has carried out a comprehensive analysis which has resulted in the following: _ 
one AT and G, are related as predicted by Eq. 30 
2. For the tests analyzed in the paper, AT/s, varied from 0.60-2.11 and > 
had an average value of 1.64. The second lowest and highest values of AT/<, a 
were 1.08 and 2.06, respectively. The theoretical value of AT/s, is equal to 
5.55 55 = 2.36, which is some 44% larger than the ‘‘observed”’ average value. - 
In terms of (A7/G,)’, this deviation should b be doubled to 88%. This result 
may be interpreted as a flaw of the theory o or simply as a reflection of the 
error associated with * ‘observed’’ values of G,. The writer’s view coincides 
with the latter interpretation, based on his experience with &, computations 
from experimental C-t curves. (Note that, in addition to G, being very Sears: 
to the tail Peer of a Ct curve, as suggested by the discusser, the tail region — 
; the prevailing C are very small and usually | 
- comparable to the lower nanautiliies limit of the instrumentation used.) Additional 
= support for this interpretation is provided in the following Raion sana | 
3. According to Eqs. 28 and 7”, the value of ATC,/(M/Q) should be equal 
)/V For AT/&, = 2.36 (= V 5.55), this works out to be 0.94. 
**Observed”’ values of ATC, /(M/Q) were calculated by substituting \oCdt 
in place of M/Q ( (to account for tracer losses), and compared with both 0. 940 


| 

i 
—— 


_and the value of ‘‘observed”’ (AT/G, )/V 2m. “Observed” values « of ATC, 
dt ranged from 0.69-1.04 with an average a 0.86. This is much closer to 0.94 
‘(predicted by the theory) than to 0.65 which is obtained using the average 
value of AT/ca,(0.65 = 1. 64/V ). Moreover, of the 51 tests 


considered, only on three occasions did the latter T approach give better a > 


‘he po concern of the discusser is the failure of the model to produce — 


results when applied to Fischer’s test om: of series 2900, and 


pk 


FIG. " —Concentration-Time lids Observed at Different Lateral Positions; | 


Series 2800, by Fischer (13); Example Ref. 33: (a)x= 16. 16.07 m; 
x = 26.11 m (1 m = 3.28 ft) 
3100. The discusser states: ““Fischer’: s experiments (15) conducted i in laboratory 
‘fam should nearly satisfy the condition of one-dimensionality required for 
_ the application of the author’s model.’’ The writer has carefully examined the 


detailed results of these test series and found that he must disagree with the 
- discusser’ S premise (note that Ref. 15 gives only cross-sectional average parame- 


ters; detailed data are given in Ref. 13). Considering, e.g., series 2800, Fig. 


11 shows C-t curves (solid lines) observed at different lateral positions within 


the same cross section for two cross sections located at x, and x, (also see 
Table 1). The curves of Fig. il can hardly be be considered to indicate “one- dimen- - 


EO 
| 


‘these were considered sufficiently of two-dimensional 
“mixing and were thus utilized to test a pertinent numerical ‘simulation algorithm, - 
described in Ref. 33—the broken lines in Fig. 11 are simulated results. Waa) 
_ Similar considerations apply to series 2900 and 3100, which the discusser ; 
may wish to verify directly using Ref. 13. For the present closure, it is considered 


sufficient to provide Table 3, in which certain characteristics of the pertinent 


_C-t curves are summarized. It is seen that, at x x=) x,, there is a threefold cy 


variation of C,, (across the channel) and a fivefold variation of AT; these figures 
a are reduced to twofold and two and a half-fold, respectively at x = x,. It 
: -_ now be recalled that the writer’s definition of one-dimensionality calls 

“for nearly identical C-t curves across any one value of x (though shifted in 


time). This implies that, to apply the model, corresponding values of i and — 
AT should not vary excessively across the channel. It is therefore concluded — : 
that the data used by the discusser are not representative of one-dimensional 


TABLE 3.—Selected Dispersion Characteristics of Fischer's Series 2800, 2900, and 


swath 
Range of C, (arbitrary units) _ Range of AT (in seconds) 
indicated from C-t as Indicated from Ct 


id 
Curves at Six Lateral Positions at Six | Lateral Positions 


ys = 16.1 m, 16.1 m, and 17.5 m for series 2800, 2900, and 3100, respectively. a 
ght = 26.1 m, 26.1 m, and 29.5 m for series 2800, 2900, and 3100, I, wee | 


- Note: 1 m = 3.28 ft. ‘ 


dispersion. As such, these data cannot provide reliable estimates of B and = 


at 


Concerning the latter, it could safely be stated that it should exceed the value ~ 
of x,and L;shouldexceed3x,, 


7 & Errata. —The following corrections should be made to the original pape paper: 


pan: 


‘by small time interval” 


160, Eq. 17: Should read = =constant” 


instead of constant’ 


Page 161, line 1: | Should read “(as will be shown” instead be 

Page 165, line ‘Should read ‘ instead of “tp Wo 


‘= 
x=X,_ atx =x, atx=x,_ 
| 
{ 
| 
| 


- Pegs 165, paragraph 4, item 2, line 2: Should read “‘e, (= transverse mixing 


coefficient) aRV,” instead of RV, ‘Heat 
Page 166, Fig. 6(a), The label of the ordinate: Should read * u’ lilies rather 

33, Beltaos, S “Transverse Mixing in in Natural Streams,” Report SWE 78-1, Transportation 
and Surface Water Engineering Division, Alberta Research Council, Edmonton, Canada, gna 

one 
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OF Hy DROGRAPH DETERMINATION: 


los by Larry w. - Mays" ‘and Lyn Coles” 


‘Th ‘The comments, suggestions, and additions , presented by Ford, Meredith an and 


Byrnes, -Wormleaton, and Diskin are pertinent and of interest to the writers. 
betas commented that practical applications of the model | seemed unlikely because 
_ procedure Then Ford goes on to present a modified form of the model that 
includes an initial loss volume and a loss rate. The modified form of the model 
is then solved by geometric programming. Basically, all this does is to include 
- the simple operations that can be done by hand in a model that is much more 
- difficult to solve, especially for large problems, with several multiperiod storms. 
Ford did not state how baseflow separation was handled. The writers feel m 
Ford complicated a simple problem with no advantages presented at all. 
stated by Wormleaton, the computational simplicity of the linear cence : 
(LP) approach and the ease with which it can be applied make it a very arate 


- been used in the objective function. In fact, the first version n of this paper 


in which | to the weight th ordinate of the. ith ‘observed. 
- runoff hydrograph. The reviewers felt that the idea of weights should be deleted 
from the paper. Meredith’s suggestion of using the weighting factor, of Eq. 
‘January, 1980, by Larry W. Mays and Lynn Coles (Proc. Paper 15129). ba 
"Asst. Prof. of Civ. Engrg., Coll. of Engrg., Environmental Health and Water Resources 


Planning, Ernest Cockrell, Jr., Hall 8.600, Univ. of Texas, Austin, Tex. bak. 
*Planner, Texas Dept. of Water Resources, Austin, Tex. Beer 


| 
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DISCUSSION 
_ Meredith and ‘Byrnes, Wormleaton, and Diskin each commented on the choice 4 


‘constraints to relate runoff to rainfall in those storms in which M,< — ’ 
‘s Meredith and Byrnes suggested using an M somewhat larger than that indicated 
by the data results in the model determining the optimal number of unit hydrograph — 
ordinates, which is simply the number of nonzero unit hydrograph ordinates. | 
| This s also is a good suggestion to enhance the capability of the model. ae 
a, Wormleaton suggested that a disadvantage of the LP method is that an initial 
guess is required for the solution (i.e., unit hydrograph solution). This is not — 
; true because most LP codes do not require an initial basic feasible solution | 
= (initial guess); in 1 fact, only the more sophisticated codes can handle a user 
writers disagree with Womileaton's s comment that the most significant 
drawback of the LP method is that the number of decision variables increases | 
with the number of storm events considered. This is a disadvantage; however _ 
LP codes do exist that can handle hundreds and ¢ even thousands of constraints 


a 
a 


writers doubt that there are few, if any, ees that have more reliable 
_ existing data in excess of what can be handled by an existing LP method. 
Meredith and of the LP approach and the least- “Squares 


the least-s -squares models m may be more. accurate. However, ‘the writers 


"The writers do agree with Diskin that the LP approach to the suultistorm 
optimization has advantages i in comparison to the least-squares solution ie 


mam unit hydrograph for several multiperiod storms when the input dota 
are observed precipitation and total runoff hydrograph ordinates. Such a peroneal 
is presently being developed at the be pritiga of Texas by the first writer. 


nonlinear constraints. The model solution requires the use of the large-scale 


= reduced gradient technique by Lasdon and Waren a at 


—. Lasdon, L. S., and Waren, A. D., * ‘Generalized Reduced Gradient Software for 
Linearly and Nonlinearly Constrained Problems,’ Paper 77-85, 


L 
¥ 
__ Suggestions concerning reduction of the problem size are applicable. sj 
q 
excess and consequently, the rainfall losses are considered as unknowns. This a 
i 
a 
q | 


Closure by Robert M. Snyder,” M. ASCE 
i writer appreciates the recognition of Kundzewicz and looks forward to 
— With respect to the discusser’s issue of ‘ “‘intuition,”” the writer must point — 
out that while the “analogy may be proved by comparing the —— 
equations describing certain electrical and hydraulic systems,”’ the discovery, | 
_ if indeed there has been one, came not from excursions into pasigraphy but, 
_ rather, from a recollection of familiar phenomena. Arguing Cantorian antinomies 

(12), however, would more properly be carried out elsewhere. 4 

‘The writer agrees that flow transformation along an estuarine channel is more a 
properly represented by a distributed circuit analogy. This has been previously | 
- pointed out in detail in a personal communication between the writer and Price, 
So was mentioned oe by the \ writer in the full version of the paper as 


felt it best to. ‘introduce. the concepts of hydraulic capacitance, hydraulic 
inductance, and their related reactances in the simplified form of ordinary 

_ differential equations. This was considered important since most hydraulic » 
engineers are not familiar with radio electronics. As pointed out by Ref. 13, ma 
breaking with pure empiricism can lead to closer ties with physical analysis 


work is just a beginning of the physical analysis of flow based on the fundamental — 
equation of forced hydraulic motion of which ‘‘steady” or “‘de”’ flow is a — 
a special case. For further discussion of “hydraulic analogs, see Ref. Bead) v 


12. Poincaré, H., Science and. Method, Dover Publications, ine. , New York, N.Y., 1952, 

13. ‘History of Hydraulics,” Rouse, Inc., Institute of Hydraulic 

Research, The University of Iowa, lowa City, Iowa, 1957, 
-. Snyder, R. M., Christensen, B. A., and Walton, R., ‘Use of Electrical Ansiogs 


Asa Design Tool for Tidal Canals,”’ Proceedings, of the shen Stormwater Management 
In Coastal Areas, ASCE, June, 1980.00 


1980, by Robert M. Snyder Paper 15202). 


e which is still unlimited in its possibilities. It is hoped that the discusser’s referenced : .f 


pon 
5 
| 
| 
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Discussion by Richard H. French 


well maar ouber det: | 


_ _ The authors are to be congratulated on their careful and important research — 

- regarding the feasibility of using triangular hyetographs to approximate the 
geometry of local rainstorm hyetographs for small drainage systems. However, 


_ the authors neglect to discuss some aspects of their work which are relevant 


storms occur in southern Nevada area. 


5B, Nevada Test ‘Site 


0.35 0.30 | 


First, no consideration is given to storms whose hyetographs ‘yield 
- extlanatee of the parameters a and b. This situation may arise when a storm _ 
hyetograph does not have a rising leg; i.e., the maximum storm intensity occurs — 
at time ¢ = 0. Using the authors’ "definitions, the writer has examined 345 _ 
‘storm records from a weighing bucket gage located at well 5B on the Nevada — 
_ test site. The time scale of the gage records was such that for storms of 2 
h or less in duration an accurate hyetograph could not be constructed. Thus, 
of the 345 storms considered, only 127 could be analyzed by the authors’ methods. _ 
Of these 127 storms, 18 resulted in negative estimates or either a or b and 
‘subsequently a° and b°. In Table 6 the results of this study y are summarized. e 
In this table, Case I denotes that storms resulting in negative estimates of either 
“June, 1980, by Ben Chie Yen and Ven Te Chow (Proc. Paper 15452). lp he Le 
Assoc. Research Prof., Desert Research Inst., Water Resources Center, Univ. of Nevada 
System, 4582 Maryland Parkway, Las Vegas, Nev. 


4 
_| 
7 
in months Case | Case Case Case | Case | Case Case 
j 0.30 0.19 | 0.86— 0.86 
23) 29) | @ | @ 
: 0.49 0.41 | 0.50 0.43 0.26 | 0.86 
0.54 | 0.57 | 0.49 | 0.49 | 0.61 | 061 | 0.37 | 0.75 
i ee Oe | 0.46 | 0.51 | 0.43 | 049 | 0.55 | 061 | 0.27 | 027 J 
The (64) (4) | (18) | 1S) | | 
I 0.50 | 0.51 | 052 | 059 | 048 | 0.47 | 0.50 | 0.43. 4 


number of storms ina a given 1 category. on 
_ Although the sample size is small, the results indicate that the effect of 
_ considering hyetographs which result in negative estimates of a° or b° can 

be significant; e.g., storms occurring in months 5-9 with durations greater than : 
4h but less than 10 h. It is further noted that hyetographs having this c characteristic — 
composed 14% of the total sample at this site. 

— Second, the authors evaluate the validity of their hypothesis by comparing 

the measured and computed hydrographs for several experiments. Although - 
method of verification has merit, it does not provide detailed information regarding © 

the validity of the triangular hyetograph hypothesis. In a number of cases at 


the well 5B site, the authors’ triangular hyetograph hypothesis yields a very 
om estimate of the measured hyetograph. The hyetograph in Fig. 10 cannot 


Hyetograph ~ 


Block Diagram 
Actuol Hystograph 


FIG. 10.—Comparison of ind for Storm of February 


be accurately approximated by the authors’ hypothess. - Although the use of 


additional moments might resolve this problem, the authors’ primary objective 7 
would be defeated. ‘This situation arises in the | case of large w winter 
- must be given to the size of the storm relative to the size. of the: drainage ~ 
~ basin. If a single storm approximately as large as the drainage basin is located ® 
i over the basin, then a very simple hyetograph might result but in the case 
of large storms, large relative to the basin, composed of many storm cells 
then additional factors must be considered. 6) ‘i. Se 
, s Acknowledgment.—The work described in this discussion was ‘sponsored by 
the U.S. of Energy wader DE-AC T6DPOI253. 


4 


| 
q 
| 
a 


Closure by Ben Chie Yen‘ and Ven Te Chow, Fellows, ASCE 
_ French’s thoughtful discussion is appreciated. Concerning the first of his two 
comments, namely on the negative ' values of the triangular hyetograph f parameters — Je 
aandb (Fig. 2), the effect was previously evaluated and reported in pp. 208-211 
_ Ref. 20. However, because of the length limitation of the paper, this as 
_ well as many other data problems was not elaborated in the paper. In fact, © 
; these problems, as they have not appeared in the literature (although many #: 
Fesearchers must have encountered and overcome them) deserve their own right 
The anomaly of negative values of a and b indeed occurs fo or some rainstorms, 
Bo when the time interval of the data is large, such as 2h used in 
Table 6. To avoid the anomaly, rainstorms having negative values of a or b 
“can be disregarded as suggested by French, resulting in a significant increase 
in the m mean value of a° and decrease in in the mean value of os Alternatively, 
one may arbitrarily set the values of a or b to zero for those rainstorms having ; 


negative a or b and include them in the statistics instead of i i them. 


DURATIONS, in hours 


of storms of ‘Storms storms 
(2) 


The result is a slightly inevenssi't mean values of a° and decreased mean value 
of b° from the case allowing negative values, as ‘shown in Case 3 for a’ * in 
- Table 7. These a’ values are computed from the June-August data given in 
- Pa 6 by F French assuming that the three disregarded rainstorms all have negative 
he . Values of a° for Cases 1 and 2 given by French are also listed in Table — 
7 o 7 for comparison. The statistics given in the writers’ paper were obtained by — 
_ Considerably more rainstorm statistics have | been obtained, using 5-min and 
: hourly rainfall data, since the publication of the paper. It has been found that 
_ in order to obtain. . reasonably reliable values of the nondimensional hyetograph j 
_ parameters, the sample size must not be small; normally there should be more 
than 100 rainstorms in each class. As the sample size becomes large, the effect 
vi of not allowing negative a and b becomes relatively small and insignificant, 
& as shown in Table 8 for the data of 469 rainstorms for the months of June-August, 


from 1959-1972, at Urbana, Ill. (20). 


“Prof. of Civ. Engrg., Univ. of Illinois at Urbana-Champaign, Urbana, Ill. 61801. = 
Prof. of Civ. and Hydrosystems Engre., Univ. of at Urbene- Urbana, 


Case 


to that of a : particular rainstorm, one e should realize that the triangular hyetograph > J 
based on the statistical mean values is an estimate of the expected hyetograph © 
in the future. It would not fit perfectly with any rainstorm of the past or future; = 
this is analogous to the case that no American family actually has 2.2 children 
although this is the statistically expected number of children in an American — 
family. As no two hyetographs are identical, obviously a synthetic hyetograph 

- that would precisely fit one actual hyetograph would not fit any other hyetograph. g 
The purpose of the synthetic hyetographs is to give a reasonable estimate of 

- the future that has not yet happened and could be in any of the infinite umber 

ol of shapes. Recently, particularly in the field of urban hydrology, too much 
TABLE 8.—Effect of Negative of a and b, Urbana, Ill., June-August, 1959-1972 

or 


3.4620 | 1.00. 


146 0.43. 


attention has been ; given to trying to duplicate what has happened in the past — 


at the expenses « fora for the unknown future. 


eel Panapoxes IN THE His History OF Hyprauics® 
_ Following his kind words about the paper, my sole discusser remarked that 
I now had the “right and privilege” i. e., obligation) to keep certain additional 
_ As recounted in my bicentennial story of American hydraulics (3), Shields’ 7 
- work on sediment movement was first brought to the profession’s attention 


the several of my early pee 4 on n the subject, but the 1 notion a that i it was = : 


in Engineering Hydraulics t to which Professor Gill refers (1) specifically en 
_ Du Boys’ original contribution. Although Shields inserted his own numerical 
values into Du Boys’ transportation equation, Shields’ real accomplishment was _ 
“June, 1980, by Hunter Rouse (Proc. Paper 15475), 


-* Carver Prof. Emeritus, The Univ. of lowa, lowa City, lowa 52242; also Visiting Prof., 
Colorado State | Univ., Fort Collins, Colo. 80523. 
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_ that of correlating the initial Stages of movement with the grain Reynolds num number, _ 
as. as Gill correctly concluded. I am not aware of any paradox in this regard. 
Until the discussion appeared, moreover, I had never seen Brown’ Ss name cor 
“mine mentioned in connection with the Einstein and Kalinske formulations. _ 


‘While bringing several chapters of Hydraulics to 


: considerable last- -moment | revision, and another was ‘completely replaced. The 


late Carl Brown, on the other hand, had consented to serve as author of a 


sediment chapter only if he received help with a few of the sections. This, 
_ several of us undertook * provide (albeit more extensively than I had anticipated), 
as acknowledged at the end of the chapter. In the bed-load section, I made 
ina comparable 1 manner; the | underlying | analyses, however, remained their own, - 
and association of either Brown’s name or mine with theirs is hardly justified. 


I thus agree fully with the discusser’s implication that there are already enough 


a special effort to Present | the Shields, Einstein, and Kalinske relationships 


Rouse, H., Hydraulics in the United Ste States, tes, 1776-1976, Yo Iowa Institute of Hydraulic | 


Estimate oF Point RAINFALL 


laos Discussion by Ismael Pagan-Trinidad,” A. M. ASCE 
and B Ben C. Yen,“ F. ASCE 
>a The authors presented an interesting study on handling missing point rainfall 
- record. It has long been known that the inverse-distance estimate of point rainfall 
is superior to the standard National Weather Service methods (5). In fact, the ; 
_ inverse-distance method has been taught in the undergraduate hydrosystems 
engineering course at the University of Illinois at Urbana-Champaign, Ill., since 
_ 1975. The authors’ result, which is based on one of the possitte inverse-distance — 
methods, further supports the advantage of this method. 
_ The writers have conducted an extensive study of spatial and temporal rainfall 
variabilities for engineering applications. The evaluation of different weighting © 
“methods to estimate missing precipitation records is a part: of the research. 
a July, 1980, by John R. Simanton and Herbert B. Osborn (Proc. Paper 15516). otal 
nil a Prof. Dept. of Civ. Engrg., Mayagiiez Campus, Univ. of Puerto Rico, Mayagiiez, 


@ - Puerto Rico; currently Grad. = Dept. of Civ. Engrg., Univ. of Illinois at Urbana- 


toa Prof., Dept. of Civ. Engrg., Univ. Of Hinois at at Urbana- -Champaign, Urbana, Ill. 61801. 


but also weighted subscribing- area and other An example 
of another inverse-distance technique is 

in which = the mean of D’. Since the complete results will appear in 


ay paper, suffice it to mention here that the selection of the preference technique 

- depends on a number of factors, including the « error measurement criteria, the 
type of rainfall record (i.e., storm event, daily, monthly, etc.), the geometry 
_ of the raingage network arrangement, and the physical constraints. The writers — 

have used over 300 rainfall events from a 12-yr record data collected at the : 
East Central Illinois | Raingage Network by the Illinois State Water Survey. The 
writers’ “results show | that the inverse- distance technique used by the authors _ 
is indeed a superior method under similar conditions, namely, single heavy 
rainstorm events and hydrometeorologically relative homogeneous areas. They 
also show that Eq. 3 does not apply for large values of b. Nevertheless, it 
should be pointed out that in all cases the isohyetal method gives equal or 
better results than any of the weighting methods evaluated, because the rainfall 
directional gradient effect can be more realistically represented. However, the 

7 es method is also more complicated and time-consuming than any of _ 

weighting methods, 

As mentioned a a a on the goodness of a method depends 

on the criteria of error measurement. The correlation Coefficient is is not necessarily — 


dependence between the t two . variates. . Moreover, the correlation coefficient — 
is a function not only of the distance between rain gages as shown in Fig. 7 
‘he but also of the orientation, synoptic pattern, type of event, and duration: 


pe used to compute the correlation coefficient. Also, the \ writers assume that 


Eq. 1, as printed, is in error, but the results. presented in Tables | and 2 and 


5. Gilman, C. S., “Rainfall,” Handbook of Applied Hydrology, Chow, ed., 
_McGraw- Hill Book Co., Inc., New York, N.Y., 1964, p. 9.28. 
% Huff, F. A., “Spatial and Temporal Correlation of Precipitation in Illinois,’ 


| 
| 
| 
| 
| | 
4 
| 
4 


DISCUSSION 


Force on Smt oF Forcep Jump 
ssion by I Karam S. Karki® and Shantanu K. Mishra* 


‘The au authors are to be commended for presenting some valuable data regarding 
the force experienced by ; a sill under various forced jump p conditions. The writers" 
would like to make a few comments and add some interesting results on sail 
_ The writers would like to know the authors’ observation about the deflection 
of the supercritical jet and the surface profile near the sill under flow condition 
I. When the downstream conditions do not affect the upstream flow, the writers _ 
have observed that the jet partly y splashes- -over and partly over turns to form 
an intense roller very close to the sill under the conditions 4 = (h/y,) = 1.25 _ 
and F > 3.5. It has been also reported that the authors could not identify 
- flow condition I for (h/y,) = 3. The writers feel that it would have been 
easier to understand this identification if it was explained in terms of surface — 
me 
di The coefficient of drag wader | flow condition I has been reported t to be constant — 
with a value approximately equal to 0.45. In this connection, the writers do 
not agree with the authors’ remark, ‘‘Although this ‘nonsplashing-over’ behavior © 
-°} does not exactly describe the flow condition I of the present investigations, — 
the Cd value and the fact that it is a constant { are consistent.’’ The writers 


“are o of the opinion that the ‘ ‘nonsplashing- over” condition for 4 = (h/y,) * 
1.25 and F > 3.5, as described in the first writer’s paper (4) is . the same as 
_ that of flow condition I. At the same time, it is not very correct to designate 
the drag coefficient under this condition of flow as the maximum drag coefficient, 
7 ha of forced jump, since there is no jump at all. Under no jump condition, 
vig 
=a the maximum drag experienced by the sill occurs under the ‘ ‘splashing over’ 
_ condition (4); the writers would prefer to call it the drag “coefficient under 
; flow condition I. Under this condition of flow, a knowledge of drag experienced — 
a - by the sill is not sufficient. It has been observed by the writers that substantial 
¢ negative pressures do occur on the downstream face of the sill. It is needless — 
to mention that these negative pressures: are of great importance in designing 
_ The authors’ second conclusion states that ‘the maximum drag coefficient — 
_ from 0.45-0.3 depending on the relative sill height,’ whereas their Fig. 
6 shows that this variation eopents on the Froude number and not on the 
i. —_ writers tried to predict pd drag coefficient of the sill at transition om 
i. teh, 1980, by Rangaswami Nerayenda and Loizois S. Schizas (Proc. Paper 15552). 


“Reader, Dept. of Civ. Engrg., Inst. of Tech., Banaras Hindu Univ., Varanasi, 221005, 
a “Reader, Dept. _ of Civ. Engrg., , Inst. of Tech., Banaras Hindu Univ., Varanasi, 221005, 


condition ition 1 t to tere momentum principle. The 
a obtained « on n the basis of the data Shown i in ‘Fig. 3 do not compare wa a 


(Gnnaie lines, distinguishing flow condition I from flow condition II are - 
“not as has also been reported by the authors. 


may be seen as the of flow II with =0, although 
this flow transition is difficult to establish but may be PO 
transition—the very first beginning of flow The 


wie t 6. 


be The values of y,/y, used for calculations are obtained by carefully extending 
back the various experimental curves of Fig. 4 for different Froude numbers 
up to the y-axis, i.e., where (X,/y,) = 0. Table 1 shows the results of these 
calculations. From this table, it is seen that the drag coefficients of the sill 
at the transition remains approximately constant with an average value of 0.28. - 
This compares well with the experimentally-observed value of 0.3. | 
_ Fiaally, the writers would like to point out that for almost all the observations , 
_ of flow condition II, the depth ratio, y, /y ,, is always higher than the correspond- 

4 ing sequence depth ratio of a free jump. This observation is contrary to all 
previous observations (2, 5, 12, 15, ‘16, 17, and the 


> 


1S. “Energy Dissipators for — and Outlet ose ’ by the Task Force on Energy 


| 
the ata of flow conditio ss mi i throu e e eauatio 
| 


DISCUSSION 
Dissipators for Spillways and Outlet Works Committee on Hydraulic neuabieds, Rn: 
__H. Berryhill, Chmn., Journal of the Hydraulics Division, ASCE, Vol. 90, No. HY1, 
Proc. Paper 3762, Jan., 1964, pp. 121-147, ooms 
16. Pillai, N. N., and Unny, T. E., “Shapes for Appurtenances in Stilling Basins,” * Journal — 
of the Hydraulics Division, ASCE, Vol. 90, No. HY3, Proc. Paper 3888, May, 1964, 
17. Weide, oad Effects of the Size and Spacing of Floor Blocks in the Control | 
in of the Hydraulic pom scat thesis presented to Colorado A. and M., at Colorado, in 


URBULENCE MEasu REMENTS IN Flow” 


_ The authors have succeeded in measuring Reynolds shear stress profiles in 
- a accelerating and decelerating flow in a flume as a simulation of a tidal flow. 
They got around the contamination problems of hot- film anemometers s by using q 
a small electromagnetic current ‘meter. 


with respect to the mean velocity. The data, however, do not seem to be entirely | 
consistent as regards these aspects. 
_ In Fig. 8 the hysteresis effect of the Reynolds stresses is significant, especially - : 


‘in the ‘more strongly accelerated and decelerated flows. If accelerated flows 


and decelerated flows are considered separately, this hysteresis effect must 

be perceptible in the difference between them. In fact, the data for T = 550 — 
“sec in Fig. 8 are based on such separate measurements. Rearranging the data 
of Fig. 6 by plotting the nondimensional Reynolds stress against 0/0, instead a 


(if apparent at all) than in . Fig. 8, although for T= 550 sec the measurements 
are the same or analogous in both cases. For T= 200 sec the time between 
the accelerating and the decelerating phase is too short in the case of Fig. 
8 to create a steady state in between. There is no reason, however, to expect — 
; an important influence of this absence of a steady state on account of Fig. 
- After all, the Reynolds stresses at the end of the accelerating phase in 
12 are comparable to the stresses at the beginning. of the decelerating 
phase or, because of the steady state in between, to the Reynolds stresses — 
in the steady state. So, the decelerating phase should also behave in accordance 
with Fig. 12, and little or none hysteresis should be measured. > Ban 
_ The data in Fig. 6 demonstrate, in conformity with some measurements in 
a tidal flow by Bowden, « et al. (20), that the location of the maximum Reynolds _ i” 
* “August, 1980, by Habib O. Anwar and Roy Atkins (Proc. Paper 15609). a oe 


*Sr. Scientific Officer, Fluid Mech. Group, Dept. of Civ. Engrg., Delft Univ. of Tech., 
Stevinweg 4, Delft, The Netherlands. 


Seti. 
= 
° 
difference between the magnitude Of the Keynoids shear stresses and their proriles 7 
: - accelerating and decelerating flow, and a hysteresis effect of these stresses _ 
— 


e.g., in flow. This decreasing Reynolds s stress near the bed does not 

_ agree with the bed shear stress obtained from Preston tube measurements. The ~~ 

data in Fig. 6, replotted in Fig. 12 na ere of the strong hysteresis 
effect i in Fig. 8, a (small) reversed effect at aty/D= 


eo Deduced fror from n Fig. 6: : (a) T = 550 sec; (b) T = 200 sec aval 


-06) occ = 000196 m/s’ 550s 
FIG. 13.—Variation of Surface Slope, i, with Dimensionless Velocity 0/0, 7 
14.—Quadratic Versus Linear Variation of Surface Slope with 
Velocity, 0/0,, tor T = 550 sec (Acceleration) “Ue 
As stronger decelerations are connected with less favorable “slopes, the 
observation that the slope for T= 200 sec remains favorable through the 
measurements, 5, does not seem correct. This is probably caused by | a possible 
; sign error in Fig. 10. After this correction the parts of Fig. 10 may be combined 


4 Contrary to the sixth conclusion in the original paper, the relation between 


« 
Q/D) 
anit 


slope i wal time or 0/0, in Fig. 10 may better by a quadratic. 
_ than by a linear relation, as shown in Fig. 14. As quadratic friction is likely — 
to dominate, this would be more acceptable (see also Eq. 5). olase 
re In Fig. 5 data by Gordon, et al. (7), found from measurements in a tidal — 
‘Cees are plotted. However, the maximum accelerations in this estuary were 
about a factor 100 smaller than 0.37 cm/s’, which means that the results of © 
the flume and those of the prototype cannot be directly compared. Therefore, 
the eighth conclusion of the original paper seems somewhat premature. 
20. Bowden, K. F., Fairbairn, L. A., and Hughes, P., ‘‘The Distribution of Shearing © 


_ Stresses in a Tidal Current,’’ Geophysical Journal of the Royal Astronomical Society, 


the 
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“Discussion H. Hawkins “M. ASCE 


a defensible procedure for determining curve numbers from data. His method © 
Satisfies the frequency- related daca in applying the curve number method 
. to peak flow calculations. 
_ An alternate acceptable nines also exists, i.e., that of Rallison and Cronshey 
- (10), which has the relative advantage of not requiring a lengthy record. Its 
application has been demonstrated by Springer, et al. (11). A useful investigation 
_ would be a comparison of the two methods. Also, the “‘visual fit’’ to the curve 


~ number frequency data tempts the suggestion of a statistically-correct fitting: 
F or a data analysis routine. Such could perhaps require satisfying a least- “squares — 

criterion for the runoff data estimated by the moment fit rainfall frequency i 

_ data. That is, what CN (or S) creates the best fit line to the observed runoff — 


frequency data using Eq. 2 with the rainfall frequency fit line estimates as_ 


4 Much of the power of the CN method lies in its adjustment to | 
_ Thus, it is of interest to consider the land condition of the five studied basins. 
= particular, was it consistent aan the sampling period, or in systems 


‘studied than smaller basins, whee a single act of random or designed land 
_ (urbanization, fire, or agriculture) could disrupt the entire anes regimen. 


“a “September, 1980, by Allen T. Hjelmfelt, Jr. (Proc. Paper 15693). 


= Dept. of Civ. and Environmental Engrg., Utah State Univ., Logan, Utah $4322. ities: a 


f 
| 
Hjelmfelt has made a useful contribution OY 
| 
a 
J 


and for the Sonoita. Creck data, a ac a of about 0. 05 seems to 0 fit.. This 
relationship, an alternate form of the deeply-entrenched ‘‘rational equation” 
has been found by the writer for other western watersheds (9), though its 
- emergence for this particular basin is surprising. Under these conditions, a 
_ relatively small consistent source area may be hypothesized, a fraction, C (from 
Eq. 4), of the watershed area. As shown in Ref. 9, this causes the observed 
curve number (as determined from storm P and Q observations) to decrease 
_ In the case of Sonoita Creek, assuming CN = 74 to be a statistically-correct 
fitting and scaling P values from the fit line in Fig. 5, the — of errors 


Probability, Error, 
as @ percentage 


enta as a percentage 


ea to” contemplate decisions made in other situations where overestimating of 
- runoff volume is not on the side of safety. An example of the latter might | 
be flood insurance studies, where a true best estimate of reality would be needed. 


9. Hawkins R. H. , “Runoff Curve Numbers from Partial Area Watersheds,” Journal 
Rallison, R. and discussion of ‘Runoff Curve Numbers with 

* by Richard H. Hawkins, Journal of the Irrigation and Drainage - : 

Division, ASCE, Vol. "105, No. IR4, Proc. Paper 15010, Dec., 1979, PP. 439-441. a 
‘ae E. P., McGurk, B. J., Hawkins, R. H., and Coltharp, G. B., ‘‘Curve Numbers 


: ‘from Watershed Data,’’ Proceedings of the Symposium on Watershed Management, 


‘The Sonoma Teek results clearly do not affirm the curve r _ 
| 
From the design safety standpoint, as with flood accommodating structures, _ 
_ it is comforting to have the CN method overestimate for extreme large events. — 
2.—Rainfall and Runoff for Sonoita Creek 

| 

| 
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ENTRAINED Air tn LineaRLy ACCELERATED Water Fiow" 


Discussion by Maurice J. Kenn? 


_ Abdin Salih will doubtless be pleased to mew that his important studies _ 
- at imperial College h have now been extended to include work on | entrained air * 
in decelerating water flows epee constant, linear, and nonlinear 
be severely shattered by the host stream. 
_ This more recent knowledge has been effectively applied at Bolarque Dam, 
in Spain, where small quantities of air have been successfully injected S| 
the decelerating water stream of the water-turbine draught tubes, in order to 
_ suppress cavitation, cavitating vortices, associated water hammer (at a frequency | 
of 3.9 Hz) with draught-tube vibrations, and the resonating vibrations of the 
-m high concrete dam, with first-mode rocking frequency also 


‘2B. Stafanakos, J. P., “‘Entrained Air in Decelerating Water Flows,”’ thesis 
to the University of London, at London, England, in 1979, in partial fulfillment of 
the requirements for the degree of Doctor of Philosophy. 
14. Kenn, M. J., Cassell, A. C., and Grootenhuis, P., “Vibrations at the Bolarque Dam,’ 
Symposium on Practical Experiences with Flow-Induced Vibrations, International 
_ Association for Hydraulic Research/International Union of Theoretical and Applied 


“October, 1980, by Abdin M. A. Salih (Proc. Paper 15728), 
a ? Acting Head, Hydr. Section, Civ. Engrg. Dept., Imperial College, London, SW7 2BU, 
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TRANSIENTS Fottowin \ VALVE Crosure 


sl Discusion Alan E. ‘Vardy* 
_ The finite element technique has justifiably gained wide acceptance in the _ 
engineering world. However, considerable difficulty has been experienced 
adapting it to the conditions pertaining in hyperbolic problems. It is to be hoped — 


work will continue towards this end. 
{ —- support of their conclusion that ‘‘both the method of characteristics and — 


the finite element method are shown to give an accurate simulation of pressure __ 

"surges produced through rapid valve closure,”’ the authors « draw specific attention 
- Fig. 5 of their paper. For reasons given in the following, the writer feels 
uneasy with this interpretation oftheir results. 

1. ‘The | term “rapid” is the valve closure time of about 


act activity normally “associated with water hammer | is not strongly it in 
_ 2. When the a:thors presented these results at a recent conference (9), several 
ry discussers drew attention to the marked differences between the magnitudes 
» a of the 1 measured and ‘Predicted pressures w wane b> > > 0. 9s sec. Significantly better = 
suspects that the descrepancy is due more to the simulation of the valve as 
a boundary condition than to the simulation of the flows in the pipe. The likelihood — 

_ of this possibility is strengthened by the existence of significant discrepancies _ 
between theory and experiment even during the preceding period of slowly 
. varying flow. As the discussion to the “earlier | paper has not been published, — 
‘the writer would like to ask again whether the authors have considered using 
the measured pressure history just upstream of the valve as the downstream — 

boundary condition for the theoretical analyses. Although a simulation of real ¥ 
— -valves is clearly advantageous in a practical environment, the elimination of - 
+ this source of error would be highly desirable for the ‘Purpose of verifying — 


the authors’ theoretical technique. 

Z The results pr presented in n Fig. | 5 can 1 be ‘interpreted | from another viewpoint 
“which encourages me to believe that the authors’ finite element technique is 
_ much better than the disappointing comparison which the experiment suggests. _ 
There is good reason to suppose that the method of characteristics (MOC) 

7 ‘is capable of predicting water hammer behavior very closely. Therefore, the 

validity of any new theoretical approach can be evaluated relatively cheaply 
_by comparing its predictions with those of the MOC. Large differences indicate - 
errors in the new approach; small differences represent open questions that — 


a “October, 1980, by C. S. Watt, J. M. Hobbs, and A. P. Boldy (Proc. Paper 15740). it 
; “Prof., Dept. of Civ. Engrg., Univ. of Dundee, Dundee DDI 4HN Scotland. Cree ac 


| 
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DISCUSSION 


must be answe were ed by precise writer ott take issue 
the authors on some of the limitations that they unreasonably attribute to the 
4 A MOC, but this would be quibbling. The point is that the agreement shown 
__ between the MOC and the authors’ finite element approach is quite close. Their 

- analysis is therefore proceeding along the right li lines and is well worth extending, 


Watt, C. S., Boldy, A. P., and Hobbs, J. M., ‘Combination of Finite Difference 
and Finite Element Techniques in Hydraulic Transient Problems,”’ Proceedings, Third 
International Conference on Pressure Surges, Mar., 1980. 


Vexociry oF SHELLS AS COASTAL SEDIMENT“ 


: - In the authors’ summary an important en is that ‘ ‘shells in stable 
eer of settling have comparatively lower drag coefficients than those settling _ 


in the unstable mode’’; this is clearly evident if Figs. § 8 and 9 are are compared. wart 


4 To judge the significance of this difference, it would be important to have ; 
a clear definition of the fall velocity, W, since it enters into the definition _ 
Of Cp; see Eq. 7. Was the fall velocity, W, obtained with the true path length © 
or the ** projected” path length; how was it measured and with what accuracy? 
Furthermore, an attempt should be made to quantify the difference between | 
= and unstable modes of settling. In this way it might be helpful to explore 
the use of a dimensionless ‘“‘stability criterion,’’ such as a Strouhal number. 7 


This might turn out to be a difficult task, if not the entire particle’s path was ; 


- While it is certainly appreciated that the shell geometries are compared with — 
shape factors, it appears that something could be le learned from a comparison Z 
with freely-falling disk data. A rather valuable summary of data for disks is 
given in Ref. 14; interestingly enough the data of the unstable settling mode, , 
as given in Fig. 9, show the same trend as freely-falling disks. In this content, 
it would be also worthwhile to regard the concept of * Wry. frequency”’ 
as put forward by Marchillon, et al. (13). “aay 
13. Marchillon, E. K., and Gauvin, W. H., “O ‘Oscillatory Motion of Freely 


_ “November, 1980, by Ashish J. Mehta, Jieh Lee, and Bent A. Christensen (Proc. Paper 


| 


os multiplier but more appropriate for quantifying closed-conduit flow, is Q /V gD° gD°, 


tsi Falling Disks, ” Physics of Fluids, Vol. 7, No. t 1, Dec., 1964. a wave 
14. Schiller, L., “‘Fallversuche mit Kugeln und Scheiben,’’ Handbuch der Experimen- 


talphysik, Vol. IV /2, » 5 Germany, 1932, pp. 376-379. 


Scour Arou ND Brince Piers at Hicu Fiow V VEL 


Discussion by Fred W. Blaisdell,’ F. ASCE 


oe This paper is exceptionally well-written, easy to read and comprehend, and “te 
complete within the scope of the research reported. The writers thank the authors - 
of this paper for their efforts on behalf of the reader. — ent a 
Since 1966, the writers have been studying local scour that has several points 
of similarity to the work of the authors. The writers’ research is on scour 
_ produced in uniform-sized sand beds by the clear-water efflux from a cantilevered 
pipe. The pipe invert is positioned from 2D below to 8D above the tailwater 
surface (D = the pipe diameter). The surface of the bed is 0.5D above the 
- tailwater level. A channel is formed in the bed downstream of the pipe exit. 
_ The discharge parameter, differing from the Froude number only by a constant 


in which QG= the discharge; and g = the acceleration due to | gravity. The 
* tested is from 0.5-5. Unlike the work of Jain and Fischer, a 
here there is no flow over the bed to transport bed material from upstream 
into the scour hole. All tests are three-dimensional, and most of them involve 
three-phase flow—air, water, and sand. Time is an additional variable. The 
— used, described after Jain and Fischer, | have the 
"Three groups of tests were made. odd sauces 
_ For the first group of tests, the scoured material deposited in and completely 
s the channel. The scoured material deposited because the channel was” 
stable and the velocities were insufficient to transport the scoured material. 
i: a result, this procedure was abandoned and a test method devised for keeping 


downstream channel clear of deposited sediment. 


4 “November, 1980, by Subhash C. Jain and Edward E. Fischer (Proc. Paper 15845). ong 


* Research Hydr. Engr., United States Dept. of Agr., Sci. and Education Administration, 
7 Agricultural Research, St. Anthony Falls Hydr. Lab., Third Avenue Southeast at Mississippi 
—.: Hydr. Engr., United States Dept. of Agr., Sci. and Education Administration, Agricul-— 
13 tural Research, St. Anthony Falls Hydr. Lab., Third Avenue Southeast at Mississippi 
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For test group, a “‘suction dredge” in n the form ofa pipe cecillating a 
_ longitudinally along a short length of the excavated channel just downstream a 
_ of the scour hole kept the channel clean of the scoured material deposited 
the ‘channel. This procedure proved to be 


‘dae explanation of s scour at bridge | piers to the writers’ scour at t cantilevered 
pipe outlets—explains what the writers observed during the second group of — 


_ Experimental observations show that the side slopes of the scour hole © 


. . are approximately equal to the natural angle of repose of the bed ‘ 


material. This indicates that the scour mechanism is in the immediate _ 
vicinity of the [point where the plunging jet impinges on the sand bed] pup 
_ and that the sediment slides in toward the bass: before it is removed. ry party 
visual observation of the scour Process . . shows that a dynamic 
equilibrium exists between the scour hole and. _ . streamflow which 
is not apparent | when the flow is stopped; the flow field near the [point __ 
of jet impingement on the bed] is strong enough to support the sides fi 
* of the scour hole at angles greater than the angle of repose of the sediment. _ 
wall periodically collapses and dumps sediment into the hole . . 
.. as the fluid forces supporting it become unstable. The maximum ‘scour > 
depth under such unstable conditions cannot be measured 


stopping the flow. 


The first method used the data obtained from the second | group of tests. 

_ The logarithm of dimensionless rate of scour (dimensionless scour depth divided 
® dimensionless time) was plotted against the logarithm of dimensionless time. | : 
A rectangular hyperbola was fitted to the data as described in Ref. 4. The 4 

“ultimate scour depth was computed from the asymptote of the hyperbola. a 


d group of tests to provide data to determine > 
the maximum m depth of scour. In these tests the maximum scour depth 2 
determined by actual measurement as described in the following. But first 

a description of the scour at the end of the time-dependent tests is in order. ee 
At the conclusion of the time-dependent scour tests, sediment was st y? 
suspended by the jet turbulence . This is shown in the authors’ ’ Fig. 9. Just 
; - upstream of the authors’ pier there is a relatively wide dark line Botrold 
~ movement of bottom material toward the periphery of the scour hole. ae 
dark line curves upward and toward the pier and there is a cap of ae 
material directed toward the upstream angle-of-repose slope of the scour hole. — 
This plume of suspended material is typical of conditions at the end of the | 
writers’ time-dependent cantilevered outlet scour tests, which usually were a 
terminated at an elapsed time of 10,000 min—seven days. Energy from the 
jet supported this suspended material —energy that otherwise be available 


_ For the writers’ third group of tests, the “dredge” * suction was extended 


J ay 
{ 
| 
| 
— 
@ 
| Like Jain and Fischer, the writers were interested in determining the maximum _ 
| 
| 
7 


material made available the energy previously used to material to further 
"scour the bed. The sides became unstable and sluffed into the base of the | 
_ hole where the sediment was suspended and “‘dredged”’ out. This process greatly 
= all scour hole dimensions. Turbulent bursts continued to move material _ 
around intermittently at the base of the scour hole even at the termination 
of the third group of tests, but the side slopes extended along a smooth line 
to the bottom of the hole; the e steep face at the base of the — in the authors’ 
Fig. 9nolongerexisted. 
_ Cantilevered outlet scour hole data are still being obtained and analyzed at 
P this writing. Presentation of quantitative data at this time is not appropriate. 
This discussion is written to emphasize the similarity between the writers’ a 
- qualitative observations and those of the authors. These two experiences show 
that, to insure the safety of structures subjected to scour, the scour tests must _ 


: e Blaisdell, F. Anderson, C. L., and Hebaus, G. G. , “Ultimate Dimensions of 


Local Scour,”” Journal of the Hydraulics Division, ASCE, Vol. 107, No. HY3, Proc. Z 
‘Paper 16144, Mar., 1981, pp. 337. Ww 
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to which the paper should be referred. Those who are planning to submit material will expedite 
the review and publication ees by omg with the following | basic requirements: — 


1. ‘Titles must have a length not not exceeding 5 50 spaces. 
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